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Abstract 

During 2004, four divisions of the American Physical Society commissioned a 
study of neutrino physics to take stock of where the field is at the moment and 
where it is going in the near and far future. Several working groups looked at 
various aspects of this vast field. The summary was published as a main report 
entitled "The Neutrino Matrix" accompanied by short 50 page versions of the report 
of each working group. Theoretical research in this field has been quite extensive 
and touches many areas and the short 50 page report [1] provided only a brief 
summary and overview of few of the important points. The theory discussion group 
felt that it may be of value to the community to publish the entire study as a 
white paper and the result is the current article. After a brief overview of the 
present knowledge of neutrino masses and mixing and some popular ways to probe 
the new physics implied by recent data, the white paper summarizes what can 
be learned about physics beyond the Standard Model from the various proposed 
neutrino experiments. It also comments on the impact of the experiments on our 
understanding of the origin of the matter-antimatter asymmetry of the Universe 
and the basic nature of neutrino interactions as well as the existence of possible 
additional neutrinos. Extensive references to original literature are provided. 
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1 Introduction 



Our understanding of neutrinos has changed dramatically in the past six years. Thanks 
to many neutrino oscillation experiments involving solar, atmospheric, accelerator and 
reactor (anti)-neutrinos [2,3], we have learned that neutrinos produced in a well defined 
fiavor cigcnstatc can be detected, after propagating a macroscopic distance, as a different 
flavor eigcnstate. The simplest interpretation of this phenomenon is that, like all charged 
fermions, the neutrinos have mass and that, similar to quarks, the neutrino weak, or 
flavor, eigenstates are different from neutrino mass eigenstates i.e., neutrinos mix. This 
new state of affairs has also raised many other issues [4] which did not exist for massless 
neutrinos: For example, (i) massive Dirac neutrinos, like charged Icptons and quarks, can 
have nonzero magnetic dipole moments and massive Dirac and Major ana neutrinos can 
have nonzero transition dipole moments; (ii) the heavier neutrinos decay into lighter ones, 
like charged leptons and quarks, and (iii) (most importantly) the neutrinos can be either 
Majorana or Dirac fermions (see later for details). 

Learning about all these possibilities can not only bring our knowledge of neutrinos 
to the same level as that of charged leptons and quarks, but may also lead to a plethora 
of laboratory as well as astrophysical and cosmological consequences with far-reaching 
implications. Most importantly, knowing neutrino properties in detail may also play a 
crucial role in clarifying the blueprint of new physical laws beyond those embodied in the 
Standard Model. 

One may also consider the possibility that there could be new neutrino species beyond 
the three known ones (z/g, ^'/i, i^r)- In addition to being a question whose answer would 
be a revolutionary milestone pointing to unexpected new physics, it may also become 
a necessity if the LSND results arc confirmed by the MiniBooNE experiment, now in 
progress at Fermilab. This would, undoubtedly, be a second revolution in our thinking 
about neutrinos and the nature of unification. 

The existence of neutrino masses qualifies as the first evidence of new physics beyond 
the Standard Model. The answers to the neutrino-questions mentioned above will add 
substantially to our knowledge about the precise nature of this new physics, and in turn 
about the nature of new forces beyond the Standard Model. They also have the poten- 
tial to unravel some of the deepest and most long-standing mysteries of cosmology and 
astrophysics, such as the origin of matter, the origin of the heavy elements, and, perhaps, 
even the nature of dark energy. 

Active endeavors are under way to launch the era of precision neutrino measurement 
science, that will surely broaden the horizon of our knowledge about neutrinos. We 
undertake this survey to pin down how different experimental results expected in the 
coming decades can elucidate the nature of neutrinos and our quest for new physics. 
In particular, we would like to know (i) the implications of neutrinos for such long- 
standing ideas as grand unification, supersymmetry, string theory, extra dimensions, etc; 
(ii) the implications of the possible existence of additional neutrino species for physics 
and cosmology, and (iii) whether neutrinos have anything to do with the origin of the 
observed matter-antimatter asymmetry in the universe and, if so, whether there is any 
way to determine this via low-energy experiments. Once the answers to these questions 
are at hand, we will have considerably narrowed the choices of new physics, providing a 



6 



giant leap in our understanding of the physical Universe. 

This review grew out of a year long study of the future of neutrino physics conducted 
by four divisions of the American Physical Society and is meant to be an overview of 
where we stand in neutrino physics today,^ where we are going in the next decades and the 
implications of this new knowledge for the nature of new physics and for the early universe. 
We apologize for surely missing vast parts of the neutrino literature in our references. We 
expect this overview to be supplemented by other excellent existing reviews of the subject 
in the literature. Regarding more references and the more experimental aspects of the 
topics under study, we refer to the other working group reports, the Solar and Atmospheric 
Experiments [5], the Reactor [6], the Neutrino Factory and Beta Beam Experiments and 
Development [7], the Neutrinoless Double Beta Decay and Direct Searches for Neutrino 
Mass [8] and the Neutrino Astrophysics and Cosmology [9] WGs. In particular, we have 
not discussed theoretical models for neutrino masses except giving a broad outline of 
ideas and getting beyond it only when there is a need to make some phenomenological 
point. Nonetheless, we hope to have captured in this study the essential issues in neutrino 
physics that will be relevant as we proceed to the next level in our exploration of this 
fascinating field. 

1.1 Our present knowledge about masses and mixings 
1.1.1 Dirac versus Majorana Neutrinos 

The fact that the neutrino has no electric charge endows it with certain properties not 
shared by the charged fermions of the Standard Model. One can write two kinds of 
Lorentz invariant mass terms for the neutrino, Dirac and Majorana masses, whereas for 
the charged fermions, conservation of electric charge allows only Dirac-type mass terms. In 
the four component notation for describing fermions, commonly used for writing the Dirac 
equation for the electron, the Dirac mass has the form iJjiIj, connecting fields of opposite 
chirality, whereas the Majorana mass is of the form ■ip'^C^^ip connecting fields of the same 
chirality, where is the four component spinor and C is the charge conjugation matrix. 
In the first case, the fermion ip is different from its antiparticle, whereas in the latter case 
it is its own antiparticle. A Majorana neutrino implies a whole new class of experimental 
signatures, the most prominent among them being the process of neutrinoless double beta 
decay of heavy nuclei, {f3j3ou)- Since f3f3ou arises due to the presence of neutrino Majorana 
masses, a measurement of its rate can provide very precise information about neutrino 
masses and mixing, provided (i) one can satisfactorily eliminate other contributions to 
this process that may arise from other interactions in a full beyond-the-standard-model 
theory, as we discuss below, (ii) one can precisely estimate the values of the nuclear matrix 
elements associated with the f3f3ou in question. 

The expressions for the Dirac and Majorana mass terms make it clear that a theory 
forbids Majorana masses for a fermion only if there is an additional global symmetry under 
which it has nonzero charge. As noted above, for charged fermions such as the electron 
and the muon, Majorana mass-terms are forbidden by the fact that they have nonzero 

^The bulk of this report was finalized at the end of January 2005. We have not included the (sometimes 
substantial) progress that has been obtained in several areas of neutrino physics since then. 
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electric charge and the theory has electromagnetic U{1) invariance. Hence all charged 
fermions are Dirac fermions. On the other hand, a Lagrangian with both Majorana and 
Dirac masses describes, necessarily, a pair of Majorana fermions, irrespective of how small 
the Majorana mass term is (although it may prove very difficult to address whether the 
fermion is of the Dirac or the Majorana type when the Majorana mass-term is significantly 
smaller than the Dirac mass term). Hence, since the neutrino has no electric charge, 
the "simplest" theories predict that the neutrino is a Majorana fermion meaning that a 
Majorana neutrino is more natural (or at least requires fewer assumptions) than a Dirac 
neutrino. In most of the discussions below we assume that the neutrino is a Majorana 
fermion, unless otherwise noted. 

1.1.2 Neutrino mixings 

We will use a notation where the electroweak-doublet neutrino eigenstate (defined as the 
neutrino that is produced in a charged-current weak interaction process associated with 
a well-defined charged lepton) is denoted by z/q,, with a = e, /i, r. We will also consider z/^ 
to include a set of possible electroweak-singlet ( "sterile" ) neutrinos. Corresponding to 
these 3 + Us neutrino interaction eigenstates are 3 -|- mass eigenstates of neutrinos, Ui. 
We will order the basis of mass eigenstates so that < m| and Am^2 < l^^isl) where 
Amfj = nij — mf. The neutrino interaction eigenstates are expressed in terms of the 
mass eigenstates as follows: = Uaii^i, where U is a {3 + Ug) x (3 -|- n^) dimensional 
unitary matrix. For the active neutrinos, with a = e, /x, r, the relevant submatrix is thus 
a rectangular matrix with three rows and 3-1-^^ columns. In seesaw models, the entries in 
the columns 4, ...3-\-ns are very small, of order mD/iTiR, where ttid is a typical Dirac mass 
and is a large mass of a right-handed Majorana neutrino. Motivated by these models, 
one commonly assumes a decoupling, so that to good approximation the electroweak- 
doublet neutrinos can be expressed as hnear combinations of just three mass eigenstates, 
and hence one deals with a 3 x 3 truncation of the full {3 + ns) x {3 + ns) neutrino mixing 
matrix. Since only the three electroweak-doublet neutrinos couple to the W, the actual 
observed lepton mixing matrix that appears in the charged weak current involves the 
product of the 3 x (3 -|- n^) rectangular submatrix of the full lepton mixing matrix with 
the adjoint of the 3x3 unitary transformation mapping the mass to weak eigenstates 
of the charged leptons. Thus, the lepton mixing matrix occurring in the charge-lowering 
weak current has three rows and 3+ns columns, corresponding to the fact that, in general, 
a charged lepton a couples to a z/q, which is a linear combination of 3+ns mass eigenstates. 
Henceforth, unless explicitly indicated, we shall assume the above-mentioned decoupling, 
so that the neutrino mixing matrix is 3 x 3, and will use U to refer to the observed lepton 
mixing matrix, incorporating both the mixings in the neutrino and charged lepton sector. 
Neutrino oscillations and the mixing of two mass eigenstates of neutrinos, Ui and z/2, to 
form the weak eigenstates z/g and z/^ were first discussed by Pontecorvo and by Maki, 
Nakagawa, and Sakata [10]. The 3x3 truncation of the full neutrino mixing matrix is 
often called the MNS, MNSP, or PMNS matrix in honor of these pioneers. 

For the case of three Majorana neutrinos, the lepton mixing matrix U can be written 
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as VK, where V will be parametrized as 

(C12C13 S12C13 si3e"*''\ 

-S12C23 - Ci2S23Sl3e*'' C12C23 - Sl2S23Sl3e"^ S23C13 , (l) 
S12S23 - Ci2C23Sl3e*'' -C12S23 - Sl2C23Sl3e"^ C23C13 / 

while X = diag(l,e^'^Se*('^2+5)) [11^12]. 

Neutrino oscillation experiments have already provided measurements for the neutrino 
mass-squared differences, as well as the mixing angles. At the Scr level, the allowed ranges 
are [13] sin^ 2^23 > 0.87; 1.4xlO~3cV^ < lAm^g] < S.SxlQ-^eV^; 0.70 < sin2 2^i2 < 0.94; 
7.1 X 10"^ eV^ < Am?2 < 8-9 x 10"^ eV^; sin2^i3 < 0.051 [14]. There is currently no 
constraint on any of the CP odd phases or on the sign of Am^g. Note that in contrast to 
the quark sector we have two large angles (one possibly maximal) and one small (possibly 
zero) angle. 

1.1.3 Matter effect on neutrino propagation 

A very important fact about neutrinos that wc seem to have learned from solar neutrino 
data is that neutrino propagation in matter is substantially different from that in vacuum. 
This effect is known as the MSW (Mikheev-Smirnov-Wolfenstein) effect [15] and has been 
widely discussed in the literature [16] . There is however an important aspect of the favored 
large mixing angle (LMA) MSW solution which needs to be tested in future experiments. 
The LMA solution predicts a rise in the survival probability in the energy region of a few 
MeV as we move down from higher to lower solar neutrino energies. Since the present data 
do not cover this energy region, new solar neutrino data is needed in order to conclusively 
establish the LMA solution [17]. 

1.1.4 Neutrino masses 

Given the current precision of neutrino oscillation experiments and the fact that neutrino 
oscillations are only sensitive to mass-squared differences, three possible arrangements of 
the neutrino masses are allowed: 

(i) Normal hierarchy, i.e. mi < m2 <S rris. In this case Am23 = "^i ~ "^2 > 0? 

rriz — \/ Am23 — 0.03 — 0.07 eV. The solar neutrino oscillation involves the two 
lighter levels. The mass of the lightest neutrino is unconstrained. If rui <^ m2, then 
we find the value of 1712 — 0.008 eV. 

(ii) Inverted hierarchy, i.e. mi ~ m2 ^ [18] with mi,2 — \/Arn|^ ~ 0.03 — 0.07 eV. 
In this case, solar neutrino oscillation takes place between the heavier levels and we 
have Am^g = ml — m\ < 0. We have no information about 777.3 except that its value 
is much less than the other two masses. 

(iii) Degenerate neutrinos [19] i.e., mi ~ 7772 — ^773. 

The behaviors of masses for different mass patterns are shown in Fig. 1. 



9 




-3 -2.5 -2 -1.5 -1 -0.5 0.5 -3 -2.5 -2 -1.5 -1 -0.5 0.5 

Log (mMiN [eV] ) Log {mmN [eV] ) 



Figure 1 : The three hght neutrino masses as a function of the hghtest mass for the normal (left 
plot) and inverted (right plot) hierarchy. 



1.1.5 Overall scale for masses 

Oscillation experiments cannot tell us about the overall scale of masses. It is therefore 
important to explore to what extent the absolute values of the masses can be determined. 
While discussing the question of absolute masses, it is good to keep in mind that none of 
the methods discussed below can provide any information about the lightest neutrino mass 
in the cases of a normal or inverted mass-hierarchy. They are most useful for determining 
absolute masses in the case of degenerate neutrinos i.e., when all > 0.1 eV. 
Neutrino mass from beta decay 

One can directly search for the kinematic effect of nonzero neutrino masses in beta- 
decay by modifications of the Kurie plot. This search is sensitive to neutrino masses 
regardless of whether the neutrinos are Dirac or Majorana particles. These may be due to 
the emission, via mixing, of massive neutrinos that cause kinks in this plot. If the masses 
are small, then the effects will occur near to the end point of the electron energy spectrum 
and will be sensitive to the quantity mp = \/^i[Uei^^rn^ ■ The Mainz [20] and Troitsk [21] 
experiments place the present upper limit on m/3 < 2.3 eV and 2.2 eV, respectively. The 
proposed KATRIN [22] experiment is projected to be sensitive to m/3 > 0.2 eV, which will 
have important implications for the theory of neutrino masses. For instance, if the result 
is positive, it will imply a degenerate spectrum; on the other hand a negative result will 
be a very useful constraint. 

1.1.6 Neutrino masses and neutrinoless double beta decay 

Another sensitive probe for the absolute scale of the neutrino masses is the search for 
neutrinoless double beta decay, PPou, whose rate is potentially measurable if the neutrinos 
are Majorana fermions and niee = U^ifT^i is large enough [23,24], or if there are new 
lepton number violating interactions [25]. In the absence of new lepton number violating 
interactions, a positive sign of /5/?oi^ would allow one to measure mee- Either way, we 
would learn that the neutrinos are Majorana fermions [26]. However, if m^e is very small, 
and there are new lepton number violating interactions, neutrinoless double beta decay 
will measure the strength of the new interactions (such as doubly charged Higgs fields 
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or R-parity violating interactions) rather than neutrino mass. There are many examples 
of models where new interactions can lead to a P^qi, decay rate in the observable range 
without at the same time yielding a significant Majorana mass for the neutrinos. As a 
result, one must be careful in interpreting any nonzero signal in /3/3oi/ experiments and 
not jump to the conclusion that a direct measurement of neutrino mass has been made. 
The way to tell whether such a nonzero signal is due to neutrino masses or is a reflection 
of new interactions is to supplement /3/3oi/ decay results with collider searches for these 
new interactions. Thus coUider experiments, such as those at LHC, and double beta 
experiments play complementary roles. 

The present best upper bounds on (3(3ou decay lifetimes come from the Heidelberg- 
Moscow [27] and the IGEX [28] experiments and can be translated into an upper limit 
on niee ~ 0.9 eV [29]. There is a claim of discovery of neutrinoless double beta decay 
of enriched ^^Ge experiment by the Heidelberg- Moscow collaboration [30] . Interpreted in 
terms of a Majorana mass of the neutrino, this imphes niee between 0.12 eV to 0.90 eV. 
If conflrmed, this result is of fundamental signiflcancc. For a thorough discussions of this 
result (see also [31]) [32], we refer readers to the report of the double beta decay working 
group [8]. 

1.1.7 Cosmology and neutrino masses 

A very different way to get information on the absolute scale of neutrino masses is from 
the study of the cosmic microwave radiation spectrum as well as the study of the large 
scale structure in the universe. A qualitative way of understanding why this is the case 
is that if neutrinos are present in abundance in the universe at the epoch of structure 
formation and have a sizable mass the formation of structure is affected. For instance, 
for a given neutrino mass m, all structure on a scale smaller than a certain value given 
by the inverse of neutrino mass is washed away by neutrino free-streaming. This implies 
a reduced power on smaller scales. Thus, accurate measurements of the galaxy power 
spectrum for small scales can help constrain or determine neutrino masses. Recent results 
from the WMAP and surveys of large scale structure have set a limit on the sum of 
neutrino masses ^rrii < 0.7 — 2 eV [33,34]. More recent results from the Sloan Digital 
Sky Survey (SDSS) place the hmit of ^ < 1.6 eV. Hannestad [34] has emphasized that 
these upper limits can change if there are more neutrino species — e.g. for 5 neutrinos, 
^rrij < 2.12 eV if they are in equilibrium at the epoch of BBN. 

A point worth emphasizing is that the above result is valid for both Majorana and 
Dirac neutrinos as long as the "right-handed" neutrinos decouple sufficiently earlier than 
the BBN epoch and are not regenerated subsequently^. 

These limits already provide nontrivial information about neutrino masses: the limit 
X^j^i = 0.7 eV, if taken at face value, implies that each individual neutrino mass is 
smaller than 0.23 cV, which is similar to the projected sensitivity of the proposed KATRIN 
experiment. PLANCK satellite observations are expected to be sensitive to even smaller 

^In the Dirac case the "right-handed" degrees of freedom are decoupled because of the smallness of 
the corresponding Yukawa couplings. However, for very small temperatures, i.e. long after BBN, it is no 
longer appropriate to describe neutrinos in terms of chiral states. This means that strictly speaking there 
is a regeneration, but this does not affect BBN (see, e.g., [35]). 
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values of 'Y^^mi, thereby providing a completely independent source of information on 
neutrino masses. These results may have implications for models of sterile neutrinos that 
attempt to explain the LSND results. 

1.1.8 CP violation 

It is clear from Eq. (0) that, for Majorana neutrinos, there are three CP-odd phases that 
characterize neutrino mixings [11,12], and our understanding of the leptonic sector will 
remain incomplete without knowledge of these [36,37]. There are two possible ways to 
explore CP phases: (i) one way is to perform long-baseline oscillation experiments and 
look for differences between neutrino and anti-neutrino survival probabilities [38]; (ii) 
another way is to use possible connections with cosmology. It has often been argued 
that neutrinoless double beta decay may also provide an alternative way to explore CP 
violation [39]. This is discussed in Sec. El 

In summary, the most important goals of the next phase of neutrino oscillation exper- 
iments are: 

(i) To determine the value of 6'i3 as precisely as possible; 

(ii) To determine the sign of Am^g, or the character of the neutrino mass hierarchy; 

(iii) To improve the accuracy of the measurement of the other angles and the mass- 
squared differences; 

(iv) To probe the existence of the three CP odd phases as best as possible. 

The discussion above assumes a minimal picture for massive neutrinos where the most 
general Majorana mass for three neutrinos has been added. While this may be the picture 
to the leading order, it is quite conceivable that there are other interesting subdominant 
effects that go beyond this. It is of utmost interest to determine to what extent one can 
constrain (or perhaps discover) these new nonstandard phenomena, since their absence 
up to a certain level (or, of course, their presence) will provide crucial insight into the 
detailed nature of the New Physics. 

1.1.9 Prospects for determining whether neutrinos are Majorana or Dirac 

As an example of what we can learn from future experiments, we focus on three experi- 
ments — searches for neutrinoless double beta decay (down to the level of 0.03 eV level), 
studies to determine the sign of Am^g = (m| — m^, and the KATRIN experiment, which 
is sensitive to the effects of a nonzero neutrino mass down to 0.2 eV in tritium beta decay. 
The interplay between the possible results of these three experiments is summarized in 
Tabled 

We see that extremely valuable information will follow from the results of these ex- 
periments. 

1.1.10 Sterile neutrinos 

A question of great importance in neutrino physics is the number of neutrino species. 
Measurement of the invisible Z-width in LEP-SLC experiments tell us that there are 
three types of light standard-model electroweak-doublet neutrinos that couple to the W 
and Z boson. These are the three known neutrinos fe,/x,T- This implies that if there are 
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Table 1: Different possible conclusions regarding the nature of the neutrinos and their 
mass hierarchy from the three complementary experiments. 
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Conclusion 
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yes 
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no 


Degenerate Hierarchy, Majorana or 








Normal Hierarchy, Majorana with heavy particle contribution 


yes 
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no 


Inverted Hierarchy, Majorana 
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yes 


Degenerate Hierarchy, Majorana 


no 
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no 


Normal Hierarchy, Dirac or Majorana 


no 
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no 


Dirac 


no 


< 


yes 


Dirac 


no 


> 


yes 


Dirac 



other neutrino-like interaction eigenstates, then they must either be sufficiently massive 
that they cannot occur in the decay of the Z or they must be electroweak singlets with 
no coupling to the W or Z. In the latter case, the interaction eigenstates are called 
sterile neutrinos. In general, a neutrino mass eigenstate will be a linear combination of 
the three electroweak-doublet neutrinos and some unknown number of electroweak-singlet 
(= sterile) neutrinos. In the presence of electroweak-singlet neutrinos, the neutral weak 
current is not, in general, diagonal [40,41]. In common parlance, the word sterile neutrino 
is often used to denote a light electroweak-singlet neutrino and hence to exclude the heavy 
electroweak-singlet neutrino- like states that may well play a role in the seesaw mechanism. 
So the question is: are there any (light) sterile neutrinos and if so, how many are they 
and do they mix with the ordinary neutrinos? 

Light sterile neutrinos have been postulated in order to explain [42] the data from the 
Los Alamos Liquid Scintillation Detector (LSND) experiment [43], where neutrino flavor 
conversion both from a stopped muon (DAR) as well as the one accompanying the muon 
in pion decay have apparently been observed. The evidence from the DAR is statistically 
more significant and is interpreted as an oscillation from to Pg- The mass and mixing 
parameter range that fits data is: 

Am^ -0.2-2 eV^ , sin^ 29 ~ 0.003 - 0.03 . (2) 

There are points at higher masses specifically at 6 cV^ which are also allowed by the 
present LSND data for small mixings. The KARMEN experiment at the Rutherford lab- 
oratory has very strongly constrained the allowed parameter range of the LSND data [44] . 
Currently the MiniBooNE experiment at Fermilab is under way to probe the LSND pa- 
rameter region [45]. 

Since this AmLgNo is so different from that Ami, ^, the simplest way to explain these 
results is to add one [42, 46] or two [47] sterile neutrinos. The sterile neutrinos raise 
important issues of consistency with cosmology as well as physics beyond the simple 
three neutrino picture and will be discussed in a subsequent section. 
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1.1.11 Neutrino electromagnetic dipole moments and neutrino decay 

A massive Dirac neutrino can have a diagonal magnetic (and a CP-violating electric) 
dipole moment. Because a Majorana neutrino is the same as its antiparticle, it has van- 
ishing diagonal magnetic and electric dipole moments. A massive Dirac or Majorana 
neutrino can have nondiagonal, i.e., transition, magnetic and electric dipole moments. 
Some discussions of diagonal and transition neutrino electromagnetic moments in renor- 
malizable electroweak gauge theories (where these can be calculated) include [48,49], 
[40], [51]- [59]. In the standard model extended to contain massive Dirac neutrinos, 
fiuj = SeGprn^. /{8n'^V2) = 1.6 x 10~^^{m„. /{I eV)) fiB for the neutrino mass eigenstate 
Uj [52], where fiB = 2m~ ^ Bohr magneton. In left-right models and others with 
new physics beyond the standard model, this may be larger (e.g., [51,53,57-59]). In 
contrast to the magnetic dipole moment, the neutrino electric dipole moment vanishes at 
one-loop order for a massive Dirac neutrino in the extended standard model [53]. How- 
ever, in a left-right model, a Dirac neutrino may acquire an electric dipole moment at 
the one-loop level [53]. In the more generic case of a Majorana neutrino, one's interest 
focuses on the neutrino transition magnetic (and electric) dipole moments. The pres- 
ence of these diagonal or transition moments allows for new electromagnetic interactions 
between neutrinos and other fermions of the Standard Model. In particular in neutrino- 
electron scattering, in addition to the usual weak interaction contribution, there will be a 
photon exchange contribution to the scattering cross section. The existing neutrino scat- 
tering measurements therefore provide an upper limit on the neutrino magnetic moment: 
/i^^ < (1 — 1.3) X 10~^^fiB- As we will discuss in more detail later, the observation of 
nonzero neutrino magnetic moment would be considered evidence of new physics at the 
TeV scale. The reason for that is that if all new physics is parameterized by (Majorana or 
Dirac) neutrino masses, or, equivalently, if all new physics effects are suppressed by the 
very large naive seesaw energy scale (close to the GUT scale) the neutrino magnetic mo- 
ments are expected to be of order 10^^^ hb (f^)- High-precision searches for a magnetic 
moment provide, therefore, complementary tools to probe the physics that is expected to 
lie just beyond the electroweak symmetry breaking scale. 

A neutrino magnetic or electric dipole moment leads to new processes that can alter 
our understanding of energy balance in astrophysical systems such as in stars and super- 
novae [60]. It can also affect considerations involving the neutrinos in the early universe 
such as the BBN. In Sec. 16. II we discuss more details on neutrino magnetic moments and 
what one can learn from various proposed experiments. 

The existence of a neutrino magnetic or electric transition moment is also related 
to neutrino decays. For instance, it would allow heavier neutrinos to decay radiatively 
to the lighter ones [40,48,50,54,61,62]. Such decays can be detectable in astrophysical 
experiments. Present upper limits coupled with the general idea about spectra of neutrinos 
from oscillation experiments, imply that lifetimes of the primary mass eigenstates in 
electroweak-doublet neutrinos are larger than lO^'' sec, much longer than the age of the 
universe [40,48,50,54,61,62]. Such decays do not therefore affect the evolution of the 
universe. 

It is however possible that there are other scalar particles to which the neutrinos 
decay; one such example is the majoron, which is a Goldstone boson corresponding to the 



14 



spontaneous breaking of a global B — L symmetry [63]. The decay to these scalar bosons 
may occur at a faster rate [64] than that to photons and may therefore have astrophysical 
and cosmological implications [65] . This will be the subject of another working group [9] ; 
so we only focus on the implications of the magnetic moment in one of the subsequent 
sections. 

1.2 Neutrino probes of other fundamental symmetries 

Neutrino experiments can also be used to probe the validity of other fundamental sym- 
metries, some of which are often commonly assumed in theoretical discussions, as well 
as the basic assumptions of local quantum field theories on which the Standard Model is 
based. Some examples of these are: 

• Violation of Lorentz invariance; 

• CPT violation; 

• Possible existence of new long range forces in nature associated with lepton number; 

• Nonstandard interactions of neutrinos such as flavor changing neutral currents in- 
volving neutrinos. 

We will explore to what extent existing limits on these departures from standard 
scenarios can be improved. 

1.3 Why neutrino mass necessarily means physics beyond the 
Standard Model? 

Neutrino oscillations are, to date, the only evidence for the existence of physics beyond the 
Standard Model (in the domain of particle physics). It is of utmost importance to decipher 
the kind of new physics indicated by the existing data and to anticipate the signals of new 
physics that might appear in future planned observations. We must understand how and 
if they fit into the different big pictures that have been advocated for independent reasons, 
including the gauge hierarchy problem and gauge coupling unification. To discuss this, 
we first introduce the Standard Model and possible ways to extend it to accommodate 
the neutrino observations. 

In the Standard Model, which is based on the gauge group SU (3)c x SU {2)l x U (l)y , 
the quarks and leptons transform as (5^,(3, 2, |), 14^(3, dij(3, 1,— |), L(l,2,— 1), 
6^(1, 1, —2). The Higgs boson H, responsible for electroweak symmetry breaking, trans- 
forms as (1, 2, +1). The electroweak symmetry SU{2)l x f/(l)y is broken by the vacuum 
expectation of the Higgs doublet (H^) = Vwk — 246 GeV, which renders the and 
gauge bosons and the electrically charged fermions massive. The reason neutrinos do not 
get mass as a result of the Higgs mechanism is that the right-handed neutrino Nr was not 
included in the list of fermions in the Standard Model; as a result there is no coupling of 
the form h^LHN^ that could have given mass to the neutrinos after symmetry breaking. 

One seemingly straightforward way to understand the neutrino mass would be to 
extend the Standard Model to include the Nr. This would also be desirable from the 
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point of view of making the model quark lepton symmetric. There are two problems with 
this naively trivial modification. One is that by quark lepton symmetry one would expect 
the neutrino masses arising from the Yukawa couphng h^LHNji to be of the same order as 
the quark and charged leptons. Observations suggest that neutrino masses are at least 10^ 
times smaller than the smallest quark and lepton masses. Therefore, a nonzero neutrino 
mass not only suggests the existence of right-handed neutrinos (of which there would be 
three if they correspond to the usual generations), but some new physics that will enable 
us to understand why M,^ <C rriq^g. The seesaw mechanism provides a plausible basis for 
this understanding, since it makes use of the fact that, among the known fermions, only 
neutrinos can have Majorana mass terms. Thus, ironically, we may have a better way 
to understand the lightness of the neutrinos than we do to understand the generational 
hierarchy factor of ~ 10^ between the masses of the top quark and the electron, for which 
there is no accepted explanation at present. 

The other problem with introducing a set of right-handed neutrino fields is the fact 
that they arc Standard Model gauge singlets. This means that, as far as the symmetries of 
the Standard Model are concerned, a Majorana mass for the fields is allowed. If such a 
mass term is present, however, the neutrino masses arc not simply given by the h^v^i^, but 
are determined by a more complicated function of h^v^i^ and the Majorana masses of the 
right-handed neutrinos. In order to avoid the presence of a Majorana mass for the right- 
handed neutrinos one is required to impose an extra symmetry to the Standard Model 
Lagrangian (say, lepton number) — a very nontrivial modification of what is traditionally 
referred to as the Standard Model of electroweak interactions. 

1.3.1 Seesaw mechanism for small neutrino masses 

A simple way to understand the smallness of neutrino mass within this minimally extended 
Standard Model is to break lepton number symmetry (or S — L symmetry) and add a 
Majorana mass for the right-handed neutrino MrA^^C~^A^r. Thus the two terms that 
give mass to the neutrinos have the form h^v^yJ^LNji + MfiN'^C~^Nji + h.c. Assuming 
n "left-handed" neutrinos ul and m "right-handed neutrinos Nr, the (n + m.) x {n + m) 
Majorana neutrino mass matrix is 



In the limit 3> hv^/k, the eigenvalues of this matrix are given by — j^^^ and M^j, 
with respective approximate eigenvectors ul and Nji. The effective active neutrino masses 
are clearly much smaller than typical charged fermion masses (which are of order h^v^^ 
as long Mr ^ f„k- This is the well-known seesaw mechanism [66-70]. If we take as a 
guide a value for h^y < 1, then atmospheric neutrino data requires that M/j < 10^^ GeV. 
It should be emphasized that there is very little concrete information or experimental 
guidance regarding the magnitude of M^, which is virtually unconstrained [71]. One 
question which arises is why this value rather than Mpi ~ lO-*^^ GeV, which, one may 
argue, would have been a more natural value? Could this be an indication of a new 
symmetry? The answer to this question is obviously of fundamental significance. 




(3) 
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An example of such a symmetry is the B — L symmetry embodied in the left-right 
symmetric models based on the gauge group SU (2) ^ x SU{2)r x U{1) b-l [72]. This gauge 
group is also a subgroup of 50(10) grand unification group. The above mentioned value 
of Mr is rather close to the conventional GUT scale of 10^^ GeV. This makes the seesaw 
mechanism a very attractive framework for discussing the neutrino mass. We will discuss 
further consequences of grand unification for neutrino masses in a subsequent section. We 
will also explore in this review unification-model independent consequences of the seesaw 
mechanism. 

1.3.2 Type I vs type II seesaw mechanism 

If there are indeed right-handed neutrinos, the most general Majorana mass matrix that 
mixes active and sterile neutrinos is given by Eq. (jH)) with the n x n zero matrix in the 
upper left-hand corner replaced by a generic (symmetric) matrix Ml- This phenomenon 
occurs, for example, when the theory containing the becomes parity symmetric as is 
the case for SU{2)l x SU{2)r x U{X)b-l or 5*0(10) based models. In this case the seesaw 
formula is modified to 

= Ml- M^M^\Mj^f , (4) 

where, in an SU{2)l x SU{2)ji x U{1)b-l symmetric model. Ml = fvL and M/j = fvR, 
where vl,r are the vacuum expectation values of Higgs fields that couple to the right and 
left-handed neutrinos. Eq. (jlj) is called the type II seesaw relation [73] . It should be noted 
that in the absence of a discrete left-right symmetry. Ml is in general not related to Mr. 

1.3.3 Triplet seesaw 

An alternative way to understand the small neutrino mass without introducing the right 
handed neutrino is the triplet seesaw mechanism. It was pointed out in various papers [74] 
in early 1980 that if the standard model is extended by the addition of a triplet Higgs Al 
with weak hypercharge Y = 2, a vev for it can lead to Majorana mass for the neutrinos 
from the interaction fu4'LC^^T2'r ■ AtpL (V'l being the lepton doublet iyL^^Lj)- However, 
one has to tune the Yukawa coupling fy by about 10~^° or so to get desirable neutrino 
masses. It has subsequently been shown [75] that in the context of grand unified theories, 

the triplet vev is given by the formula < A^, >~ where Mu is close to the grand 
unification scale and corresponds to the physical mass of the triplet Higgs field. Since 
Mu ^ Vwk, this provides a natural suppression of the triplet vev and the right order for 
the neutrino mass emerges. Note also that this can also emerge from the type II seesaw 
formula in the limit of M^r^ —>■ oo. In this case, the neutrino mass matrix is directly 
proportional to the coupling matrix f^. 

1.3.4 Seesaw with triplet fermions 

Yet another possible extension of the standard model without right handed neutrinos 
which leads to small neutrino masses is to postulate the existence of triplet vectorlike 
fermions: A [76]. Since a vectorlike triplet can have an arbitrary mass, it also leads to 
seesaw mass formula. 
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1.3.5 Understanding Icirge mixings within the seesaw mechanism 

A major puzzle of quark-lepton physics is the fact that the quark mixing matrix and the 
leptonic one are quahtatively different. In order to understand the mixing angles [4,77], 
we have to study the mass matrices for the charged leptons and neutrinos. 
A general approach 

To see the possible origin of neutrino mixings, one can start with the following form 
for the mass part of the neutrino Lagrangian: 

jO-mass = i^LM^NR + eLM,eR + NIMrNr + h.c. (5) 

Using the seesaw mechanism one can derive from this equation, the formula for neutrino 
masses can be written as for the case of type I seesaw: 

Ml = -M^M^'iM^f, (6) 

To obtain the lepton mixing matrix, one can diagonalize the charged lepton mass matrix 
by Me = t/^M/Ft and = {U*)^mt{W)^ and find that U = U}U^. 

With this theoretical preamble, understanding of neutrino mixings can proceed along 
two paths. In theories where quark and lepton mixings are disconnected (such as many 
weak scale theories), one may like to pass to a basis where the charged lepton masses are 
diagonal. In that case, all the neutrino mixing information is in the effective neutrino 
mass matrix. One can then look for the types of mass matrices for neutrinos that can 
lead to bi-large mixings and try to understand them in terms of new physics. Here we 
give a brief overview of some generic structures for M,y that do the job. 
(i) The case of normal hierarchy: In this case, one neutrino mass matrix that leads to 
"bi-large" mixing has the form: 

= mo I e 1 + e 1 I (7) 




where mo is a/ Am\. We have omitted order one coefficients in front of the e's. This 
matrix leads to tan 9 a — 1, AuIq/ Am\ ~ and also to a large solar angle. For the LMA 
solution, we find the interesting result that e ~ A where A is the Cabibbo angle (~ 0.22). 
This could be a signal of hidden quark lepton connection^. In fact we will see below that 
in the context of a minimal SO{10) model, this connection is realized in a natural manner, 
(ii) The case of inverted hierarchy: The elements of the neutrino mass matrix in this case 
have a pattern 



mj, = mo c e e . (8) 




where c = cos6 and s = sm6 and it denotes the atmospheric neutrino mixing angle. 
An interesting point about this mass matrix is that in the limit of e ^ 0, it possesses 

^Alternatively, the relation 0q + X = 7r/4, nowadays known as quark-lepton-complementarity, can 
also be interpreted as such a connection [78]. 
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an Lf. — — symmetry [79]. One therefore might hope that if inverted hierarchy 
structure is confirmed, it may provide evidence for this leptonic symmetry. This can be 
an important clue to new physics beyond the Standard ModeL This issue of leptonic 
symmetries will be discussed in the main body of this report. 

(iii) Degenerate neutrinos: One may either add a unit matrix to the just mentioned mass 
matrices and look for new physics models for them; alternatively, one may look for some 
dynamical ways by which large mixings can arise. It turns out that if neutrinos are 
mass degenerate, one can generate large mixings out of small mixings [80-84] purely as a 
consequence of radiative corrections. We will call this possibility radiative magnification 
and will discuss it in a future section. 

In grand unified theories, quark and lepton mass matrices are connected. One may 
therefore lose crucial information about symmetries if one works in a basis where the 
charged leptons are diagonal. Furthermore, if either of the quark (up or down) mass 
matrices are chosen diagonal, it may not even be possible to go to the diagonal charged 
lepton basis. Thus in this case, we have U = UlUu- So one may seek an understanding 
of large mixings in the charged lepton sector. For example in SU{b) type theories, Ut is 
related to the mixings of right-handed quarks which are unobservable in low energy weak 
interactions and can therefore be the source of large mixings. Models of this type are 
called lopsided mixing models [85] . 

The basic strategy then would be to look for clues for new symmetries in the structure 
of the mass matrices, which could then provide information about the nature of physics 
beyond the Standard Model. The symmetries of course may become obscured by our 
choice of basis where the charged leptons are diagonal. It is this which gives different 
possibilities for arriving at the bi-large mixings and the hope is that different strategies 
will lead to different predictions for observables, which can then be put to experimental 
test. 

1.3.6 Alternatives to high-scale seesaw 

While the high-scale seesaw mechanism is the simplest and perhaps the most elegant 
way to understand the small neutrino masses and become couched in a quark leptonic 
and parity symmetric framework leading to simple grand unification theories, there are 
alternatives to seesaw which can also explain the small neutrino masses [86-95]. In such 
a case the neutrinos can either be Dirac or Majorana fermions depending on the theory. 

Unlike the non-supersymmetric seesaw models, alternatives such as the one presented 
in [90] often predict observable charged lepton lepton-fiavor violating signals [96], e.g., 
II —>■ 67, II, T ^ eee, etc. More generically, searches for charged-lepton flavor violation 
crucially help distinguish among the several theoretical interpretations of the origin of 
neutrino masses. 

1.4 Summary of the Introduction 

Some of the questions that we would like to answers in the course of this work are: 
• Can we decide whether the neutrino is a Dirac or Majorana particle? 
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To what extent can the planned neutrino experiments pin down the structure of the 
three neutrino mass matrix? This involves such questions as determining the sign 
of Am23, higher precision measurement of mixing parameters, etc. 

What is the impact of a ^^13 measurement (and the improved determination of the 
other elements of the lepton mixing matrix) on the general landscape of physics 
beyond the Standard Model? We find that ^13 is a powerful discriminator of models. 

Can we test the seesaw hypothesis and discriminate between different types of seesaw 
using lepton flavor violation and other "non-neutrino probes?" 

How can one experimentally discover or limit physics beyond the standard scenario? 
This will address such aspects as: 

(1) Flavor changing neutral currents for neutrinos; present limits and future prospects; 

(2) Admixtures of sterile neutrinos, both heavy and light; 

(3) magnetic moments of neutrinos. 

What can we learn about CP violation in the lepton sector and how can we connect 
it to the question of the origin of matter via leptogenesis. Given what we know about 
the neutrino masses, assuming thermal leptogenesis, do we have an explanation of 
the observed baryon to photon ratio? 
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2 What can we learn about neutrino mass matrices 
from experiments? 

In this section we briefly review our ability to reconstruct the neutrino mass matrix. We 
will also discuss (from "the bottom up") what we hope to learn from the neutrino mass 
matrix itself, instead of trying to quantify what different models predict for the neutrino 
mass matrix. See, for example, [4,77] for reviews of a few different models. In a subsequent 
section, we will discuss the connection of neutrino masses to GUTs, and will spend a little 
more time on "top-down" predictions for neutrino masses and mixing angles. 

As mentioned earlier, we will assume that the neutrinos are Majorana fermions. While 
there is no experimental evidence that this is the case, the majority of the theoretical HEP 
community considers it more likely that the neutrinos are Majorana fermions, and a larger 
amount of phenomenological research effort has gone into understanding and interpreting 
Majorana neutrino mass matrices than Dirac mass matrices. For some discussions of 
Dirac neutrino mass matrices and how they are related to the large mixing in the leptonic 
sector and the neutrino mass-squared differences, see, for example, [97]. 

Below the electroweak phase transition, the Majorana neutrino mass matrix m^, is the 
coefficient of the operator (using four-component-spinor notation) 



where a, f3 = e, ij,,t, . . . are flavor indices, and m"^ are the components of the neutrino 
mass matrix (note that rrij, is symmetric, i.e., m^^ = m(^"). In this section we will 
concentrate on a purely active 3x3 mass matrix. A detailed discussion of 4 x 4 (and 
larger) mass matrices, which also allow for the existence of fourth generation and/or sterile 
neutrinos is the subject of subsequent sections. Note that Eq. is not sensitive to the 
mechanism that generates neutrino masses. These will be discussed in more detail in a 
later section. 

In general, one cannot work back from a knowledge of the observed lepton mixing 
matrix to the individual nondiagonal mass matrices in the charged lepton and neutrino 
sectors. It is, indeed, the diagonalization of both of these mass matrices that gives rise 
to the observed lepton mixing, and models exist where the mixing in the charged lepton 
sector is large. One can always choose to work on the weak basis where the charged lepton 
mass matrix is diagonal — the price one pays for doing this is that the flavor structure of 
the theory may not be manifest. In this case, one can calculate the neutrino mass matrix 
in terms of the observed lepton mixing matrix as 



We choose sign conventions such that the neutrino mass eigenvalues are real and positive. 
By choosing to write U = VK, where V and K are given by Eq. (1) we have removed 
all of the redundancy contained in rriy associated with re-defining the neutrino fields by 
a complex phase. Hence, rriu as defined by Eq. ()10|) is only a function of observable 
parameters. The phases in K are the so-called Majorana phases [11,12]. They can be 



(9) 




(10) 
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redefined away by allowing the neutrino mass eigenvalues to be complex. In this case, 
U = V, mi is real and positive, and m2 = |m2|e^^*'^i, = \m3\e~'^^'^^ . 

In the near future, we hope to significantly improve the determination of the elements 
of the neutrino mass matrix, although some uncertainty will still remain (for a detailed 
discussion, see, for example, [98]). Through neutrino oscillation experiments, all three 
mixing angles 612,023, and 9i^, are expected to be determined with good precision (this 
is one of the main goals of next-generation neutrino oscillation experiments, discussed in 
great lengths in this report), while there is hope that the "Dirac phase" S can be probed 
via long-baseline oscillation searches. Neutrino oscillation experiments will also 

determine with good precision the neutrino mass-squared differences (Am^g at the 5%- 
10% level, Am^3 [including the sign] at the few percent level). In order to complete the 
picture, three other quantities must also be measured, none of which is directly related 
to neutrino oscillations. 

One is the overall scale for neutrino masses. As already briefly discussed, this will 
be probed, according to our current understanding, by studies of the end-point spectrum 
of beta-decay, searches for neutrinoless double beta decay, and cosmological observations 
(especially studies of large-scale structure formation). Note that neutrinoless double beta 
decay experiments are sensitive to |m^^|, i.e., they directly measure the absolute value of 
an element of m^,. 

The other two remaining observables are the "Majorana" phases. Neutrinoless double 
beta decay experiments are sensitive to a particular combination of these, the so-called 
effective Majorana mass, 

\nf^\ = {m)eff = |cos2 ^13 (|mi| cos^ 612 + \m2\e^''^' sin^ ^12) + sin^ 9i3\ms\e^''^^ \ . (11) 

With present uncertainties in the nuclear matrix elements, however, it seems at least very 
challenging [99] to obtain any information regarding Majorana phases from neutrinoless 
double beta decay. For a detailed study, see, for example, [100]. 

A few comments are in order. First, the relation between the rate for neutrinoless 
double beta decay and the Majorana phases and neutrino masses only holds under the 
assumption that the neutrino masses are the only source of lepton-number violation (as 
far as neutrinoless double beta decay is concerned). Second, only one or a combination 
of the two independent Majorana phases can be determined in this way. It is fair to 
say that there is no realistic measurement one can look forward to making in the near 
future that will add any information and help us disentangle the "other" Majorana phase. 
Third, it is curious to note that the effect the Majorana phases have on the rate for 
neutrinoless double beta decay is CP even [37] . While Majorana phases can mediate CP 
violating phenomena [37], it seems unlikely that any of them can be realistically studied 
experimentally in the foreseeable future. For further discussion of CP violation among 
neutrinos see Ref. [101]. 

In spite of all the uncertainty due to our inability to measure Majorana phases, it 
is fair to say that we expect to correctly reconstruct several features of the neutrino 
mass matrix [98], especially if the overall mass scale and the neutrino mass hierarchy 
are determined experimentally. What do we hope to accomplish by reconstructing the 
neutrino mass matrix? The answer is that we wish to uncover whether there are new 
fundamental organizing principles responsible for explaining in a more satisfying way 
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the values of the neutrino masses and the leptonic mixing angles. In other words, we 
would like to establish whether there is a fundamental reason behind the fact that the 
state contains almost the same amount of i/^ and u^-, while at the same time containing 
a relatively small amount of u^. Are there flavor (or family) symmetries, capable of 
dynamically distinguishing the different generations of quarks and leptons and, we hope, 
explaining why there are three quasi-identical particles for each matter field? 

In the neutrino sector, we are only getting started. We have, for example, identified 
several textures for the neutrino mass matrix that lead to the currently observed mass- 
squared differences and mixing angles, and have identified some of the measurements 
that will allow us to identify which textures best describe Nature. As has been already 
pointed out, it is not clear whether this is the best avenue to pursue as far as identifying 
whether there is a deep explanation for the patterns we observe in experiments. For 
example, it may turn out that we have made a weak-basis choice that renders the job 
more complicated (it is possible that the mixing angles are "determined by the charged 
lepton sector" [102-104]), or that all the structure contained in the neutrino sector is 
obscured after heavy degrees of freedom are integrated out (as may happen in type-I 
seesaw models). Nonetheless, we will discuss a few of these textures in order to exemplify 
some of the measurements (and how precise they should be) that will shed a significant 
amount of light in the issue of interpreting neutrino masses and mixing angles. 

Arguably the simplest assumption one can make is that there is no symmetry or 
dynamical principle that explains why leptonic mixing angles are large [105]. This "fla- 
vorless" neutrino flavor model is often referred to as "neutrino mass anarchy" and is, 
currently, compatible with data [105,106]. Curiously, the anarchical hypothesis is not 
without predictions: it requires that the unobserved magnitude of the Ues element of the 
leptonic mixing matrix is |f/e3p > 0.01 at the 95% confidence level (see [106] for details 
and a proper definition of this bound). This means that after the next-generation of re- 
actor and/or long-baseline experiments analyze their data we will know whether we can 
afford a "random" leptonic mixing matrix or not. It should be noted that this model 
applies only for the leptonic mixing matrix — it has nothing specific to say about the 
order of magnitude of neutral and charged lepton masses, or their hierarchies. 

If one assumes that there is a nontrivial texture to the neutrino mass matrix, and 
that this texture "explains" the observed values of the mixing parameters, there are 
several completely different options. Some are tabulated in Table El and will be discussed 
briefly. Before proceeding, however, it is important to explain how these textures should 
be interpreted. The hypothesis is that, at leading order, the neutrino mass matrix can 
be parametrized by far fewer than the usual six complex coefficients. These are chosen 
in such a way that the dominant features of neutrino masses and mixings are explained. 
These are: (i) mg is either much larger or much smaller than m\,m^. This splitting 
determines the atmospheric mass-squared difference, (ii) the Ve content of the z/3 state 
is zero, (iii) the and v^- contents of the z/3 state are similar (or, perhaps, identical). 
In order to accommodate the other observed features (like If/gsl, the solar mass-squared 
difference and the solar angle) one includes sub- leading effects that violate the leading- 
order structure. The structure of these sub-leading effects determines the "predictions" 
for the observables that are not determined by the leading order mass-texture. In Table |2l 
we list the predictions obtained in the case of a structureless sub- leading mass matrix, i.e., 
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Case 


Texture 


Hierarchy 


\Ue3\ 


COS 2^23 

(n.s.) 


cos 26123 


Solar Angle 
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/ \ 

1 1 
lo 1 1 J 
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m — o~ 

/ Amf2 


0{1) 


m — — 

/Amfj 


0{1) 
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y Ami, 
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V ^'""is 1 


1 u 
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-i 1^ 
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Inverted 


Amjj 
|Am^3l 


— 


Amfj 
|Am^3l 


0{1) 


C 




/ \ 

1 
\l 0^ 




Inverted 


Amjj 
\Am-i,\ 


0{1) 


Amfj 
|Am?3| 


1 COS Mi2\ 

~ 4^ 

|Am^3l 


Anarchy 




{\\\. 


) 


Normal 


> 0.1 


0{1) 




0{1) 



Table 2: Different leading-order neutrino mass textures and their "predictions" for various 
observables. The fifth column indicates the "prediction" for | cos 2^23] when there is 
no symmetry relating the different order one entries of the leading-order texture ('n.s.' 
stands for 'no structure', meaning that the entries of the matrices in the second column 
should all be multiplied by and order one coefficient), while the sixth column indicates 
the "prediction" for | cos 26*23! when the coefficients of the leading order texture are indeed 
related as prescribed by the matrix contained in the second column. See text for details. 
One may argue that the anarchical texture prefers but does not require a normal mass 
hierarchy. 

one proportional to the anarchical texture [107]. In the case of a structured sub-leading 
mass matrix, expectations may vary significantly from these quoted in Table El 

Case A is characterized by large entries in the "/i — r" sub-matrix, and small entries 
in the "e" column and row. The determinant of the "/i — r" sub-matrix is constrained 
to be small in order to guarantee a hierarchy between the two independent mass-squared 
differences. The hierarchy of the neutrino masses is predicted to be normal (mg > m"^ > 
ml). Maximal atmospheric mixing can be imposed at the leading order by requiring that 
the "/X — r" sub-matrix is democratic. The introduction of sub-leading effects leads to a 
"large" \Ue3\ and cos 26*23, of order the square-root of the ratio of mass-squared differences, 
which is (9(0.1). If this texture is indeed realized in nature, we expect to observe a 
nonzero jf/esl and a deviation of the atmospheric mixing from maximal at next-generation 
experiments. It may prove difficult to distinguish between case A and the anarchical 
texture via neutrino oscillation measurements alone. One potential discriminant seems to 
be the expected rate for neutrinoless double beta decay. 

Case B is characterized by small "e — /i" and "e — r" entries, a small determinant 
of the "/i — r" submatrix, and the constraint that the trace of m,^ is close to 2m^'^. In 
this case, one predicts an inverted mass hierarchy (mg > mf ^ ^3), and both {Ue^l and 
cos 26*23 are constrained to be of order the ratio of the mass-squared differences (C(O.Ol)). 
The system is constrained enough that it is hard to obtain a much larger deviation of 
the atmospheric angle from maximal or a much larger {Uesl, while much smaller ones 
are, of course, obtainable if the sub-leading contributions are structured. If case B is 
indeed realized in Nature, there is a good chance that no {Uesl effects will be observed in 
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next-generation oscillation experiments, while precise measurements of the atmospheric 
mixing angle will remain consistent with 6*23 = 7r/4 (equivalently, if a large deviation of 
the atmospheric angle is detected this texture will be ruled out) . On a more positive note, 
one should expect a "large" rate for neutrinoless double beta decay (m^^ ~ \J Amfg) . A 
texture which is naively similar to case B is to change the sign of the "/i — r" sub-matrix, 
such that the trace of the leading order mass matrix is close to zero. This case, however, 
is disfavored by solar data, as the solar angle is constrained to be too small (for a more 
detailed discussion see, for example, [107]). 

Case C is characterized by "e— and "e— r" entries which are much larger than all the 
other ones (set to zero at leading order) and thus corresponds to the case of approximate 
Le — L^j_ — Lr symmetry [79]. It leads to an inverted mass hierarchy, and a close to bi- 
maximal [108] leading order mixing matrix. The solar angle is (at leading order) exactly 
maximal, while the atmospheric angle is generically large, becoming maximal in the limit 
when ml'^ — ±.rrfj (for real m^"). Sub-leading corrections to case C which are responsible 
for splitting the two heavy leading-order mass eigenstates will induce a \Uez\-, cos 2^23 and 
cos 2^12 of order the ratio of the mass-squared differences ((9(0.01)), or smaller. Hence, 
similar to case B, it seems unlikely that t/es-^ffects will be measured at next-generation 
experiments. This scenario is currently disfavored, as it also predicts a solar angle 612 
very close to 7r/4 [102, 107, 109]. One should not conclude, however, that scenarios based 
on "perturbations" around bi-maximal mixing are ruled out. A related issue has been 
discussed in detail recently by the authors of [102, 103]. 

A realistic three generation extension of the mass matrix in case (A) that leads to 
large but not maximal solar neutrino mixing is given by 



= he 1 + ce -1 (12) 




Note that iih = d and c = 1, the atmospheric neutrino mixing is maximal and the mixing 
parameter ^13 — [110]. In this case the mass matrix has n — t interchange symmetry. 

Depending on how this symmetry is broken, the parameter ^13 is either of order ^J^^- 

(for c 7^ 1) or . jf (for h ^ d) [111]. Therefore search for ^13 down to the level of 0.01 
will be a big help in determining the structure of the neutrino mass matrix for the case 
of normal hierarchy. In both these cases, there is an important correlation between ^13 
and Oa-I [111]. 

In general, a neutrino mixing matrix originating from a, /ji — r symmetric mass matrix 
has the following structure (for simplicity, we did not include here the Majorana phases) 
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(13) 



Note that the mass spectrum of the neutrinos is not predicted by the — r symmetry. 
Depending on the value of ^12, several interesting mixing schemes can arise: if ^12 = 7r/4 
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Figure 2: Departure from /i — r symmetry and correlation between ^13 and 9a (Case (i)): 
opposite to case (ii) (figure by H. Yu) 



then we have bi-maximal mixing. However, as mentioned, the observed deviation from 
7r/4 is rather large. 

Much closer to current data is the so-called tri-bimaximal mixing scheme [112], cor- 
responding to sin^ 9i2 = 1/3 and leading to the following often-studied mixing matrix: 



U 



( 



\ 







(14) 



Models which give rise to such a matrix (for some recent attempts see [113]) are typically 
quite intricate and not as straightforward to construct as models leading to bi-maximal 
mixing. 

There are other viable neutrino mass textures, including some that lead to degenerate 
neutrino masses. We refer readers to the literature for a more thorough discussion ( [4, 77] 
and references therein). The point we wish to emphasize here is that the amount of 
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information we have concerning neutrino masses and leptonic mixing is still very limited. 
This is reflected in the fact that too many different hypothesis can be raised in order 
to "explain" the same set of observables. The situation is bound to change in the near 
future, and there is hope that the data will "select" one specific neutrino mass matrix. 
Our job will then be to interpret what Nature is trying to "say" through rrii,. A more 
accurate determination of a few observables will already shine a significant amoimt of 
light in the currently obscure picture we arc trying to obtain: (i) what is the neutrino 
mass hierarchy? (ii) is l^/gsl larger than 0.1? (ui) is |cos26'23| > 0.1? All three of these 
can be answered in a next- generation i/^ — > i/g long-baseline accelerator experiment, while 
the second one can be addressed by a next-generation reactor neutrino experiment (with 
a baseline of 0{1 km)). Of course, in order to be sure we are on the right track we also 
need to (iv) determine that lepton number is not a conserved symmetry. 
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3 /3/3ojy-Decay and CP Violation 



3.1 (3f3ov decay 

In this section, we focus our attention on what we can learn from neutrinoless double 
beta decay experiments. As already alluded to in the introduction, in a given theory, 
neutrinoless double beta decay can arise from two sources: (i) neutrino Majorana mass 
or/and (ii) lepton number violating interactions. While the absence of a signal in a P(3ou 
experiment will constrain both sources (and associated theories), a positive signal cannot 
necessarily be considered as evidence for one or the other exclusively. For instance, one 
must supplement the results from /3/3oi/ experiment with collider experiments such as 
from LHC or a possible e^e~ type to get a definitive information regarding the source. 
Alternatively, one may conduct the double beta decay experiment with different nuclei 
and if the matrix elements happen to differ substantially for the two sources, one may 
be able to disentangle the different sources. Therefore as we interpret any positive signal 
for PPqi, decay one must keep this in mind. Below, we discuss what we can learn about 
neutrino masses and mixings, once it is established that the source of the positive signal 
is the Majorana mass for the neutrino. The experiments with solar and atmospheric 
neutrinos and with reactor antineutrinos have provided data on 6*12, 6^23 and 6*13, and on the 
neutrino mass squared differences driving the solar and atmospheric neutrino oscillations, 
Am^2 Am^g. Future oscillation experiments will improve considerably the precision 
on these basic parameters. However, these experiments are insensitive to the nature of 
massive neutrinos uj which can be Dirac or Majorana particles [11,114] (see also, e.g., [23]). 
They cannot give information on the absolute scale of neutrino masses (i.e., on the value 
of mi), and on the two Majorana CP violation phases — the latter do not enter into 
the expressions for the probabilities of flavor neutrino oscillations [11]. Determining the 
nature of massive neutrinos and obtaining information on the absolute neutrino mass 
scale is one of the fundamental problems in the studies of neutrino mixing. 

Neutrinos are predicted to be Majorana particles in the seesaw model of neutrino mass 
generation. This model gives a natural explanation of the smallness of the neutrino masses 
and, through the leptogenesis theory, provides an explanation of the observed baryon 
asymmetry in the Universe, which thus is linked to the existence of neutrino mixing. The 
only experiments which have the potential of establishing the Majorana nature of massive 
neutrinos are the neutrinoless double beta decay experiments searching for the process 
{A, Z) — > (A, Z + 2) +e~ + e~ (see, e.g., [23,24]). The observation of the /5/5o!/-decay and 
the measurement of the corresponding /9/9o!/-decay rate with a sufficient accuracy, would 
not only be a proof that the total lepton charge is not conserved in nature, but might 
provide also a unique information on i) the type of the neutrino mass spectrum, ii) the 
absolute scale of neutrino masses, and on iii) the values of the Majorana CP violation 
phases. Let us add that in supersymmetric theories with seesaw mechanism of neutrino 
mass generation, the rates of lepton flavor violating decays /x— s>e + 7, r— >-/i + 7 can be 
interestingly large (e.g., [116]) and may depend on the Majorana CP violating phases in 
the lepton mixing matrix (see, e.g., [117,118]). Furthermore, the values of the Majorana 
phases can be important for the stability under RGE running of the neutrino mass and 
mixing parameters, see Sec. 15.31 
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Let us recall that the SK atmospheric neutrino and K2K data do not allow one to 
determine the sign of Am\. This implies that if we identify Am\ with Amf^ in the case of 
3-neutrino mixing, one can have Am^g > or Amfg < 0. The two possibilities correspond 
to two different types of neutrino mass spectrum: with normal hierarchy, mi < m2 < m^, 
and with inverted hierarchy, ma < mi < m2. In the case of strong inequalities between 
the masses, the spectra are called normal hierarchical (NH) and inverted hierarchical (IH). 
The NH and IH spectra correspond to mi -C 0.02 eV and -C 0.02 eV, respectively. If 
mi = m2 = m3 = mo and m^ 3> |Am^|, Amg, the spectrum is quasi-degenerate (QD). 
The QD spectrum is realized if mi, 2,3 > 0.20 eV roughly requiring that the largest mass 
difference is about 10% of the common mass. 

Under the assumptions of (1) 3-neutrino mixing, for which we have compelling ev- 
idence from the experiments with solar and atmospheric neutrinos and from the Kam- 
LAND experiment, (2) massive neutrinos Uj being Majorana particles, and (3) (3 jSQi^-decay 
generated only by the (V-A) charged current weak interaction via the exchange of the three 
Majorana neutrinos Uj, the effective Majorana mass in /?/5o,y-decay of interest is given by 
(see, e.g., [23,119]): 

(m)e// = |mi|f/eip+m2|t/e2|' C^'^^ + mg | ?7e3 6^"*"] , (15) 

where Uej, j = 1,2,3, are the elements of the first row of the lepton mixing matrix U, 
rrij > is the mass of the Majorana neutrino uj, and 0i and 02 are the two Majorana CP 
violating phases [12,114]. In the case of CP conservation we have e^*''^^'^ = 7721(31) = ±1, 
rjij being the relative CP parity of the neutrinos z/j and uj. 

One can express [120] two of the three neutrino masses, say, m2,3, in terms of the 
third mass, mi, and of Am^ and Am^, while the elements \Uej\ can be expressed in 
terms of 6q and 6 (a concise discussion of the relevant formalism can be found, e.g., in 
Refs. [119,121,122]). Within the convention employed in the present study in both cases of 
neutrino mass spectrum with normal and inverted hierarchy one has: Ami, = Amf2 > 0, 
and m2 = \/mf + AttIq. In the case of normal hierarchy, Am^ = Am^g > and 
ms = ^ym^~+~Am^, while if the spectrum is with inverted hierarchy, Am^ = Am23 ^ 
and thus mi = ^/mg + Am^ — Am|. For both types of hierarchy, the following relations 
hold: \Uei\^ = cos2^0(l - It/esH, |f^c2p = sin2^0(l - If/esH, and |f/e3p = sin^e. We 
denote the smallest neutrino and we have mmin = "mi (3) for the case of 

normal (inverted) hierarchy. The problem of obtaining the allowed values of (m)e// 
given the constraints on the relevant parameters following from the neutrino oscillation 
data has been first studied in Ref. [120] and subsequently in a large number of papers, see, 
e.g., [39,119,123-126]. Detailed analyses were performed, e.g., in Refs. [29,100,119,122], 
and in particular in Ref. [127], where the allowed values of Am^, Am^, ^0 and 9, obtained 
from the most recent neutrino oscillation data, were used. The results are summarized in 
Tabs. El and H and in Fig. El 

In Fig. El (taken from Ref. [127]) we show the allowed ranges of values of {m)eff as 
a function of rrimm for the cases of NH and IH spectrum The predictions for {m)eff are 
obtained by using the allowed at 90% C.L. (Fig. EJ, values of Am|, 9q and Am^ from 
Refs. [128] and [130] and for three fixed values of sin^ 9. 

The existence of significant lower bounds on {m)eff in the cases of IH and QD spec- 
tra [131], {m)Yff > lOmeV and (m)^^^ > 43meV, respectively [127], which lie either 
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0.0 


2.6 (2.6) 
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4.6 (4.3) 


19.1 (16.6) 


48.5 (42.4) 


68.5 



Table 3: The maximal values of {m)eff (in units of meV) for the NH and IH spectra, 
and the minimal values of {'m)eff (in units of meV) for the IH and QD spectra, for 
the best fit values of the oscillation parameters and sin^ 9 = 0.0, 0.02 and 0.04. The 
results for the NH and IH spectra are obtained for |Am^| = 2.6 x 10"^ eV (2.0 x 
10~^ eV^ — values in brackets) and mi = 10~^ eV, while those for the QD spectrum 
correspond to rriQ = 0.2 eV. (From Ref. [127]). 
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11.4 (9.9) 


54.0 (48.6) 
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Table 4: The same as in Table El but for the 90% C.L. allowed regions of AttIq and 9q 
obtained in [128], and of Am\ obtained in Ref. [129] (Ref. [130] — results in brackets). 
(From Ref. [127]). 

partially (IH spectrum) or completely (QD spectrum) within the range of sensitivity of 
the next generation of /3/3ojy-decay experiments, is one of the most important features of 
the predictions of (m)ejj. The indicated minimal values are given, up to small corrections, 
by Am^ cos 2^0 and mo cos 2^©. According to the most recent combined analyses of the 
solar and reactor neutrino data, including the latest SNO and KamLAND results (see, 
e.g., [132]) i) the possibility of cos 26^0 = is excluded at more than 6 s.d., ii) the best 
fit value of cos 26^0 is cos 26^0 = 0.40, and ni) at 95% C.L. one has for sin^ 6^ = (0.02), 
cos 200 > 0.27 (0.32). The quoted results on cos26'0 together with the range of possi- 
ble values of Am^ and mg, lead to the conclusion about the existence of significant and 
robust lower bounds on (m.)e// in the cases of IH and QD spectrum [127,133]. At the 
same time one can always have {'m)eff = in the case of neutrino mass spectrum with 
normal hierarchy [134]. It follows from Tabs. IHl and ID that in this case (m)e// cannot 
exceed 5.5 meV. This implies that the maximal value of (m)e// in the case of neutrino 
mass spectrum with normal hierarchy is considerably smaller than the minimal values 
of (m)e// for the inverted hierarchy and quasi-degenerate neutrino mass spectrum. This 
opens the possibility of obtaining information about the type of neutrino mass spectrum 
from a measurement of (m)e// 7^ 0, or from obtaining a sufficiently stringent upper bound 
on (m)e//. In particular, a positive result in the future generation of PPou-decay exper- 
iments with {m)eff > 10 meV would imply that the NH spectrum is excluded. The 
uncertainty in the relevant nuclear matrix elements^ and prospective experimental errors 
in the values of the oscillation parameters, in (m)e//, and for the case of QD spectrum — 
in mo, weaken, but do not invalidate, the reported results (see, e.g., Refs. [136]). If the 

^Recently, encouraging results in what regards the problem of the calculation of the nuclear matrix 
elements have been obtained [135]. 
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Figure 3: The dependence of {m)eff on m^m in the case of the LMA-I solution, for 
normal and inverted hierarchy, and for the 90% CL. allowed regions of AttIq and siii^ 9q 
found in Ref. [128] and of Ami in Ref. [130], and a fixed value of sin^^ = 0.0(0.02) [0.04] 
in the upper (middle) [lower] panel. In the case of CP conservation, the allowed values 
of {m)eff are constrained to lie: for i) normal hierarchy and the middle and lower panels 
(upper panel) - in the medium-gray and light-gray regions a) between the two lower thick 
solid lines (between the two lower thick solid lines) if 7721 = 7731 — 1, b) between the two 
long-dashed lines (between the two lower thick solid lines) if 7721 = —Vsi = between 
the three thick dash-dotted lines and the axes (between the dash-dotted lines and the axes) 
if V21 — —Vsi — d) between the three thick short-dashed lines and the axes (between 
the dash-dotted lines and the axes) if 7721 = r^si = — 1; and for ii) inverted hierarchy and 
the middle and lower panels (upper) - in the light-gray regions a) between the two upper 
thick solid lines (between the two upper thick solid lines) if 7732 = 7/31 = ±1, b) between 
the dotted and the thin dash-dotted lines (between the dotted and the thick short-dashed 
lines) if 7732 = —7/31 = 1, c) between the dotted and the upper thick short-dashed lines 
(between the dotted and the thick short-dashed lines) if 7^32 = —7731 = — 1. In the case of 
CP violation, the allowed regions for {m)eff cover all the gray regions. Values of {m)eff 
in the dark gray regions signal CP violation. (From Ref. [127]). 
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neutrino mass spectrum turned out to be of the QD type, a measurement of {m)eff in 
PPo^-decaj experiment and of mg in the KATRIN experiment [22] could be used, in par- 
ticular, to check the validity of the light Majorana neutrino exchange mechanism for the 
PPou-decaj and search for indications for contributions from other types of mechanisms 
(see, e.g., [25,137]). 

It follows from Fig. IHl that a measurement of (m)e// > 10 meV would either i) 
determine a relatively narrow interval of possible values of the lightest neutrino mass 
rrimin, or ii) would establish an upper limit on the possible values of rrimin- If a sufficiently 
stringent upper limit on {■m)eff is experimentally obtained below 100 meV, this would 
lead to a significant upper limit on the possible value of rrimin- 

The possibility of establishing CP violation in the lepton sector due to Majorana CP 
violating phases has been studied in [126,134] and in much greater detail in Ref. [100]. It 
was found that it is very challenging^: it requires quite accurate measurements of (m)e// 
and of mi, and holds only for a limited range of values of the relevant parameters. For 
the IH and the QD spectra, which are of interest, the "just CP violation" region [119] — 
an experimental point in this region would signal unambiguously CP violation associated 
with Majorana neutrinos, is larger for smaller values of cos2^0. More specifically, proving 
that CP violation associated with Majorana neutrinos takes place requires, in particular, 
a relative experimental error on the measured value of (m)e// not bigger than (15 - 20)%, 
a "theoretical uncertainty" in the value of {m)eff due to an imprecise knowledge of the 
corresponding nuclear matrix elements smaller than a factor of 2, a value of tan^ 6q > 0.55, 
and values of the relevant Majorana CP violating phases (20i.2) typically within the ranges 
of ~ {tt/2 - 37r/4) and ~ (57r/4 - 37r/2) [100]. 

3.2 The MNSP Lepton Mixing Matrix and CP Violation in the 
Lepton Sector 

It is well known that, in general, in gauge theories with massive neutrinos the MNSP 
lepton mixing matrix results from a product of two matrices: 

U = Uj U, . (16) 

where Ui and Uu are two 3x3 unitary matrices: Ue arises from the diagonalization of 
the charged lepton mass matrix, while Ut, diagonalizes the neutrino Majorana mass term. 
Any 3x3 unitary matrix contains 3 moduli and 6 phases and can be written as [138] (see 
also [139]): 

U = e''^PUQ, (17) 

where P = diag(l, e**^, e*"^) and Q = diag(l, e^^, e''") are diagonal phase matrices having 2 
phases each, and f/ is a unitary "CKM-like" matrix containing 1 phase and 3 angles. The 
charged lepton Dirac mass term, M^, is diagonalized by a bi-unitary transformation: 

Me = Ul Mf f/t ^ (18) 

very pessimistic conclusion about the prospects to establish CP violation in the lepton sector due 
to Majorana CP violating phases from a measurement of {m)eff and, e.g., of mo, was reached in [99]. 



32 



where Ul^r are 3x3 unitary matrices and M^"^^ is the diagonal matrix containing the 
masses of the charged leptons. Casting Ul^r in the form (fT7|) . i.e., Ul,r = e**^-^ Pl,r Ul,rQl,r, 
we find 

Me = e^(*--*«) Pl Mf^^ Q^j, Uj, Pj, . (19) 

The term Pl M^^"^^ Q^j contains only 2 relative phases, which can be associated with the 
right-handed charged lepton fields. The three independent phases in e^^'^'^"'^^^ Ql can 
be absorbed by a redefinition of the left-handed charged lepton fields. Therefore, Ue is 
effectively given by Ul and contains three angles and one phase. 
The neutrino mass matrix M^, is diagonalized via 

= Ul Mf''^ ul . (20) 

The unitary matrix Ui, can be written in the form ()17|) . It is not possible to absorb phases 
in the neutrino fields since the neutrino mass term is of Majorana type [12, 114]. Thus, 

U = ul U, = e^*^ ill P, U, Q, . (21) 

The common phase has no physical meaning and we will ignore it. Consequently, in 
the most general case, the elements of U given by Eq. ()21|) are expressed in terms of six 
real parameters and six phases in Ui and Ui, . Only six combinations of those — the three 
angles and the three phases of U, are observable, in principle, at low energies. Note that 
the two phases in Q^, are "Majorana-like" , since they will not appear in the probabilities 
describing the flavor neutrino oscillations [12,114]. Note also that if Ug = 1, the phases 
in the matrix P,^ can be eliminated by a redefinition of the charged lepton fields. 
If one assumes that, e.g., Uu is bimaximal. 
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which permits a rather simple explanation of the smallness of sin ^13 and the deviation of 
^0 from 7r/4, then U^, is real . In this case the three angles and the Dirac phase in the 
neutrino mixing matrix U will depend in a complicated manner on the three angles and 
the phase in Ue and on the two phases in P^. The two Majorana phases will depend in 
addition on the parameters in Q^. See [102] for details. 

It should be emphasized that the form of U given in Eq. (j2ip is the most general one. 
A specific model in the framework of which Eq. (PT|) is obtained might imply symmetries 
or textures both in me and M^, which will reduce the number of independent parameters 
in ul and/or U^. 

In the scheme with three massive Majorana neutrinos under discussion there exist three 
rephasing invariants related to the three CP violating phases in U, 6 and 0i_2 [140-144]. 
The first is the standard Dirac one Jcp [140], associated with the Dirac phase S: 

JcP = lm{U,,U,2U:,U;,} . (23) 

It determines the magnitude of CP violation effects in neutrino oscillations [141]. Let us 
note that if L^^ = 1 and U^, is a real matrix, one has Jcp = 0- 
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The two additional invariants, 5*1 and 5*2, whose existence is related to the Majorana 
nature of massive neutrinos, i.e., to the phases a and /5, can be chosen as [142, 144] (see 
also [119])^: 

= Im {C/ei Ul^} , -52 = Im {U,-, \Jl^} . (24) 

If 5"! 7^ and/or ^2 7^ 0, CP is not conserved due to the Majorana phases 0i and/or 
02- The effective Majorana mass in /3/3oiy-decay, (m)e//, depends, in general, on 5"! and 
5*2 [119] and not on Jcp- Let us note, however, even if 5'i^2 = (which can take place if, 
e.g., If/esI = 0), the two Majorana phases </)i^2 can still be a source of CP non- conservation 
in the Icpton sector provided Im {f/ei 7^ and Im {f/^2 f^^s} 7^ [144]. 

Let us denote the phase in [/^ by t/). We will include it in Ui in the same way this is 
done for the phase 5 in the standard parametrization of the U matrix. If we write — 
diag(l, e"^, e*'^) and Qjj = diag(l, e*'', e'*^), the Dirac phase 8 in [/, which has observable 
consequences in neutrino oscillation experiments, is determined only by the "Dirac phase" 
in and the phases ip, and uj. The Majorana phases in U receive contributions also 
from the two remaining phases p and a. Allowing the phases 5 and 0i^2 to vary between 
and 27r, permits to constrain (without loss of generality) the mixing angles in C/g, 6ij, 
to he between and 7r/2. 

There arc interesting specific cases in which there are direct relations between all 3 
CP violating phases in the U matrix [102]. 



^We assume that the fields of massive Majorana neutrinos satisfy Majorana conditions which do not 
contain phase factors. 
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4 Testing Seesaw Models 



Although it is far from clear that the seesaw mechanism [67-70] is responsible for neutrino 
masses, most physicists consider that it is by far the most elegant mechanism. It fits very 
well into the big picture of other areas of particle physics such as super symmetry, grand 
unification etc. It is therefore important to discuss how we can test the seesaw models. 
Evidently, such tests are indirect, since the right-handed electroweak-singlet neutrinos are 
much too heavy to be produced at colliders. In this section, we consider two aspects of 
the seesaw mechanism: (i) indirect signals of seesaw mechanism in lepton flavor violating 
processes, which can be measured in near future; (ii) leptogenesis, which can give us a 
hint as to the CP violating phases in the lepton sector as well as perhaps the spectrum of 
the RH neutrinos. The presence of CP violating phases needed for leptogenesis (see Sec. 
14. 2j) in turn can lead to CP violating low energy observables in the seesaw models. We 
explore these probes of the seesaw mechanism in this section. 

4.1 Lepton Flavor Violation and Lepton Electric Dipole Mo- 
ments 

Neutrino oscillation experiments have revealed that the violation of flavor symmetry is 
much greater in the lepton sector than in the quark sector. We will discuss how this flavor 
violation manifests itself via the seesaw mechanism in other observable quantities where 
lepton flavor and/or CP are violated. As we are going to discuss, among the laboratory 
observables particularly interesting are Lepton Flavor Violating (LFV) decays, like — 67 
and r — * fij, and Lepton Electric Dipole Moments (LEDM), like d^. and d^. 





BR{fi 67) 


BR(r ^ /i7) 


de e cm 


d^ e cm 


present 145 


< 1.2 10-'' 


< 1.1 10-*^ 


< 1.5 10-^'' 


< 10-^8 


planned [146] 


< 10-14 


< 10-8 


^ -^Q-29(-32) 


^ lQ-2^-26) 



Table 5: Present status and future prospects for LFV decays and LEDM. 



Searches for LFV decays and for LEDM are experimentally very promising, since the 
present upper bounds could be strengthened by many orders of magnitude, as summarized 
in Table El Their impact on theory is also very promising: finding LFV and LEDM means 
discovering new low energy physics beyond the SM supplemented with the seesaw [61,147], 
like e.g. supersymmetry. This can be easily understood by identifying the SM with 
the operators of dimension c? < 4 of a low energy effective theory valid up to a cutoff 
A. Flavor and CP are accidentally conserved in the leptonic sector of the SM, hence 
there is no room for neutrino oscillations nor for LFV decays and LEDM. Their possible 
sources have to be found among the operators of d > 4: neutrino masses arise from 
the d = 5 operator i>^C-'h'{H^)'^ /A, while LFV decays and LEDM from the d = 6 
operator la^^{l + '-y^)lF^y{H^) / . In the seesaw lepton flavour is no longer conserved 
but A ~ Mr ^ 10^^ GeV and, as a consequence, the d = Q operator above is so strongly 
suppressed that its effects are negligibly small [61,147]. However, if additional physics is 
present at smaller mass scales and if this additional physics violates LF and/or CP, the 
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suppression is milder and LFV decays and LEDM could be raised up to the experimentally 
interesting range. 

This enhancement due to new low-energy physics is precisely what happens in low- 
energy supersymmetry ® where, due to loops with sleptons and gauginos, the d = & 
operator is suppressed by powers of msusY- The experimental limits on LFV decays and 
LEDM then imply such severe constraints [150] on the amount of LP and CP violations in 
slepton masses (defined in the basis where charged fermions are diagonal), that one would 
expect LP and CP to be exact symmetries of the supersymmetry breaking mass terms 
defined at the appropriate cutoff scale (the Planck scale for supergravity, the messenger 
mass for gauge mediation, etc). It is important to stress that, if below this scale there are 
LP and CP violating Yukawa interactions, in the running down to msusY they nevertheless 
induce a small amount of LP and CP violations in slepton masses. 

It is well known that this is the case for the seesaw interactions of the right-handed 
neutrinos [151] and/or the GUT interactions of the heavy triplets [152]. Remarkably 
enough, this radiative contribution to the LPV decays and LEDM, which essentially 
depends on the supersymmetric spectrum and on the pattern of the Yukawa interactions, 
might be close to or even exceed the present or planned experimental limits. Clearly, this 
has an impact on seesaw models, possibly embedded also in a GUT framework. In the 
following we will discuss separately the case of type I and type II seesaw. 

4.1.1 Type I Seesaw: LFV 

Por type I seesaw, in the low energy basis where charged leptons are diagonal, the ij mass 
term of L-sleptons, rn?ij , is the relevant one in the decay ^i — > £^7. Assuming for the sake 
of simplicity the mSUGRA [153] spectrum at Mpi, one obtains at the leading log [151] 
(see also [154]): 

= ^(3m^ + 2Ag)a, , Q, = (^^)^^ (^^).*^ 1^ ^ ' (25) 

where Yy = / w^k and uiq and Aq are the universal scalar masses and trihnear couplings 
at Mpi, respectively, and we have chosen the basis where Mr is diagonal. Por the full KG 
results, see [155]. The seesaw model dependence thus resides in |Cjj|, and an experimental 
limit on BR(£j ij'y) corresponds to an upper bound on \Cij\. Por /x — 67 and t ^ fij 
this bound [156] is shown in Pig. 0] as a function of the lightest charged slepton and 
gaugino masses. 

It has been shown that many seesaw models predict |C^e| and/or \Crfj.\ close to the 
experimentally accessible range [157, 158] and, in particular, this might be the case for 
models based on U{1) flavor symmetries [159]. To reduce the uncertainty due to the 
supersymmetric spectrum, it is interesting to exploit the correlation between LPV decays 
and muon g — 2 [160] or neutralino dark matter [161]. 

Planned searches could also help in discriminating between categories of seesaw models 
[156]. To give some hints on the latter issue consider, e.g., hierarchical eigenvalues of Y^. 
The different N'^ thresholds can then be neglected and in first approximation one has 

*It can also happen in quite different theories, such as extended technicolor [148,149]. 
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Figure 4: Upper bounds on |Ct/^|, |C^e| in the plane (Mi,me^), respectively the bino 
and right slepton masses expressed in GeV. We choose as reference values the planned 
experimental sensitivities BR{r — > /i7) < 10~^, BR{fj, — > ej) < 10~^^, with tan/3 = 20, 
reminding how to adapt the numbers in the plots for different values of these parameters. 
We also assume mSUGRA with mg = Aq and fi from radiative electroweak symmetry 
breaking. The plots are adapted from those in the last Ref. of [156]. 



\Cij\ ~ I^LisI l^jsbs log(Mpi/M3), where Vl is the lepton analog of the CKM mixing 
matrix. In S'O(10)-inspired models = yt ^ 1, and the model dependence essentially 
resides in the magnitude of iV^jsl | Vl^sI, namely on the amount of LF violation present in 
the left-mixings of Yy. Under the above assumptions |Ct-^| = 0(10 x |Vl32|). If at high 
energy LF is strongly violated in the r — /i sector (models with 'lopsided' as the one 
studied in [162]) planned searches for r — > could be successful for a significant region 
of the supersymmetric parameter space. If this violation is on the contrary tiny like in 
the quark sector — in which case the large LF violation observed at low energy purely 
arises from a magnification effect of the seesaw [163] — r — > /i7 would not be observed. 
Progress in the experimental sensitivity to the latter decay would thus offer precious 
information. The prediction for yU — > 67, linked to the product | ^2,23! I^Lial, is more model 
dependent but, on the other hand, the present experimental bound is already very severe. 
For instance, simple U{1) flavor symmetries, those with all lepton charges of the same 
sign, predict |C^e| = 0{1Q x Amg/Am^) and for LMA are already in crisis [156,159]. 
Since the present limit corresponds to a significantly smaller degree of LFV at high energy, 
this means that a much richer flavor symmetry has to be at work. Notice also that in the 
future we could test |C^e| up to the CKM-level [164]; in fact, if y^ = 0{1) and Vl ~ Vckm, 
then Cue = O{10~^) which is well inside Fig. El 



37 



4.1.2 Type I Seesaw: EDM 



Let us now discuss the consequences of type I seesaw models for lepton EDM. It is well 
known that in the simplest supersymmetric models (with or without neutrino mass) the 
dipole moments of electrons and muons obey a simple scaling law de/d^ ~ me/m^. Given 
the present bound on d^, this implies d^ < 10~^^ e cm, which is at the level of the best 
experimental prospects. 

Things can change in seesaw models due to the fact that interactions involving right- 
handed neutrinos via radiative corrections can affect the scaling law. In type I seesaw 
with degenerate N'^, the radiative contributions to de and still preserve the scaling 
law. However, with hierarchical N'^ this proportionality is broken due to threshold effects 
arising from both the flavor conserving A-term contribution [165] and by the FV double- 
insertion contributions [166,167], which dominate for tan/3 > 10. Nevertheless, if only 
the type I seesaw radiative contributions are taken into account, d^ and d^ turn out 
to be barely at hand of future experimental searches (but can be for very particular 
textures [165]). 

Discovering LEDM would then suggest the presence of additional particles and inter- 
actions beyond those of the supersymmetric type I seesaw. The heavy color triplets of 
GUT theories are excellent candidates for this [152,166,168]. In particular, it turns out 
that the limits on dg are competitive with those on proton lifetime in constraining the pat- 
tern of GUT theories where heavy triplets and right-handed neutrinos are simultaneously 
present [169]. 

4.1.3 Type II Seesaw 

We will now consider a class of models where the right-handed neutrino mass arises from 
a renormalizable coupling of the form f'NNAji, where is a right-handed neutrino, / is 
a coupling constant and is a Higgs field whose vacuum expectation value (vev) gives 
mass to the right-handed neutrino. This is a natural feature of models with asymptotic 
parity conservation, such as those based on SU{2)l x SU{2)r x U{1)b-l or any higher 
gauge group such as 50(10), where the A^j field is part of an SU{2)r triplet field. Parity 
invariance then implies that we also have an fui/Ai coupling term as a parity partner of 
the NNAfi^ coupling. In this class of theories, whenever A/j acquires a vev, so does A^ 
and they are related by the formula (A/,) = vl = where f„k is the weak scale, vr is 
the A^^ vev and 7 is a coupling constant in the Higgs potential. The Ax, vev contributes 
a separate seesaw suppressed Majorana mass to neutrinos leading to the type II seesaw 
formula (see Eq. (@))) [73]. In the case where right-handed Majorana masses are heavy 
enough, the second term in the type II seesaw formula can be negligible, and the first 
term. Ml = fvi, is dominant. We will call this type II seesaw. The type II seesaw gives 
rise to the most simple explanation of the neutrino sector and is phenomenologically very 
successful especially when we try to construct realistic models. 

The simplest model that can be imagined for type II seesaw has just MSSM and right- 
handed neutrinos below the GUT scale and hence there is no new symmetry breaking 
scale. The right-handed neutrino masses can have hierarchies and therefore get decoupled 
at different scales below the GUT scale. Due to the radiative corrections from the RGEs, 
the flavor- violating pieces present in Y^, and / get transmitted to the flavor universal scalar 
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masses and produce lepton flavor violation. The / term contribution is the additional 
contribution that is typical of type II seesaw models [170]: 



dY^jdt 



dYe/dt 




(26) 



drv?ij^/dt 



Here, Tn\i^ represents the soft left-handed scalar masses. The flavor non-diagonal pieces 
generate the lepton flavor violation. In the type II seesaw the structure / coupling gen- 
erates the neutrino mixing parameters. Both / and are determined by the particular 
model which explains the quark and lepton masses. 

In order to calculate the BRs of ^ 67, r — > yU7 and the electric dipole moments 
for the electron and muon, we use the minimal SUGRA universal boundary conditions 
at the GUT scale. The unifying framework of S'O(IO) has been chosen and the values of 
quark masses and the CKM CP violation are satisfied. The values of the universal scalar 
mass mo, universal gaugino mass mi/2, universal trilinear term Aq, tan/? and the sign of 
/i as free parameters determine the final result. The assumption of universality allows us 
to probe the flavor violation originating from the neutrino sector. We also assume that 
there is no phase associated with the SUSY breaking. The Yukawa and/or the Majorana 
couplings are responsible for CP violation in these models. 

The mSUGRA parameter space is constrained by the experimental lower limit on mu-, 
the measurements of 6 — >■ S7 and the recent results on dark matter relic density [33]. For 
low tan/3, the parameter space has a lower bound on mi/2 stemming from the light Higgs 
mass bound of mh > 114 GeV. For larger tan (3 the lower bound on mi/2 is obtained by the 
CLEO constraint on BR(6 — *■ 57). The lightest neutralino is the dark matter candidate 
in this model and we satisfy the 2a range of the recent relic density constraint ^cdm = 
0.1126]'^o;QQg [33] in the parameter space. The allowed parameter space of mSUGRA mostly 
reduces to the neutralino-stau co-annihilation region for mo, mi/2 < 1000 GeV and when 
we satisfy the relic density constraint, mo gets determined within a very narrow band. 
For example, mo varies between 60-100 GeV for Ao = line in the graph. In Figs. EHHl 
we show BR[/i — >• 67] and BR[r — >• /i7] as a function of mi/2 for different values of Aq. 
We find that the BR is large in most of the parameter space and can be observable. In 
addition, BR[r /i7] can also be observable in the near future. The figures demonstrate 
that lepton flavor violation typically increases with increasing tan (3. 

The electron EDM is plotted in Fig. [7| We find that the maximum value of EDM is 
~ 10~^^ e cm. The muon EDM is shown in Fig.|Hland the maximum value shown is about 
10~^^ e cm. The scaling is broken in this model. We do not assume any new CP phases in 
SUSY parameters, hence all CP phases arise from the Yukawa and Majorana couplings. 

It is clear that if the seesaw mechanism eventually turns out to be the explanation 
of small neutrino mass, the crucial question becomes whether it is of type I or type II. 
One may then use lepton flavor violation as a way to discriminate between these two 
possibilities. 
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Figure 5: BR[/i 67] is plotted as a function of mi/2 for different values Aq and 
tan/3 = 10, 40 and 50 in pure type II seesaw. 
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Figure 7: The electron EDM is plotted as a function of mi/2 for different values Aq and 
tan j3 = 40 and 50 in pure type II seesaw. 
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Figure 8: The muon EDM is plotted as a function of mi/2 for tan/5 = 10, 40 and 50 in 
pure type II seesaw. 
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4.2 Leptogenesis in the Type I Seesaw 



The origin of matter is a fundamental puzzle of cosmology and particle physics. The 
seesaw provides many mechanisms to generate this excess; we discuss what we can learn 
about neutrino physics, as well as the pattern of right-handed neutrino masses, from the 
observed baryon asymmetry. 

Three ingredients are required to generate the observed Baryon Asymmetry of the 
Universe [171]: baryon number violation, C and CP violation and some out-of-thermal 
equilibrium dynamics. The seesaw model [66-70], which was introduced to give small 
neutrino masses, naturally satisfies these requirements, producing the baryon asymmetry 
by the "leptogenesis" mechanism [172]. It is interesting to investigate the relation be- 
tween the requirements of successful leptogenesis and the observable neutrino mass and 
mixing matrices. In particular, does the CP violation that could be observed in neutrino 
oscillations bear any relation to leptogenesis? Do the Majorana phases that appear in 
neutrinoless double beta decay experiments do so? 

The next subsection reviews the thermal leptogenesis scenario, focusing on the Type I 
seesaw, with hierarchical RH neutrinos (Mi ^ M2,3/10). The relation with light neutrino 
parameters in type I seesaw models with three generations is discussed in Subsection l4.2.21 
The situation in type II seesaw is given in Section 14. 3t and the case of quasi-degenerate 
Nji masses is discussed in subsection 14.51 



4.2.1 Thermal Leptogenesis 

The idea of leptogenesis is to use the lepton number violation of the Nfi Majorana masses, 
in conjunction with the B + L violation contained in the Standard Model [173], to gener- 
ate the baryon asymmetry. The most cosmology-independent implementation is "thermal 
leptogenesis" [172, 174, 175, 182, 183] which will be reviewed in the following paragraph. 
Other leptogenesis scenarios, where the initial number density of (s)neutrinos is produced 
non-thermally (by inflaton decay, scalar field dynamics,. . .) depend on additional param- 
eters of the cosmological model. 

If the temperature Trh of the thermal bath after inflation is ^ Mjn, the lightest 
will be produced by scattering. If the Nr subsequently decay out of equilibrium, a CP 
asymmetry ei in the decay produces a net asymmetry of Standard Model leptons. This 
asymmetry is partially transformed into a baryon asymmetry by the non-perturbative 
B + L violation [184]. Thermal leptogenesis has been studied in detail [174, 175, 177, 182, 
183, 185]; the baryon to entropy ratio produced is 

71 

Ys-Ct^-er , (27) 

s 

where k < 1 is an efficiency factor to be discussed in a moment, n/s ^ 10~^ is the ratio of 
the equilibrium number density to the entropy density, and ei is the CP asymmetry 
in the Njn decay. C ~ 1/3 tells what fraction of the produced lepton asymmetry is 
reprocessed into baryons by the B + L violating processes. depends largely on three 
parameters: the Nm mass Mjn, its decay rate Fi, and the CP asymmetry ei in the 

decay. The decay rate of N^j can be conveniently parametrized as F^ = ^g|^" = 
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Figure 9: The efficiency parameter k as function of (mi, M/ji/Trh) for the Type I Seesaw 
with hierarchical Nr, in the SM and MSSM. In this plot Mm = 10^° GeV; the plot would 
only change slightly for Mjn <^ 10^'' GeV. 

171 ■ ? 

g-^^ , where rrij is often of order of the elements of the mass matrix, although it is a 

wk 

rescaled Nr decay rate. The requisite CP violating decay asymmetry ei is caused by the 
interference of the tree level contribution and the one-loop corrections in the decay rate 
of the heavy Majorana neutrinos. For hierarchical neutrinos it is given by: 

^ T{Ni i+) - T{Ni ^ r ) 

" r(iVi^$-r) + r(iVi-.$+r) ^^^^ 

where $ and ^ indicate the Higgs field and the charged leptons, respectively. 

Eq. (P7j) can be of the order of the observed ~ 3 x 10~^^ when the following 
conditions are satisfied: 

1. Mri should be ^ Trh^. This temperature is unknown, but bounded above in 
certain scenarios {e.g. Trh ~ 10^ GeV due to gravitino overproduction in some 
supersymmetric models, see 14. 2. 6|) . This can be seen in Figure El where the efficiency 
factor K falls off rapidly for Mri > Trh. 

2. The Nri decay rate oc fhi should sit in a narrow range. To be precise, fhi must 
be large enough to produce an approximately thermal number density of A^ri, and 
small enough that the Nri lifetime is of order the age of the Universe at T ~ Mri 
(the out-of-equilibrium decay condition). These two constraints are encoded in the 
efficiency factor k, plotted in Figure IHl 

3. ei must be ^ 10^^. 

In Figure ^1 is plotted the baryon asymmetry, produced by thermal leptogenesis, as a 
function of Mri and mi, for Trh ^ Mri, and ei = 10~^. To reproduce the observations, 
Mri and mi must be inside the three neighboring (blue) lines. 

^In the so-called 'strong washout' regime, Trh can be an order of magnitude smaller than Mi [185]. 
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Figure 10: Contour plot, from [186], of the baryon to photon ratio produced in thermal 
leptogenesis, as a function of Mpti and rhi. The decay asymmetry ei was taken to be 
10~®. The three (blue) close-together lines are the observed asymmetry. The horizontal 
contours, for small rhi assume a thermal Nji abundance as initial condition. 

4.2.2 Parametrizing the type I seesaw 

Twenty-one parameters are required to determine the three generation lepton sector of 
the type I seesaw model. This includes the charged lepton masses, and a mixing matrix 
(with three complex parameters, e.g. the MNSP matrix) in the left-handed lepton sector. 
The remaining 9 real numbers and 3 phases can be chosen in various ways: 

1. " top-down" — input the Nfi sector: the eigenvalues of the mass matrix and of 
Y^, and a matrix transforming between the bases where these matrices are diagonal 
[138] (see also Refs. [187,188]). 

2. " bottom-up" — input the sector: the eigenvalues of the mass matrix and of 
Y,^, and a matrix transforming between the bases where these matrices are diagonal. 

3. "intermediate" — the Casas- Ibarra parametrization [158]: the and eigen- 
values, and a complex orthogonal matrix R which transforms between these two 
bases. 

To relate the RH parameters relevant for leptogenesis to the LH ones, many of which are 
accessible at low energy, it is useful to consider the first and second parametrization. 

4.2.3 Implications for CP conserving observables 

The second requirement (i.e., the range of mi) sets an upper bound on the mass scale 
of light neutrinos. The scaled decay rate rhi is usually ~ m2,ms; for hierarchical light 
neutrinos, it naturally sits in the desired range. One can show [189,190] that mi < mi, 
so that mi ^ 0.15 eV [183,185,191,192] is required for thermal leptogenesis in the type I 
seesaw, with hierarchical N^. This is shown in Figure ITTl 

In type I seesaw models with hierarchical Nji, the third condition (ei ^ 10~^) im- 
poses Mri ^ 10® GeV, because ei < 3Mii;i(m3 — mi) / {Snv'^^^) in most of parameter 
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Figure 11: Upper bound on the light neutrino mass scale, assuming hierarchical A^r, taken 
from [192]. The plot shows the measured baryon asymmetry (horizontal line) compared 
with the maximal leptogenesis value as function of the heaviest neutrino mass m^. 



space [190, 192]. For three Nji, the value of Mm has little implication on low energy 
neutrino observables. If ei is maximal — that is, M/ji close to its lower bound, — this 
sets one constraint on the 21 parameters of the type I seesaw. This has no observable 
consequences among Standard Model particles, because at most 12 masses, angles and 
phases are measurable, and ei can be maximized by choice of the nine other parameters. 
The situation is more promising [138, 193] in SUSY models with universal soft terms, 
where some of the 9 additional parameters can contribute to slepton RGEs and thereby 
to lepton flavor violating branching ratios, as discussed in Section 14.11 



4.2.4 Relations between Leptogenesis and leptonic CP violation 

The leptogenesis parameter ei is a C/P asymmetry, suggesting a possible correlation with 
CP violation in u oscillations (the phase 6), or to the low energy Majorana phases (0i,2)- 
Let us assume that ei is large enough — so thermal leptogenesis works — and concentrate 
on the implications for low-energy CP violation. 

The first thing that must always be said, in discussing potential connections between 
phases in the MNSP lepton mixing matrix and leptogenesis, is that there is no linear 
relation: leptogenesis can work when there is no C/P in MNSP, and measuring low energy 
leptonic phases does not imply that there is CP violation available for leptogenesis. This 
was clearly and elegantly shown by Branco, Morozumi, Nobre and Rebelo in [194]. The 
problem is that six phases are included in the general three neutrino seesaw scenario — 
it would be astonishing if the C/P parameter we are interested in (ei) is proportional to 
the low energy phases (5, 0i,2)- 

Nonetheless, some sort of relation between ei and the low energy phases would be 
interesting — so what can we say? Needing more inputs that the data provides is a familiar 
problem for extensions of the SM. The usual solutions are to scan over unknowns, or to 
fix them. But this is subtle: any relation depends on the choice of "independent" phases. 
For instance, if ei and S are chosen as inputs, then it follows that they are unrelated 
A choice of parametrization is not obvious: the C/P parameter ei is a function of 
phases, masses, mixing angles, whereas the observable phases are those of U. 

A useful step is to write ei oc a Jarlskog invariant, which can be done for thermal 

-'^"This arises in the Casas-Ibarra seesaw parametrization. 
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leptogenesis with hierarchical N^f [195]: 

ei (X Q{Tt[MIM,MI{Y,YJ)-'M,{Y:yJ)-']}. (29) 

The advantage is that Jarlskog invariants can be evaluated in any basis/parametrization. 
It is easy to see, evaluating Eq. in the ul mass eigenstate basis, that the C/P 
for leptogenesis is controlled by a matrix W which transforms between the bases were 
Yi, and Mj^ are diagonal. This matrix is unobservable, verifying the no-go theorem of 
[194]. However, in many popular/common Yukawa texture models, where Y^, and 
are almost simultaneously diagonalizable, W ^ U and the MNSP phases are relevant for 
thermal leptogenesis. For S13 larger than the mixing angles between diagonal Yy and M^, 
ei cx sin 2(01 — (1)2 + 5). 

4.2.5 Model dependent approaches 

Within specific models interesting links between the phase relevant for leptogenesis and 
the phase 5 measurable in neutrino oscillation experiments have been made. The precise 
link depends on how many "texture" zeroes are assumed to be present in the neutrino 
Dirac mass matrix. For example within the class of two right-handed neutrino models, if 
two texture zeroes are assumed then there is a direct link between 5 and the leptogenesis 
phase, with the sign of 5 being predicted from the fact that we are made of matter 
rather than antimatter [196,197]. Two right-handed neutrino models can be obtained as 
a limiting case of sequential dominance models, and in such models if only the physically 
motivated texture zero in the 11 entry of the Dirac mass matrix is assumed, then the link 
is more indirect [198]. Other approaches which give rise to a link between leptogenesis 
and CP violation include GUT models [199], or textures [117,138,187,200] or left-right 
symmetric models, to be discussed in the next Subsection. On the other hand, if the 
charged lepton sector contributes significantly to the lepton mixing 6'i3 and therefore also 
to 5, such links may be spoiled [201]. 

4.2.6 Leptogenesis in supersymmetric scenarios 

In supersymmetric scenarios, the history of the early universe is subject to various con- 
straints. Many of them are associated to the gravitino problem [202,203]. In short, 
unstable gravitinos are notoriously in conflict with nucleosynthesis (see [204] for a recent 
analysis), while stable ones may 'overdose' the universe. If the gravitino is very light 
^ keV [205]) or very heavy (7/13/2 ^ lOTeV), these bounds disappear, and thermal 
leptogenesis works (see, e.g., [206]). For all other masses, nucleosynthesis or 'overclosure' 
constraints translate into bounds on the gravitino abundance at T ~ 1 MeV or today, 
respectively. Assuming that gravitinos are not produced by infiaton decays (see [207]), 
this gravitino abundance is linear in the reheat temperature [208]. 

Unstable gravitinos with masses below ~ 10 TeV lead to severe constraints on the 
reheat temperature Trh [204] which are in conflict with thermal leptogenesis where Trh ~ 
10^ GeV. There are, however, various alternative leptogenesis scenarios such as non- 
thermal leptogenesis [209] where the heavy neutrinos are directly produced by infiaton 
decays, or mechanisms using the superpartner of the neutrino, the sneutrino [210], or 
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sneutrino oscillations (see Subsec. 18 .31) . These scenarios can be consistent with unstable 
gravitinos. 

Stable gravitinos, on the other hand, may evade the constraints from nucleosynthesis 
provided that the decays of the next-to-lightest superpartner into the gravitino are harm- 
less [211]. However, the 'overclosure' constraint leads to Trh ~ (10^ — 10^°) GeV. Such an 
upper bound on the reheat temperature is suggested independently by string-theoretical 
arguments where Tjih ^ ■^m3/2 Mp [212]. Stable gravitinos are thus (marginally) con- 
sistent with thermal leptogenesis, and provide a natural dark matter candidate [213]. It 
is clear that the neutrino mass bound, as discussed in 14.2.31 becomes much tighter now 
since fhi ~ 10~^eV (cf. Fig. irUj) and mi ^ rhi [189,190]. This scenario therefore predicts 
hierarchical light neutrinos as well as gravitino cold dark matter. These predictions will 
be tested in neutrino experiments and at future colliders [214]. 

4.3 Leptogenesis and type II seesaw mechanism 

In type II seesaw scenarios the neutrino mass matrix Mi^ reads 

Mj^ = Ml- M^' {Mj^f = Ml + Mi , (30) 

where we divided the mass matrix in the conventional type I part Ml = —Mj^ M^^ [Mj^Y' 
and the part characteristic for type II, Ml. A type II seesaw term can, e.g., be present in 
5*0(10) models, in which the B — L symmetry is broken by a 126 Higgs fields. Depending 
on the parameters of the model, either Ml or Ml can be the dominating source of Ml^ . As 
mentioned above, in case of the conventional type I seesaw mechanism with three families 
of light and heavy Majorana neutrinos, there are six phases. As already discussed, in this 
case there is in general no relation between the PMNS phase and leptogenesis phase. 

For the type II case the phase counting gives the result of 12 independent CP phases 
and there is no connection between low and high energy CP violation either. The number 
of CP phases can be obtained by going to a basis in which both Ml and Mr are real and 
diagonal. Any CP violation will then stem from the matrices M^f and M^Mj (with M^ 
being the charged lepton mass matrix). Those two matrices posses in total 9 + 3 = 12 
phases. 

The term Ml is induced by a SU{2)l Higgs triplet, whose neutral component acquires 
a vev vl oc "W^k/MA^, where Ma^ is the mass of the triplet and f^k the weak scale. 
Consequently, the triplet contribution to the neutrino mass matrix is 

Ml = VLfL, (31) 

with Jl a symmetric 3x3 coupling matrix. The magnitude of the contribution of Al 
to Ml^ is thus characterized by its vev vl- In left-right symmetric theories the left-right 
symmetry necessarily implies the presence of a SU{2)fi triplet, whose coupling matrix 
is given by /j? = //, = / and its vev is given by vr, where vl vr = j v'^^^ with 7 a 
model dependent factor of order one. The right-handed Majorana neutrino mass matrix 
is thus given by Mr = vr f = vr/vl Ml- Before acquiring its vev, the presence of the 
doubly charged Higgs and the coupling of the Higgs triplet with the doublet introduces the 
possibility of additional diagrams capable of generating a lepton asymmetry. First, there 
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Figure 12: Scatter plots of the effective mass against tan6'i3 and Amg against Jcp for 
the normal (left) and inverted (right) neutrino mass spectrum. Taken from [220]. 



is the possibility that in the decay Ni ^ L a virtual Higgs triplet is exchanged in the 
one-loop diagram, which contributes to the decay asymmetry in the decay of the heavy 
Majorana neutrinos [215-218]. The corresponding term adds up to the conventional 
term ei whose properties were discussed in the previous Subsection. 

The second new diagram possible is given by the decay of the doubly charged Higgs 
triplet in two charged leptons. One-loop exchange of a heavy Majorana neutrino gives 
rise to the decay asymmetry [215] 

^ r(A^ r n - TjAl I I) 

~ t{Al I' /^) + r(A* ^11) ' ^ ' 

If Ml -C (Ml ^ Ma^) the decay of the Majorana neutrino (Higgs triplet) will govern 
the baryon asymmetry. Thus, depending on which term dominates M^^ , four different 
situations are possible [217]. The discussion of 3 of the cases has so far not been discussed 
in as much detail as the conventional leptogenesis in type I seesaw mechanisms. 

If Ml -C Ma^ and the conventional term M^ dominates M^^ , we recover the usual 
seesaw and leptogenesis mechanisms and the statements given in Sec. 14.2.11 applv. 

In situations in which Mi <^ Ma^, the heavy Majorana neutrinos display a hierarchical 
structure and Ml dominates M^' , it has been shown in [217,218] that one can rewrite the 
decay asymmetries such that depends on Ml and ei on M^. However, since matrices 
are involved, ei can still be the dominant contribution to the decay asymmetry, a situation 
which in the context of left-right symmetry has intensively been investigated in [219,220], 
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see also [221]. Calculating in this framework the baryon asymmetry in terms of light 
neutrino parameters (a bottom-up approach) leads typically to a main dependence on the 
Majorana phases in the PMNS matrix. If is given by the up-quark mass matrix and 
the light neutrinos display a normal hierarchal spectrum, one of the low energy Majorana 
phases has to be very close to zero or 7r/2 [219]. For given by the down quark or 
charged lepton mass matrix one finds that the in general unknown mass spectrum of the 
heavy Majorana neutrinos is exactly given by the measurable mass spectrum of the light 
Majorana neutrinos. In case of a normal hierarchy, the Majorana phases should lie around 
7r/4 or 57r/4. Both values give comparable results for the rate of neutrinoless double beta 
decay. Thus, measuring neutrinoless double beta decay fully determines the neutrino 
mass matrix in this scenario. It is also possible to set limits on the lightest neutrino mass 
mi because the baryon asymmetry is proportional to mi. It should be larger than 10~^ 
eV in order to produce a sufficient baryon asymmetry [220]. For an inverted hierarchy 
of the neutrinos it turns out that rather sizable values of ^13 are required. Thus, sizable 
effects of CP violation in future long-baseline neutrino oscillation experiments are possible. 
The preferred value of the Majorana phase implies in addition a rather sizable rate of 
neutrinoless double beta decay. Furthermore, the lightest neutrino mass should be heavier 
than 10~^ eV. Figure ^1 shows for the normal and inverted hierarchy typical examples 
for the expected values of 6'i3, the effective mass and the CP violating parameter Jcp in 
neutrino oscillations. A similar example within a framework incorporating spontaneous 
CP violating is discussed in Section 15.2.41 

Finally, if neutrinos possess a quasi-degenerate mass spectrum, one of the Majorana 
phases is required to lie around tt or 7r/2. A measurement of neutrinoless double beta 
decay can resolve this ambiguity. 

The other possible scenarios have not been discussed in detail in the literature so far 
(see, e.g., [222]). General statements are however possible. If, e.g.. Mi -C and 
the term dominates the decay asymmetry, the limits on the light neutrino masses of 
order 0.1 eV (see Sec. 14.2.^ no longer apply [217], since the couplings responsible for the 
neutrino masses do not influence the wash-out processes. For hierarchical light neutrinos 
the upper bounds on the decay asymmetry ei and are identical. In case of quasi- 
degenerate neutrinos, however, the limit in case of type II seesaw is weaker by a factor 
of 2m'^/ ^m\ [218], where mo is the common neutrino mass scale. Along the same lines, 
the lower limit of order 10^ GeV on the lightest of the heavy Majorana neutrinos can be 
relaxed by roughly one order of magnitude [218], thereby making thermal leptogenesis 
less in conflict with the gravitino problem. 

Consider the case when Mi S> M/^^ and Ml dominates Ml^ . Then eA will produce 
the baryon asymmetry and again the limits on light neutrino masses do not apply. The 
same is true when Mi ^ and Ml is the main contribution to Ml^ . A smaller range 
of allowed parameters is expected in this case [217]. 

Therefore, given the fact that quasi-degenerate light neutrinos are hard to reconcile 
with standard thermal leptogenesis in type I seesaw models, our discussion implies that 
if we learn from future experiments that neutrinos are indeed quasi-degenerate, triplet 
induced leptogenesis represents a valid alternative (this was first noted in [223]). 

Also possible is — in inflationary scenarios — that the decay of the inflaton into light 
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particles together with interference of one-loop diagrams with exchanged SU{2) triplets 
and heavy Majorana neutrinos generates a lepton asymmetry. Various slepton decays in 
future coUiders are expected to be observable [224] . 

Up to now the discussion was constrained to the presence of only one triplet. If only 
one triplet is present, right-handed Majorana neutrinos are necessary to produce a de- 
cay asymmetry. Introducing two or more triplets allows self-energy diagrams which can 
produce a decay asymmetry without heavy Majorana neutrinos [225]. The possibility 
of lowering the triplet mass scale due to a resonance effect of close-in-mass triplets is 
possible [226] , giving the prospect of collider phenomenology. The presence of light and 
detectable Majorons is also possible. There are models implementing this kind of triplet 
self-energy scenarios with light left-handed neutrinos [226] and with quasi-degenerate 
ones [227]. The latter also predicts a stable proton due to R parity conservation. Intro- 
ducing the triplet induced neutrino mass matrix of the type II seesaw mechanism along 
the lines of [228], i.e., by a conjunction of flavor and permutation symmetries will typically 
include many additional Higgs fields. Rich phenomenology in form of rare charged lepton 
decays or charged lepton EDMs will be among the interesting consequences. 

To sum up, the type II seesaw mechanism displays the most general but more compli- 
cated framework of neutrino mass and leptogenesis. Nevertheless, richer phenomenology 
is expected, most of which remains to be explored. 

4.4 Dirac Leptogenesis 

Since the seesaw mechanism coupled with existing data on neutrino masses and mixings 
does not give complete information about the RH neutrino sector, one must consider lep- 
togenesis within various scenarios for RH neutrino masses that correctly explain neutrino 
observations. In this section, we discuss a possibility that the leptogenesis occurs with 
Dirac neutrinos. The conventional leptogenesis [172], which we call Majorana leptogene- 
sis for definiteness, was based on the fact that the standard model violates B + L [173], 
while the Majorana neutrinos violate L, and hence both B and L are violated. Therefore 
it is possible to create L from the decay of right-handed neutrino that is subsequently 
converted to B [184]. On the other hand, Dirac neutrinos conserve L and hence B — L 
is an exact symmetry. Therefore B — L stays vanishing throughout the evolution of the 
universe and it appears impossible to generate non- vanishing baryon asymmetry^^. 

Dirac leptogenesis overcomes this problem by the following simple observation [229]. 
Recall that the Dirac neutrinos have tiny Yukawa couplings, = Fy^wk, — 10~^^. If 
this is the only interaction of the right-handed neutrinos, thermalization is possible only by 
processes hke NL HW and they do not thermahze for T > g^Y^M-pi ~ 10 eV.^^ At this 
low temperature, obviously both H and W cannot be produced and the thermalization 
is further delayed until ~ M^, when neutrinos become non-relativistic. Therefore the 

^^An obvious exception is electroweak baryogenesis, where B — L = while B = L ^ after the 
electroweak phase transition. 

^^In contrast, the Yukawa couphng of the right-handed electron cr is large enough to equilibrate the 
eijs before sphaleronic processes switch off [230]. 
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Figure 13: The evolution of the lepton asymmetry in Dirac leptogenesis models. At the 
first stage, an asymmetry between the ordinary leptons and the right-handed neutrinos is 
created without lepton-number violation. Then the asymmetry in the ordinary leptons is 
partially converted to the baryon asymmetry. Finally, the right-handed neutrinos come 
in thermal equilibrium with the other leptons. The net baryon and lepton asymmetries 
remain while the overall B — L vanishes. 

number of left-handed and right-handed neutrinos are separately conserved practically up 
to now. We call them L and N , respectively, and the total lepton number is L + A^. The 
combination L + N — B is strictly conserved. 

Suppose the decay of a heavy particle produced an asymmetry L = —N ^ 0. The 
overall lepton number is conserved (see Fig. ITS|) . N is frozen down to Ty. On the other 
hand, the lepton asymmetry L is partially converted to the baryon asymmetry via the 
Standard Model anomaly. After the electroweak phase transition T ^ 250 GeV, the 
anomaly is no longer effective. Finally at T^,, L and equilibrate. In the end there is a 
baryon asymmetry i? = — (L + A^). 

The original paper [229] introduced new electroweak doublet scalar $ that has the 
same quantum numbers as the Higgs doublets and Yukawa couplings $LA^ and ^*LE. If 
there are two sets of them, there is CP violation and their decays can create the asymmetry 
L = —N 7^ 0. However, these doublets are there just for this purpose and have no other 
motivations. 

On the other hand, light Dirac neutrinos are natural in models where the neutrino 
Yukawa couplings are tied to the small supersymmetry breaking effects [231,232]. The 
Dirac neutrino mass is due to the effective operator 



where x is the hidden sector field which acquires a vacuum expectation value (x) — 
^3/2 + ^^^3/2Afpi and M is the heavy mass scale. The neutrino Yukawa coupling is 
Yy ~ 771^/2/ M, and is naturally small. The operator can be obtained by integrating out 
(two sets of) new doublets (p + (j)'^ that couple as W = 4>NHu + 4>'^Lx + Mcpcj)'^. The 
asymmetries L = —N 7^ are created by the decay of (p [233]. Then the origin of 
small neutrino mass and the origin of the lepton asymmetry are tied in the same way 
as the Majorana leptogenesis. Also concerning the gravitino problem (cf. I4.2.(ij] . Dirac 
and Majorana leptogenesis are on the same footing [233,234]. Such a scenario may be 
supported by the lack of neutrinoless double beta decay as well as the existence or right- 
handed sneutrino at LHC and Linear Collider. 




(33) 
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4.5 Resonant Leptogenesis 

The right handed neutrino sector of generic seesaw models is almost entirely unconstrained 
by existing data on neutrino masses and mixings. It is therefore necessary to consider all 
the various possibilities for the structure of the RH neutrino sector which are compatible 
with current experimental data. In this section, we consider the case where two or more 
right-handed neutrinos are nearly degenerate in mass. 

An important, further motivation for this possibility comes from the severe limits 
on the right handed neutrino sector that exist in models of thermal leptogenesis with 
hierarchical right handed neutrinos. In particular there exists a bound on the mass of 
the lightest right handed neutrino, Mjn < Trh, discussed in section 4.2.1. In supersym- 
metric theories with unstable gravitinos, this bound can be in conflict with the bound 
Trh ~ 10^ GeV, coming from nucleosynthesis considerations (see 4.2.6). This motivates 
us to consider scenarios where the scale of the right handed neutrino masses can be low- 
ered whilst still being compatible with thermal leptogenesis [175,176,239]. This may 
be achieved naturally in scenarios with nearly degenerate right handed neutrinos [175], 
in complete accordance with current neutrino data, and with the advantage that the fi- 
nal baryon asymmetry generated is independent of the initial lepton, baryon or heavy 
neutrino abundances [177]. 

If the mass difference between two heavy Majorana neutrinos happens to be much 
smaller than their masses, the self-energy (e-type) contribution to the leptonic asymmetry 
becomes larger than the corresponding (e'-type) contribution from vertex effects [174, 
175]. Resonant leptogenesis can occur when this mass difference of two heavy Majorana 
neutrinos is of the order of their decay widths, in which case the leptonic asymmetry 
could be even of order one [175,177]. As a result, one can maintain the RH neutrino 
masses around the GUT scale [242] or one can contemplate the possibility that the heavy 
neutrino mass scale relevant to thermal leptogenesis is significantly lower, for example 
in the TeV range [175]. This of course requires a different realization of the seesaw 
mechanism [90] but it can be in complete accordance with the current neutrino data [177]. 

The magnitude of the e-type CP violation occurring in the decay of a heavy Majorana 
neutrino iVj is given by [175], 



where F^i* is the tree level total decay width of Ni. It is apparent that the CP asymmetry 
will be enhanced, possibly to e ~ 1, provided 



l-p(o) 
2 

im(i;tn)2, 



rriN^ - niN^ ~ -F^^^ , 



{Yl, Yy)ii{Yu Yy 



1 . (35) 



3] 



It is important to note that Eq. (jH^ is only valid for the mixing of two heavy Majorana 
neutrinos. Its generalization to the three neutrino mixing case is more involved and is 
given in [177]. 
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Figure 14: Numerical estimates of the lepton to photon, and neutrino to photon ratios, r^x, 
and ?7Ari 2 as functions oi z = mN^/T for scenarios with mjv^ = 1 TeV. The model is based 
on the Froggatt-Nielsen mechanism and is completely consistent with all current neutrino 
data (see [177] for details). It naturally provides a degeneracy of I!!^ _ 1 = 9.2 X 10-i\ 

with i^T-factors and CP asymmetries of Ki = K2 = 6570 and 5^^ = = —0.003. ^ is a 
free parameter. The horizontal dotted line shows the value of t]l needed to produce the 
observed tib- The vertical dotted line corresponds to T = Tc = 200 GeV. 



Successful leptogenesis requires conditions out of thermal equilibrium. To quantify 
this, we introduce the parameter, K''^^ = r^./if(T = mjvj where H{T) is the Hubble 
parameter and F^y. is the total decay width of Ni into a lepton species / (/ = e,fi,T). In 
typical hierarchical leptogenesis scenarios K^^. is small, usually K^^. ~ 1. This constraint 
can be translated directly into an upper bound on the Yukawa couplings of the neutrinos 
which can be expressed in terms of effective light neutrino masses, rhi, 

m = "^^^^^ ^ lO-^i^^.eV, (36) 

where K]si^ = -^tv- However, resonant leptogenesis can be successful with values of K^^ 
larger than 1000 [177] (see Figure IT^ . This has implications for leptogenesis bounds on 
the absolute mass scale of the light neutrinos. If a large, ^ 0.2 eV, Majorana mass was seen 
in neutrinoless double beta decay, this could be naturally accommodated with resonant 
leptogenesis, whereas thermal leptogenesis models based on hierarchical heavy neutrinos 
would be strongly disfavored, as they naturally require smaller values of Kj^. [186,235]. 

Conditions close to thermal equilibrium (with large -R'atJ endow resonant leptogenesis 
models with another particularly attractive feature; the final baryon asymmetry generated 
is almost independent of the initial baryon, lepton or heavy neutrino abundances [177,179]. 
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For Kj^. > 1, and under the assumption that the neutrino Yukawa couphngs are the 
same for each lepton flavour, an order of magnitude estimate for the baryon to photon 
ratio may be obtained by [177] 

r 

Ni 

where 5]^^ is the leptonic CP asymmetry in the decay of Ni, K = Ki and 

Ki = ^ e-^™^--™^!)/™^! i^jv^. (38) 

N, 

It is apparent that if the CP-asymmetry is enhanced, for example through resonant effects, 
then K can be increased without an impact on the final baryon asymmetry. 

In resonant leptogenesis scenarios with neutrino Yukawa couplings that are not uni- 
versal for each lepton flavour, the effects of individual lepton flavours on the resultant 
baryon asymmetry may become very important [178,179]. This also applies to scenarios 
with mildly hierarchical RH neutrinos [179]. These effects may result in an enhancement 
of the baryon asymmetry predicted using ()37|1 by a factor as large as 10^ in some models 
of resonant leptogenesis. An order of magnitude estimate for the final baryon asymmetry, 
when the neutrino Yukawa couplings are not flavour- universal, may be obtained by [179] 

3 j^l 

r]B 10-2 X V y e-^'"^.-™^!)/'"^! 5L ^ , (39) 

1=1 Ni ' 

where (5jv^ is the CP asymmetry in the decay of Ni to leptons of flavour /. For a more 
precise computation of the resultant baryon asymmetry, a network of Boltzmann equations 
must be solved, one for each heavy Majorana neutrino species and one for each lepton 
flavour [179]. By adding a further Boltzmann equation for the baryon abundance, and 
including effects due to the rate of -B -|- L violating sphaleron transitions, it can be shown 
that successful resonant leptogenesis is possible with heavy Majorana neutrinos as light 
as the electroweak scale [179]. 

In particular, models of resonant r-leptogenesis [178], where a lepton asymmetry is 
generated predominantly in the r-family, can allow large Yukawa couplings between the 
e and n lepton families and some of the RH neutrinos. These couplings, in conjunction 
with the low RH neutrino scale, lead to a signiflcant amount of accessible phenomenology, 
such as potentially observable neutrinoless double beta decay, /i 67, fi — > eee and 
coherent n ^ e conversion in nuclei, and the possibility of the collider production of 
heavy Majorana neutrinos [178,179]. 

The conditions for resonant leptogenesis can be met in several ways. Models based on 
the Froggatt-Nielsen mechanism can naturally provide nearly degenerate heavy Majorana 
neutrinos satisfying Eq. (|35p and provide a light neutrino spectrum fulfllling all experimen- 
tal constraints. It can be shown that a model like this can produce the observed baryon 
asymmetry by solving the network of Boltzmann equations - including gauge mediated 
scattering effects [177] (see Figure IT^ . In this model the 'heavy' Majorana neutrinos can 
be as light as 1 TeV. 
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5*0(10) models with a "type III seesaw mechanism" naturally predict pairs of nearly 
degenerate heavy Majorana neutrinos suitable for resonant leptogenesis [236]. In addition, 
a model of neutrino mass from SUSY breaking has been shown to naturally lead to 
conditions suitable for resonant leptogenesis [238] . 

In the radiative leptogenesis mechanism [181], small mass differences arise through 
renormalization group corrections between RH neutrinos which are exactly degenerate in 
mass at some high scale. The leptonic CP asymmetry induced in this scenario is sufficient 
to produce the observed baryon asymmetry. 

In soft leptogenesis [180, 237] , soft SUSY breaking terms lead to small mass differences 
between sneutrinos. Resonant effects allow sneutrino decay to generate the required CP 
asymmetry. 

Several other mechanisms for leptogenesis where the right handed neutrinos can be 
at a TeV scale have been suggested [239] . Clearly for the seesaw mechanism to operate 
in such models, the Dirac mass must be constrained, e.g., by a leptonic global U(l) 
symmetry [240]. 
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5 The Heavy Majorana Mass Matrix 



5.1 General Considerations 

We have already seen from the Introduction that the type I seesaw mechanism requires 
the existence of right-handed neutrinos N^, and then the hght Majorana mass matrix is 
given as 

= -Mj^M^^Mf (40) 

where M^f is the Dirac neutrino matrix (to be thought of as perhaps similar to the quark 
and charged lepton mass matrices) and Mr is the heavy Majorana mass matrix. While 
the elements of must be at or below the electroweak scale, the characteristic scale 
of right-handed neutrino masses can and must be much higher. Having introduced right- 
handed neutrinos into the Standard Model for the purpose of accounting for light physical 
neutrino masses via the type I seesaw mechanism^^ it is clearly an important question to 
understand the mass spectrum and couplings of the right-handed neutrinos. Since their 
only couplings are their Yukawa couplings to Higgs and left-handed neutrino fields, it will 
clearly not be an easy task to answer this question. However there are three areas where 
important clues may emerge: the light-neutrino mass matrix M,^; the baryon asymmetry 
of the universe; and (assuming supersymmetry) lepton flavor violation. Taken together 
with other theoretical ideas, we shall show that it may be possible to shed light on the 
right-handed neutrino sector of the theory. 



5.1.1 The Three Right-Handed Neutrino Paradigm 

It is most common to assume that there are exactly three right-handed neutrinos. Such 
an assumption is motivated by unified theories such as S'O(IO) which predicts that the 
number of right-handed neutrinos is equal to the number of quark and lepton families, 
since a single right-handed neutrino makes up each 16-plet of the theory. In fact this 
prediction also follows more generally from any theory which contains a gauged right- 
handed group SU{2)r, such as left-right symmetric theories, Pati-Salam and so on. 

Assuming three right-handed neutrinos one can ask whether their mass spectrum is 
hierarchical, or contains an approximate two or three-fold degeneracy. From the point of 
view of the type I seesaw mechanism in Eq. fl4()j) it is clear that the answer to this question 
depends on the nature of the Dirac neutrino mass matrix . For example, suppose that 
the right-handed neutrinos had a three- fold degeneracy M^i = diag(M, M, M), then the 
Eq. (jinj) would predict 

M. = "-jf^. (41) 

Then if the Dirac neutrino mass matrix were hierarchical and approximately proportional 
to the up-type quarks, for example, then Eq. PT|) would imply 

nil ■ ^2 ■ TTT'S ~ ^1 '■ TTT'I '■ (^2) 
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Here we shall assume M^, = Ml. We shall later comment also on the type II seesaw mechanism. 
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which is much too strong a mass hierarchy compared to the rather mild experimentally 
measured ratio 0.1 < 1712/ < 1. 

The remaining two possibilities are that the three right-handed neutrinos are either 
hierarchical or contain an approximate two-fold degeneracy. In either case it is convenient 
to work in a basis where their mass matrix is diagonal: 

( X' \ 
= X . (43) 

V r y 

The neutrino Dirac mass matrix M,f in this basis can be written as 



Ml^ =\ h' h e \ (44) 



a' 


a 


d 


b' 


b 


e 


d 


c 


f 



where in this convention the first column of Eq. ()44j) couples to the first right-handed 
neutrino, the second column of Eq. couples to the second right-handed neutrino and 
so on. Note that in the hierarchical case in Eq. ()43|) we do not specify which of the three 
right-handed neutrinos X',X,Y is the lightest one, which is the intermediate one and 
which is the heaviest one, since the columns of and the eigenvalues X', X, Y of Mji 
may simultaneously be re-ordered without changing M^. 

Having displayed the unknown Yukawa couplings associated with the Dirac neutrino 
mass matrix in Eq. (jl^ it is clear that without further input it is not possible to say 
anything about the right-handed neutrino masses or couplings from the experimentally 
determined light Majorana mass matrix M^,. On the other hand, rather natural theoretical 
assumptions can lead to a great deal of information about the unknown masses and 
couplings of the right-handed neutrinos, as we now discuss. 

Regarding the implementation of the type I seesaw mechanism there seem to be two 
possible options: either all the right-handed neutrinos contribute equally (democratically) 
to each element of M^, or some right-handed neutrinos contribute more strongly than 
others. In the second case, called right-handed neutrino dominance [241], a rather natural 
implementation of the seesaw mechanism is possible. For example if the right-handed 
neutrino of mass Y contributes dominantly to the physical neutrino mass m^, and the 
right-handed neutrino of mass X contributes dominantly to the physical neutrino mass 
777.2, while the right-handed neutrino of mass X' contributes dominantly to the physical 
neutrino mass mi, then a sequential dominance of these three contributions leads to a 
neutrino mass hierarchy mi <^ m2 ^ m^. With such sequential dominance the mixing 
angles are then given as simple ratios of Dirac neutrino mass matrix elements: tan 623 ~ 
e//, tan^i2 ~ ^/2a/{b — c), which can be naturally large independently of the neutrino 
mass hierarchy. 

The physical neutrino masses are given by ms ^ (e^ -|- p)/Y, 1712 ~ Aa^/X, mi ^ 
(a', b', d)'^/X' and the mass ordering of the right-handed neutrino masses X', X, Y is not 
determined unless further information is specified about the Dirac neutrino masses. In 
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general there are six possible mass orderings of the three right-handed neutrinos: 



Y <X <X' 

Y <X' <X 
X <Y <X' 
X' <Y <X 
X' <X <Y 
X <X' <Y 



(45) 
(46) 
(47) 
(48) 
(49) 
(50) 



The dominant right-handed neutrino of mass Y (the one mainly responsible for the mass 
ms) may thus be the lightest one as in Eqs. PH| the intermediate one as in Eqs. 
()47[ HHj) . or the heaviest one as in Eqs. ()491 ISUj) . The neutrino of mass X' is essentially 
irrelevant from the point of view of the light Majorana mass matrix M^, since the lightest 
physical neutrino of mass nii is approximately zero in the hierarchical case. Thus X' 
cannot be determined from any low energy experiments such as neutrino oscillations or 
neutrinoless double beta decay. If X' happens to be the heaviest right-handed neutrino, 
as in Eqs. (^31 ITTj) then if its mass is above the GUT scale then it completely decouples 
from observable physics. In this case the three right-handed neutrino model becomes 
effectively a two right-handed neutrino model [241]. However even in this case there is 
a remaining ambiguity about whether the dominant right-handed neutrino is the lightest 
one or the next-to- lightest one as in Eqs. (jlHt ITTjl . 

5.1.2 Grand Unification 

It is clear that further theoretical input is required in order to elucidate the nature of 
the masses and couplings of the right-handed neutrinos. In some GUT model, one can 
expect generically that the Dirac neutrino masses are related to the other quark and 
charged lepton masses and this additional information about M^P can then be input into 
the type I seesaw formula Eq. pUj) to help to yield information about the right-handed 
neutrino mass matrix Mr- For example assuming that M", the up-type quark 

mass matrix, and inputting the approximately determined light Majorana mass matrix, 
the seesaw formula can then be rearranged to yield right-handed neutrino masses with a 
very hierarchical mass spectrum [242]: 



which can be compared to the naive expectation for the physical neutrino masses in 
Eq. (021) • Numerically Eq. ()5H1 yields the order of magnitude estimates Mi ~ 10^ GeV, 
M2 ~ 10^° GeV, Ml ~ 10^^ GeV, with an uncertainty of a least one or two orders 
of magnitude in each case. In addition there may be special cases which completely 
invalidate these estimates. 

In specific GUT models the above expectations can also be very badly violated. For 
example in the S'O(IO) model with S't/(3) family symmetry [243], although the neutrino 
Dirac mass matrix is strikingly similar to the up-type quark mass matrix, a very different 
pattern of right-handed neutrino masses emerges: 



Ml : M2 : M3 ^ ml : 




(51) 



Ml : M2 : M3 ^ e^f : e^f : 1 



(52) 
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where e ~ 0.05, e ~ 0.15. In this model the dominant right-handed neutrino is the hghtest 
one Y = Ml, with X = M2, while the heaviest right-handed neutrino is decoupled 
X' = M3 as in Eq. ()45|) . This model therefore acts effectively as a two right-handed 
neutrino model, with the two right-handed neutrinos being very similar in mass, with 
interesting implications for leptogenesis, to which we now turn. 

5.1.3 Leptogenesis 

Leptogenesis and lepton flavor violation are important indicators which can help to resolve 
the ambiguity of right-handed neutrino masses in Eqs. ()45ll50p . In the simplest case of two 
right-handed neutrino models, leptogenesis has been well studied with some interesting 
results [196-198]. In general, successful thermal leptogenesis for such models requires the 
mass of the lightest right-handed neutrino model to be quite high, and generally to exceed 
the gravitino constraints, required if supersymmetry is assumed. Such a strong bound is 
also at odds with the strong right-handed neutrino mass hierarchy expected from GUTs 
as in Eqs. ^ET\ 13^ . In unified theories with type II see-saw, this potential problem can 
be resolved (see e.g. [244]). 

In three right-handed neutrino models with sequential dominance, if the dominant 
right-handed neutrino is the lightest one, then the washout parameter mi ~ (^(ms) is 
rather too large compared to the optimal value of around 10~^ eV for thermal leptogenesis. 
However, if the dominant right-handed neutrino is either the intermediate or the heaviest 
one then one finds mi ~ 0{m2) or arbitrary mi, which can be closer to the desired 
value [245]. 

5.1.4 Sneutrino Inflation 

It has been suggested that a right-handed sneutrino, the superpartner to a right-handed 
neutrino, could be a candidate for the inflaton in theories of cosmological inflation. This 
has interesting consequences for the masses and couplings of the right-handed neutrinos. 
For example in the case of chaotic sneutrino inflation, the mass of the right-handed sneu- 
trino inflaton must be about 10^^ GeV [246], while in sneutrino hybrid inflation its mass 
could be considerably lighter [247] . In both scenarios, the decaying sneutrino inflaton can 
be responsible for non-thermal leptogenesis, and can give a reheat temperature compatible 
with gravitino constraints providing its Yukawa couplings are sufficiently small. This typi- 
cally implies that the associated right-handed neutrino must be effectively decoupled from 
the see-saw mechanism, so that it corresponds to the decoupled right-handed neutrino of 
mass X' in sequential dominance discussed above. 

5.1.5 Type II Seesaw Models 

Once the more general type II seesaw framework is permitted [73], then it apparently be- 
comes more problematic to determine the properties of the right-handed neutrinos which 
contribute to via the type I part of the seesaw mechanism. On the other hand, the 
type II seesaw mechanism provides the most direct way of raising the neutrino mass scale 
to a level that will be observable in neutrinoless double beta decay. Furthermore, the 
difficulties of providing consistency of leptogenesis scenarios with the gravitino bound in 
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supcrsymmetric theories, or simply with the strong right-handed neutrino mass hierarchy 
expected from GUT models, motivates a more general type II seesaw framework which 
can in principle help resolve some of these difficulties. It has recently been shown [104,248] 
how to construct natural models for partially degenerate neutrinos by using an SO{3) fam- 
ily symmetry to add a type II contribution to the light neutrino Majorana mass matrix 
proportional to the unit matrix, with large neutrino mixing originating from sequential 
dominance. Compared to the pure type I limit, the masses of the right-handed neutrinos 
become larger if the mass scale of the light neutrinos is increased via the type II contribu- 
tion. This can also help to resolve the potential conflict between the typical predictions 
for Ml as in Eqs. (47,48) and thermal leptogenesis [244]. In addition, increasing the neu- 
trino mass scale has interesting phenomenological consequences, such as a decreasing CP 
violating phase 6 and a decreasing mixing angle 6*13, testable in future experiments. 

5.1.6 Right-Handed Neutrinos in Extended Technicolor 

In the mechanism that has been constructed for producing light neutrinos in extended 
Technicolor [92,93,441,442], there are two right-handed neutrinos. Interestingly, this 
mechanism involves a seesaw, but one in which the relevant Dirac neutrino mass matrix 
elements are greatly suppressed down to the level of a few keV, and the Majorana masses 
are also suppressed, of order 100 MeV to 1 GeV. The origin of this suppression stems 
from the fact that the left- and right-handed chiral components of neutrinos transform 
differently under the ETC gauge group. Although the mechanism does involve a seesaw, it 
does not involve any GUT-scale masses. This is clear, since extended Technicolor models 
do not contain any such mass scales. It serves as an existence proof of how a seesaw 
mechanism can work with much lower Dirac and Majorana mass scales than the usual 
GUT-scale seesaw. 

5.2 Seesaw Neutrino mass and Grand unification 

One of the major ideas for physics beyond the Standard Model is supersymmetric grand 
unification (SUSY GUT) [249]. It is stimulated by a number of observations that are 
in accord with the general expectations from SUSY GUTs : (i) A solution to the gauge 
hierarchy problem i.e why t^^k <^ -^Ph (h) unification of electroweak, i.e. SU{2)l x 
U{1)y and strong SU{3)c gauge couplings assuming supersymmetry breaking masses are 
in the TeV range, as would be required by the solution to the gauge hierarchy; (iii) a 
natural way to understand the origin of electroweak symmetry breaking. Supersymmetric 
grand unified theories generically predict proton decay via dimension-five operators, which 
typically give large branching ratios for modes like p ^^iK^ ■ The current lower limits on 
proton decay modes place significant constraints on these theories and probably rule out 
a number of simpler SUSY GUT models [250]. Nevertheless, the idea of grand unification 
is so attractive that we will proceed on the basis that appropriate modifications allow 
supersymmetric GUTs to evade current nucleon decay limits. 

Gauge coupling unification leads to a unification scale of about 10^^ GeV and simple 
seesaw intuition leads to a seesaw scale of 10^^ GeV to fit atmospheric neutrino data. This 
suggests that seesaw scale could be the GUT scale itself; thus the smallness of neutrino 
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mass could go quite well with the idea of supersymmetric grand unification. However, in 
contrast with the items (i) through (iii) listed above, the abundance of information for 
neutrinos makes it a highly nontrivial exercise to see whether the neutrino mixings indeed 
fit well into SUSY GUTs. In turn, the freedom in constructing reahstic GUT models 
allows many different ways to explain current neutrino observations. Thus even though 
neutrino mass is a solid evidence for physics beyond the Standard Model, the true nature 
of this physics still remains obscure. The hope is that the next round of the experiments 
will help to narrow the field of candidate theories a great deal. 

To see how this is likely to come about, the first point is the choice of the grand 
unification group. Even though attempts to implement the seesaw mechanism using an 
extension of the SU(5) with the addition of a right-handed neutrinos have been made, a 
more natural GUT gauge group from the point of view of neutrino mass is SO (10) since its 
basic spinor representation contains the right-handed neutrino automatically along with 
the other fifteen fermions of the Standard Model (for each family). Thus in some sense 
one could argue that small neutrino masses have already chosen 5*0 (10) GUT as the most 
natural way to proceed beyond the Standard Model. 5*0(10) has therefore rightly been 
the focus of many attempts to understand neutrino mixings. It turns out that within the 
5*0(10) SUSY GUTs there are many ways to understand large mixings. We outline below 
only the major differences among the different ideas. The hope is that they differ in their 
predictions sufficiently so that they can be tested by planned experiments. 

One of the features that distinguishes 5*0(10) from SU(5) is the presence of local 
B — L symmetry as a subgroup and the 5*0(10) models divide into two classes depending 
on whether B — L symmetry is broken by a 16 Higgs field or an 126. In the first case the 
right-handed neutrino mass necessarily arises out a nonrenormalizable coupling whereas 
in the second case it arises from a renormalizable one. Secondly, the breaking of S — L by 
16 Higgs necessarily leads to low energy MSSM with R-parity breaking so that the model 
cannot have cold dark matter without additional assumptions, whereas 126 breaking of 
B — L preserves R-parity at low energies so that the low energy MSSM that derives from 
such an 5*0(10) has a natural dark matter candidate i.e. the lightest SUSY particle. 

As noted in the Introduction, the 5*0(10) model has in general the type II seesaw 
formula for neutrino masses which can reduce to type I for some range of parameters. For 
instance, in the 16 based models, the first term in type II seesaw formula is negligible 
and therefore the neutrino masses are dictated by type I seesaw formula. In contrast in 
126 Higgs models, the neutrino mass can be given either by the first term or the second 
term in the type II seesaw formula or both. 

5.2.1 A minimal 126-based 50(10) model 

As mentioned, in 5*0(10) models where a 126 Higgs breaks B — L symmetry, the right- 
handed neutrino masses can arise from renormalizable couplings. A minimal model of 
this type based on a single 10 and a single 126 field has a number of attractive features 
[251,252]. Since the 126 field also contributes to charged fermion masses through the 
MSSM doublets in it, this model unifies the fiavor structure in the quark and the neutrino 
sector thereby increasing the predictivity of the model in the neutrino sector. In fact in 
the absence of CP violation, the model has only 12 free parameters all of which are 
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Model 


^13 


126 based models 
Goh, Mohapatra, Ng 

V^iit;!!, iVlcLlicLllLiicljJjJcL 


0.18 
n 1 ^ 


16 based models 
Babu, Pati, Wilczek 

Albright, Barr 
Ross, Velasco-Sevilla 
Blazek, Raby, Tobe 


0.0005 
0.014 
0.07 
0.05 



Table 6: The table lists some typical predictions for ^13 in different 50(10) models and 
shows how the next generation of experiments can narrow the field of possible SO(IO) 
unification models. 

determined by the quark masses and mixings along with the charged lepton masses. As 
a result all mixings and masses are predicted by the model up to an overall scale. It has 
been shown that if one uses the type I seesaw formula, the model fails to reproduce the 
observed solar mixing angle and also the solar mass difference squared and is therefore 
ruled out [253]. It has however been shown recently that if one uses the type II seesaw 
mechanism with the first term dominating the neutrino mass matrix, the large mixings 
come out due to b-tau mass convergence in a very natural manner [254]. In particular, 
an interesting prediction of this model is that 6^13 ~ 0.18 making it quite accessible to the 
next generation of experiments such as long-baseline, off-axis and the reactor experiments. 

5.2.2 16-based models 

The main characteristic of the S'O(IO) models where a 16 Higgs breaks B — Lis that right- 
handed neutrino masses arise from nonrenormalizable couplings in the superpotential, the 
implicit assumption being that there is a high-scale theory (perhaps string theory or a 
renormalizable high scale theory with heavier fields) that below the heavy scale leads to 
this version of S'O(IO). This means that without additional symmetry restriction, there 
are more parameters than the physical inputs. Often in these models symmetries that 
tend to explain quark mixings restrict the number of couplings somewhat and one can 
make predictions in the neutrino sector. There exist several interesting examples of this 
kind of models [255]. Several of these models tend to give values for ^13 which are much 
below the range that can be probed by the next generation of planned experiments. We 
give a very small sample of the different predictions for in models with both 16 and 
126 in Table El 

5.2.3 Summary of what we can learn about SO{lQ) 

A review on different neutrino mass models based on 5*0(10) can be found in Ref. [256]. 
From these models, we learn that 

1. First a very generic prediction of all 5*0(10) models is that the neutrino mass hier- 
archy is normal. The basic reason for this is the quark lepton symmetry inherent in 
the model, which tends to make the neutrino Dirac mass to be of similar hierarchy 
as the quarks, which via seesaw mechanism implies normal hierarchy for neutrinos. 
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leading to ^m^^ > 0. Again this is a result that can be probed in the long-baseline 
oscillation or neutrinoless double beta decay experiments; 

2. The second point about the SO{10) models is that they make definite predictions 
about the mixing angle ^13 as given in Table El and often for the other mixing 
angles. The planned experiments will therefore considerably narrow the field of 
viable 5*0(10) models through their measurement or upper limit on these mixing 
parameters. 

5.2.4 Implications of Models with Spontaneous CP Violation 

Relations between leptogenesis and CP violation in low energy processes generally do not 
exist due to the presence of un-known mixing angles and phases in the heavy neutrino 
sector, as discussed in Sec. 14.2.41 In models with spontaneous CP violation, all Yukawa 
coupling constants are real. CP violation occurs due to the presence of the phases in the 
expectation values of the scalar fields, which break the gauge symmetry spontaneously. 
Recently it has been shown that [257] in the minimal left-right SU (2)^ x SU {2)r symmetric 
model [72] with spontaneous CP violation, there exist very pronounced relations between 
the CP violation in low energy processes, such as neutrino oscillation and neutrinoless 
double beta decay, and leptogenesis, which occurs at a very high energy scale. The 
minimal left-right symmetric model contains a bi-doublet and a pair of triplet Higgses. 
Using the gauge degrees of freedom, one can rotate away all but two phases present in 
the expectation values of the scalar fields. Thus there are only two intrinsic phases, the 
relative phase between the two vevs in the bi-doublet and that between the left- and 
right-handed triplets, to account for all CP violation in the quark sector and in the lepton 
sector. The relative phase between the two vevs in the bi-doublet is responsible for CP 
violation observed in the quark sector, while CP violation in the lepton sector dominantly 
comes from the relative phase between the vevs of the two triplet Higgses. The relative 
phase between the two vevs in the bi-doublet appears in the lepton sector only at the 
sub-leading order due to the large hierarchy in the bi-doublet vevs required by a realistic 
quark sector. As a result, the relation between CP violation in the quark sector and CP 
violation in the lepton sector is rather weak. Due to the left-right parity, the RH and 
LH neutrino Majorana mass terms are proportional to each other, which further reduces 
the unknown parameters in the model. In this model, both leptogenesis and the leptonic 
Jarlskog invariant are proportional to the sine of the relative phase between the vevs of the 
two triplet Higgses. Using the experimentally measured neutrino oscillation parameters 
as inputs, to obtain sufficient amount of lepton number asymmetry, the leptonic Jarlskog 
invariant has to be larger than ~ 10~^. As the Type-II seesaw mechanism is at work, the 
hierarchy in the heavy neutrino sector required to obtain the observed neutrino oscillation 
parameters is very small, leading to a heavier mass for the lightest RH neutrino, compared 
to the case utilizing the Type-I seesaw mechanism. As a result, the requirement that the 
decay of the lightest RH neutrino is out-of-equilibrium in thermal leptogenesis can be 
easily satisfied. 

Similar attempts have been made to induce spontaneous CP violation from a single 
source. In one such attempt SM is extended by a singlet scalar field which develops a 
complex VEV which breaks CP symmetry [258]. Another attempt assumes that there is 
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one complex VEV of the field which breaks the B — L symmetry in SO (10) [259]. Unlike 
in the minimal left-right symmetry model described above [257], there is no compelling 
reason why all other vevs have to be real in these models. 



5.3 Renormalization group evolution of neutrino parameters 

Neutrino masses and mixing parameters are subject to the renormalization group (RG) 
evolution or running, i.e. they depend on energy. As theoretical predictions for these 
quantities typically arise from models at high energy scales such as the GUT scale, this 
implies that in general RG corrections have to be included in the testing of model predic- 
tions. In the case of leptonic mixing angles and CP phases the changes can be large for 
partially or nearly degenerate neutrino masses. On the other hand, strongly hierarchical 
masses bring about very small RG corrections for the mixing angles, while the running of 
the mass squared differences is sizable even in this case. 



5.3.1 Running masses, mixings and CP phases below the seesaw scale 

At energies below the seesaw scale Mr, the masses of the light neutrinos can be described 
in a rather model-independent way by an effective dimension 5 operator if they are Ma- 
jorana particles. The RG equation of this operator in the SM and MSSM [80,260-263] 
leads to differential equations for the energy dependence of the mass eigenvalues, mixing 
angles and CP phases [264-268]. Up to 0(^13) corrections, the evolution of the mixing 
angles is given by [266] 

327r^ AuIq 

^13 = :^7r2 sin 26*12 sin 26*23-7 — ^ /1 , .-a ^ Hmi,(pi,S) , (53) 
327r"^ Am^ (1 + Q 



n _ 1 [ 2 |_ , ^ ,2 , „2 \m,e^^^+m,\'] 



C12 \iTT'2 e^^ + m^\ + s 



12 



1 + C 



where the dot indicates the differentiation d/dt = ^d/dfx (/x being the renormalization 
scale), Sij := sin6'jj, Cij := cos 9ij, ( := AttIq/ Am\, C = —3/2 in the SM and C = 1 in 
the MSSM, and /(m^, ipi, S) := [mi cos((y9i — 5) — (1 + C) ^2 cos(v92 — 6) — (m^ cos 5]. 
denotes the r Yukawa coupling, and one can safely neglect the contributions coming from 
the electron and muon. For the matrix K containing the Majorana phases, we use the 
convention K = diag (e~*'^^/^, e"""^^/^, 1) here. For a discussion of RG effects in the case of 
exactly degenerate neutrino masses, where the above expressions cannot be applied, see 
e.g. [269-277]. iFrom Eqs. (jSHjl one can easily understand the typical size of RG effects 
as well as some basic properties. First, in the SM and in the MSSM with small tan/?, 
the RG evolution of the mixing angles is negligible due to the smallness of the r Yukawa 
coupling. Next, the RG evolution of the angles is the stronger the more degenerate the 
mass spectrum is. For a strong normal mass hierarchy, it is negligible even in the MSSM 
with a large tan/?, but for an inverted hierarchy a significant running is possible even if 
the lightest neutrino is massless [278,279]. Furthermore, non-zero phases 6, ipi and ip2 
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can either damp or enhance the running. For instance, the running of 612 can be damped 
by non-zero Majorana phases [81,280,281]. Typically, 6*12 undergoes the strongest RG 
evolution because the solar mass squared difference is much smaller than the atmospheric 
one. Finally, in the MSSM, 612 runs from smaller values at high energies to larger values 
at low energies [282]. 

The RG equations for the CP phases [265-268] show that their changes are propor- 
tional to l/ArriQ. Therefore, whenever the mixing angles run sizably, the same happens 
for the CP phases. This is very important for the relation between the phases relevant for 
high-energy processes like leptogenesis and those appearing in neutrino oscillations and 
neutrinoless double beta decay. The evolution of the Dirac CP phase can be especially 
drastic for a small CHOOZ angle, since 6 contains a term proportional to 1/^13. It is 
also possible to generate a non-zero value of this phase radiatively if at least one of the 
Majorana phases is non-zero [265]. An exception is the CP conserving case where all 
phases are or tt and do not change with energy. 

Finally, the neutrino masses always change significantly with energy due to flavor-blind 
terms in the RG equations which contain large quantities like gauge couplings and the top 
Yukawa coupling. For strongly hierarchical masses and small tan/3, these terms dominate, 
so that the masses experience a common rescaling which is virtually independent of the 
mixing parameters [281]. 

Radiative corrections for Dirac neutrino masses have also been studied [283,284]. 
Roughly speaking, the RGBs for the Dirac case are obtained from Eqs. by averaging 
over the Majorana phases. 

5.3.2 Details of the running in seesaw models 

In order to obtain precise results, one has to go beyond the simple approximations listed 
above and solve the RG equations numerically. This involves solving a rather complex 
system of coupled differential equations, as all parameters of the theory have to be evolved 
from high to low energy. 

A further complication arises in seesaw models with heavy singlet neutrinos which are 
in general non-degenerate in mass. The running above their mass thresholds is typically 
at least as important as the evolution below both in the SM and in the MSSM unless the 
neutrino Yukawa couplings are tiny [82,270,271,285-289]. This part of the RG evolution 
depends on many parameters of the model. An analytic understanding has been obtained 
only recently [290,291]. If the singlet masses are non-degenerate, one can calculate the 
evolution of the neutrino mass parameters by considering a series of effective theories 
arising from integrating out the singlets successively at the respective thresholds [285,286]. 
In general, it is not a good approximation to integrate out all singlets at the same energy 
scale, since the threshold corrections can be very large. 

5.3.3 Implications for model building 

As discussed above, predictions of high-energy mass models can differ substantially from 
low-energy experimental results due to the running. Therefore, RG corrections have to 
be included in the analysis. The RG evolution also opens up new interesting possibilities 
for model building, like the radiative magnification of mixing angles. In particular, small 
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or vanishing solar mixing at high energy can be magnified to the observed large mixing at 
low energy (see e.g. [82,84,292,293]). Vice versa, the large but non-maximal solar mixing 
9i2 can also be reached starting from bimaximal mixing at the GUT scale [287, 294, 295] 
(for examples, see Fig. ll5|) . It is, however, important to stress that large mixing is no fixed 
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Figure 15: Examples for the RG evolution of the lepton mixing angles from the GUT 
scale to the SUSY-breaking scale (taken to be ~ 1 TeV) in the MSSM extended by 3 
heavy singlets (right-handed neutrinos) [82,287]. The masses of the lightest neutrinos for 
these examples are around 0.05 eV. The figures illustrate how the large but non-maximal 
value of the solar mixing angle 612 is reached by RG running if one starts with bimaximal 
lepton mixing or with vanishing solar mixing at the GUT scale. The kinks in the plots 
correspond to the mass thresholds at the seesaw scales, where the heavy singlets are 
successively integrated out. The gray-shaded regions mark the various effective theories 
between the seesaw scales. 



point under the RGE in the usual see-saw framework. It has been observed that in SUSY 
models large mixing can be a fixed point for different (i.e. non-seesaw) types of neutrino 
mass operators [296,297]. In addition, the small neutrino mass squared differences can be 
produced from exactly degenerate neutrino masses at high energy (see e.g. [275-277,298- 
300]), if the neutrino masses are nearly degenerate. For further specific models where the 
RG evolution is relevant for neutrino masses and mixings see, for example, [248,267,301- 
304]. 



5.3.4 High-scale mixing unification and large mixings 

Another question one can ask is whether starting with small mixing angles at the seesaw 
scale (as would be naively expected in models with quark-lepton unification) can one 
get large mixings at the weak scale due to RG extrapolation. In a specific model where 
neutrino masses are quasi-degenerate, starting with neutrino mixings that are equal to 
quark mixings at the GUT scale i.e. 6*12 — Vus, O23 — Vcb and 6'i3 ~ Vuh, one can indeed get 
mixing angles at the weak scale which are consistent with present observations as shown 
in Fig. El below. We have chosen tan /5 = 55 in this calculation. An interesting point 
is that this mechanism works only if the common mass of the neutrinos is bigger than 
0.1 eV, a prediction which can easily be tested in the proposed neutrinoless double beta 
decay experiments [81,83,305-307]. 
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Figure 16: Evolution of small quark-like mixings at the seesaw scale to bi-large neutrino 
mixings at low energies for the seesaw scale Mr = 10^^ GeV with tan/? = 55, Msusy = 1 
TeV. The solid, long-dashed and short-dashed lines represent sin 6*23, sin 6^13, and sin 6*12, 
respectively. The evolution of the sines of quark mixing angles, sin 9gij{i,j = 1,2,3), is 
presented by almost horizontal lines. 

5.3.5 Deviations of ^13 from and of 623 from 7r/4 due to RG effects 

At present, observations are compatible with ^13 = and ^23 = vr/4. New experiments are 
being planned to lower the limits on deviations from these values. Additional motivation 
for this kind of measurements is provided by the RG: although it is possible that a 
symmetry produces an exactly vanishing 61^ and exactly maximal atmospheric mixing, 
this symmetry would typically operate at a high scale, and therefore its predictions would 
be subject to the RG evolution. Hence, without fine-tuning, one expects non-zero values 
of ^13 and 023 - 7r/4 at low energy [266,277,288,308]. For example, in the MSSM one 
finds a shift Asin^ 2^13 > 0.01 for a considerable parameter range, i.e. one would expect 
to measure a finite value of 613 already in the next generation of experiments [266]. On 
the other hand, there are special configurations of the parameters, especially the phases, 
where RG effects are suppressed. Furthermore, there may be symmetries which stabilize 
some mixing angle completely against radiative corrections. For instance, for an inverted 
hierarchy with = 0, 6*13 = is stable under the RG [309] (see also Eqs. (j^Sl))- Hence, 
if future precision measurements do not find ^13 and 6'23 — 7r/4 of the size of the generic 
RG change, one can restrict parameters, or even obtain evidence for a new symmetry. 

5.3.6 Implications for leptogenesis 

As has been discussed in Sec. 14. 2| the requirement of successful baryogenesis via leptogen- 
esis places an upper bound on the mass of the light neutrinos. It is important to note that 
in order to relate constraints on the neutrino mass spectrum coming from physics at Mi 
to observation, one has to take into account radiative corrections. It turns out that there 
are two effects operating in opposite directions [183, 185, 266, 310] which partially cancel 
each other: since the mass scale is increasing, the washout driven by Yukawa couplings 
is stronger. On the other hand, larger Am^s allow for a larger decay asymmetry. Taking 
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into account all these effects, one finds that the upper bound on the neutrino mass scale 
becomes more restrictive [183,185,266]. The RG evolution, together with thermal effects 
or spectator processes [311], gives rise to the most important corrections to the mass 
bound [185]. 

6 Non- Standard Neutrino Interactions 

6.1 Neutrino magnetic moments 

Once neutrinos are massive, they can have transition magnetic and electric magnetic 
dipole moments and, in the Dirac case, also diagonal magnetic and electric dipole mo- 
ments [48,49], [40], [51]- [59]. The lepton-number conserving (diagonal and transi- 
tion) magnetic and electric dipole moment operators are given by (l/2)z/jO'^^z/jF^j, and 
{—i/2)uia'^^j5i'jF^i,. Analogous AL = 2 expressions hold (with i ^ j because of their 
antisymmetry under i <-> j) for Majorana neutrinos. Therefore, a magnetic or electric 
dipole moment always connects one species of neutrino with another. These moments are 
defined for mass eigenstates. In the case of a Dirac mass eigenstate, one sees that the op- 
erator connects a left-handed electroweak-doublet neutrino to a right-handed electroweak- 
singlet (sterile) neutrino. In the case of Majorana mass eigenstates, the (transition) dipole 
moment operators connect two neutrino fields of the same chirality. The two have funda- 
mentally different physical implications. We have reviewed some basic properties of these 
magnetic and electric dipole moments above. 

Neutrino magnetic moments can be directly measured in terrestrial experiments using 
the neutrino beam from the Sun as in Super-K [312] or with neutrinos from close by 
nuclear reactors as in the MUNU [313] and in the Texono [314] experiments. These 
experiments have put upper bounds of the order of 10~^°/xb on the effective neutrino 
magnetic moment (defined below). It is interesting that in models involving right-handed 
charged currents the diagonal and transition neutrino magnetic moments [40,48-51,53] 
are not suppressed by the neutrino mass (as in models with just the Standard Model 
interactions [52]) and could thus be somewhat larger. However, in general, the same 
interactions that can enhance neutrino magnetic moments can give corrections enhancing 
neutrino masses, so in a particular model one must be careful to avoid excessive loop 
contributions to the latter. The Borexino prototype detector has recently been utilized 
to put a bound of \Heff\MSW < 5.5 x IO'^^/ib at 90% C.L. using the elastic scattering 
of electrons by the solar pp and '^Be neutrinos [315]. At these sub-MeV energies, the 
solar neutrino beam contains roughly equal proportion of and 1/2 (-Pi — P2 — 0.5). As 
a result, for the same u^ff the bounds on ^11,^13 would be much better in experiments 
utilizing the ppj Be neutrinos compared to the bounds that can be obtained by using the 
^B neutrinos (which are predominantly 1/2)- 

Neutrinos with non-zero magnetic moments contribute to the elastic scattering of 
electrons in water Cerenkov detectors [316-318]. The effective neutrino magnetic moment 
Heff (we neglect the contribution of the electric dipole moment here) responsible for the 
scattering event Ui + e~ Vj -\- e~ is proportional to the incoherent sum of outgoing 
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neutrino states as follows: 



(54) 



where Ai{L) is the probability amplitude of a neutrino produced as a flavor eigenstate 
(lets say z/g or Uf.) to be in the i'th mass eigenstate on propagating over the source-detector 
distance L. For vacuum oscillations Ai{L) = Uei exp{—iEiL) and the effective magnetic 
moment depends upon the Am^ and the mixing angles as follows 



In the expression for \Heff\vo above we have assumed that is negligibly small and the 
atmospheric mixing angle is maximal (^23 = 1/2). We have also dropped CP violating 
phases. In MUNU [313] and Texono [314] , Uf. from nuclear reactors were detected by 
the elastic scattering with e~. The source-detector distance is small (L = 18m in MUNU 
and L = 28m in TEXONO) compared to the ui — 1/2 oscillation length so that the Cosine 
term in is unity. The magnetic moment matrix fiij is symmetric for Dirac neutrinos 
and anti-symmetric if the neutrinos are Majorana. It is clear from (j55|) that there exists 
a possibility that there may be a cancellation between the last (interference) term which 
could be negative and the first two terms. So experimental upper bounds on \fieff\vo 



which is 0.9 x IQ-^Vi? (MUNU) [313] and 1.3 x IQ-^Vs (Texono) (both at 90% C.L.) 



do not constrain the elements of the /ijj matrix without making added assumptions that 
there is no cancellation between the different terms in ()55|) . 

This problem does not arise for solar neutrinos as the interference term averages to 
zero since 2£'y/Am^2 << Learth-sun- For the solar neutrinos the expression for |/ie//P 
reduces to a sum of two positive definite quantities 



where Pi = |y4ei(L)p and P2 = \Ae2{L)\'^ = 1 — Pi are the probabilities of the solar 
neutrinos to be in the mass eigenstate z/i and 1^2 respectively at the earth. The recent upper 
bound \^eff\Msw < 1-1 x 10~^°yUB at 90% C.L. established by Super-Kamiokande [312] 
can be translated into bounds on individual elements of Hij without extra assumptions. 
The neutrinos which are detected by electron scattering at Super-K are predominantly 
z/2 state (P2 = 0.94 and Pi = 0.06). The Super-K bound on \fieff\Msw therefore implies 
Ifiul < 1.1 X 10~^°/iB ; |/i22|, |/i23| < 1-13 X IO'^^hb and \nn\, l/iial < 4.49 x lO'^^fis- 

It is also possible to put bounds on /ijj from SNO-NC data using the fact that neutrinos 
with non-zero magnetic moments can dissociate deuterium [319] in addition to the weak 
neutral currents. The bounds established from SNO-NC data do not depend upon the 
oscillation parameters unlike in the case of Super-K. However the bounds are poorer due 
to the large uncertainty in our theoretical knowledge of the theoretical ®P flux from the 
sun [320]. 

We will see in a subsequent section (the one on extra dimensions) that the effective 
magnetic moment of the neutrinos can get substantially enhanced in a certain class of 



f^efflvo 



+ 2Ci2 Si2 (yUll/i21 + yUl2/i22 + fJ'13fJ'23) COS 




1 2 / 2 2 2 \ / 2 2 2 \ 

MflMSW = Pi (/^ll + /^12 + /^is) + (/i21 + /^22 + /^23) 



(56) 
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extra dimensions models. Searching for /ij^ can therefore be used to put hmits on theories 
with extra dimensions. 

6.2 Flavor changing and conserving nonstandard neutral cur- 
rent interactions 

The latest results of neutrino oscillation experiments indicate that the conversion mech- 
anism between different neutrino flavors is driven by a non-vanishing mass difference 
between mass eigenstates together with large mixing angles between families. These anal- 
yses are done supposing that no non-standard neutrino interactions (NSNI) are present. 
In the presence of electroweak-doublet and electrowcak-singlet neutrinos, the neutral weak 
current is, in general, nondiagonal in mass eigenstatc neutrino fields [40,41]. This is the 
same type of nondiagonality in the neutral weak current that was present in the original 
Weinberg electroweak model before the advent of the Glashow-Iliopoulos-Maiani (GIM) 
mechanism. It will be recalled that in this original Weinberg model the d and s quarks 
were assigned to a left-handed 811(2)^ doublet {u,dcos9c + s sin 9c) l and a left-handed 
SU(2)i singlet —dsinOc + ssin^^c*. The necessary condition for the diagonality of the 
neutral weak current is that all of the fermions of a given charge and chirality must 
transform according to the same weak T and T3. Alternatively, NSNI in the form of 
nondiagonal couplings between different neutrino flavor eigenstates [321] and/or neutrino 
flavor-diagonal but fiavor nonsymmetric couplings, can exist [322]. Including NSNI can 
modify the characteristics of neutrino conversion, and in general large values of NSNI pa- 
rameters worsen the quality of the fit to data. We summarize here the present limits that 
can be obtained to NSNI parameters, using the result of neutrino oscillation experiments. 

6.2.1 Atmospheric neutrinos 

As repeatedly mentioned, atmospheric neutrino data are well described by the oscillation 
driven by one mass scale, Am|2, and with maximal mixing between second and third fam- 
ilies. One important prediction for these numbers is that the high-energy neutrino events 
that generate the through-going muon data are well described together with the low- 
energy neutrino events, due to the energy dependence of the Hamiltonian that describes 
the neutrino evolution. 

Assuming a non-vanishing NSNI acting together with mass and mixing, the solution 
to the atmospheric neutrino discrepancy can be spoiled if the NSNI parameters have too 
large values. This happens because the NSNI entries in the Hamiltonian that describes 
the neutrino evolution are energy independent [322]. Since a simultaneous explanation of 
low-energy and high-energy neutrino events requires a strong energy dependence in the 
u^, Vt conversion probability, inclusion of energy independent terms in the Hamiltonian 
tends to decrease the quality of the theoretical predictions fit to atmospheric neutrino 
data. 

The NSNI can be parametrized as a relative strength of such interactions to the e^- = 
where / stands for the fermion involved in the new interaction: 
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The solar neutrino data is able to stablish new limits in the flavor- diagonal NSNI ()57p . 
With the definitions written above, these limits stands for: 

In [323] an analysis of atmospheric neutrinos and NSNI is performed, and found the 
following limits, at 3a: 

M < 0.03 
|6;^-e;,| < 0.05 (57) 

These bounds refer to NSNI with d-quarks, and were obtained assuming a two-fiavor 
(z/^ and Ur) system. A three family analysis significantly relaxes some of the bounds 
above [324], such that order one values for some of the |e|s are not currently ruled out. 
For many more details, see [324]. 

For u-quarks the bounds are expected to be of the same order, and for NSNI with 
electrons we expect bounds looser by a factor of ~ 3. 



6.2.2 KamLAND and solar neutrinos 

The excellent agreement between the LMA parameter region that provides a solution 
to the solar neutrino problem and the parameter region compatible with KamLAND 
dataprovides us with an opportunity to use these data sets to establish a limit on non- 
standard neutrino interactions (NSNI). 

The effect of NSNI is negligible in KamLAND due to the short distance traveled inside 
the earth, but due to the high density and long travel distances in the sun, the presence of 
a NSNI could displace the best fit point of solar neutrino analyses, and spoil the agreement 
between solar and KamLAND allowed regions. 

The oscillation of solar neutrinos is driven by only one mass scale, Am^i. The higher 
mass scale relevant for atmospheric neutrino oscillations, decouples, and the mixing 

between the first and third family is very small, and will be set to zero in what follows. 

In this approach, after rotating out the third family from the evolution equation, we 
can write the 2x2 Hamiltonian that describes the neutrino evolution as 



H 



MSW 



+ 



Ne{r) - ^cos2e ^ sin 20 
^ sin 20 ^ cos 20 



V2GFefNf{r) 
y/2GFefNf{r) ^/2GFt'fN f{r) 



(5^ 
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where Nf{r) is an effective density felt by the neutrino, given by Nf — Np + 2A^„ for 
d-quarks, Nf = 2Np + for u-quarks and Nf = Np for electrons. 
The NSNI parameters can be written as: 



^/ ^ 2^23 - sin 26123 - ^ee 

ef = ee;, cos 6*23 - Cer sin 6I23 ■ (59) 



The atmospheric neutrino analyses have put strong bounds on e^r and e^^ — err, and 
we expect a near to maximal mixing between second and third families (cos 2^23 — 0). By 
further assuming that ef — (which would be the case if eg^u and were negligible) we 
can write e'f as 

^f ^i^l^ — = Ctt — (60) 

With the present data set and the assumptions listed above, we are able to establish 
the following limits to the NSNI parameters, at la {2a): 

- 0.20 < e' < 0.12 (e' < 0.30) d-quarks 
-0.18 < e' < 0.10 (e' < 0.30) u-quarks 

-0.55 < e' < 0.25 (e' < 0.86) electrons . (61) 

The limits obtained at 2a reflect the weak bounds on Am^ obtained by KamLAND. At 
present there are a number of possible "islands" in the parameter region compatible with 
KamLAND data, so the displacement in Am^ could make the best fit point of a combined 
analysis jump between consecutive islands. The increase of statistics at KamLAND will 
determine in which of these islands the true values of neutrino parameters lie, avoiding 
this kind of jump and improving the limits on the NSNI parameters. 

Simulating 1 kton-year of data at KamLAND, the allowed range of Am^ would be 
reduced, and in the case of an agreement between the neutrino parameters coming from 
KamLAND and solar neutrino analysis (that would indicate a vanishing NSNI) new limits 
could be obtained. Although the la regions do not change significantly, at 2a we have: 

- 0.42 < e' < 0.24 (d - quarks) 
-0.40 < e' < 0.18 (u - quarks) 

e' < 0.40 (electrons) . (62) 

Further increase of KamLAND statistics will not improve these bounds, which are now 
determined by the uncertainty in Am^ in the solar neutrino analysis. 

Details of the analysis presented here can be found in [325] . Similar analyses were also 
done in [326]. It should be pointed out that more severe constraints on (or a positive hint 
of) NSNI can be obtained in next-generation solar neutrino experiments, especially those 
sensitive to neutrino energies below 6 MeV [326]. 

6.2.3 Bounds from non-oscillating phenomena 

Apart from phenomena that involve neutrino oscillations, bounds on NSNI can also come 
from the effects of such non-standard interactions on the charged leptons [327,328]. We 
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should only be careful in translating such bounds to the neutrino sector, since usually 
this translation can only be possible with a few assumptions on the model that generates 
the non-standard interactions. Some bounds can also be found using neutrino scattering 
experiments [329,330]. We quote here some of the numbers obtained by these analyses. 
Details of the calculations can be found in the references. The following tables should be 
read as limits on ej^-, where i and j stand for e, /i and r, and are the lines and columns 
of the tables. 

The values quoted here depend on details of the models that generate the NSNI, 
such as SU{2)l breaking effects, absence of fine-tuning cancellations and the scale of new 
physics. Also, since neutrino oscillations are sensitive only to the vector coupling constant 
of the NSNI, correlations between the limits m and should be taken into account 
in order to compare the numbers presented here with the ones coming from neutrino 
oscillation experiments. 



Electrons: 
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/i 


r 
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-0.07 <eL< 0.1, 
-1< efl < 0.5 [330] 


e < 10-^ [328] 
eL,R < 5 X 10"^ [330] 


e < 4.2 X 10-^ [328] 
le^l <0.4, le^l <0.7 [330] 






|eL,i?| < 0.03 [330] 


e < 3.1 X 10"^ [327] 
eL^R < 0.1 [330] 


T 






[-0.05; 0.05]^,,^ [330] 



d-quarks: 
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/i 
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IclI < 0.3, 
-0.6 <tR< 0.5 [330] 


e < 10"^ [328] 
eL^R < 7.7 X 10-^ [330] 


e < 10-^ [328] 
eL^R < 0.5 [330] 






e < 0.1 [328] 
le^l < 0.003, -0.008 < eR < 0.015 [330] 


e < 1.2 X 10-^ [327] 
eL,R < 0.05 [330] 


T 






jeil < 1.1, 
|ei?| < 6 [330] 



u- quarks: 





e 
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-1 < eL < 0.3, 
-0.4 < e < 0.7 [330] 


e < 10"^ [328] 
eL,R < 7.7 X 10-4 [330] 


e < 10-^ [328] 
eL,R < 0.5 [330] 






e < 0.1 [328] 
le^l < 0.003, -0.008 < < 0.003 [330] 


e < 1.2 X 10-^ [327] 
eLR< 0.05 [330] 


T 






\eL\ < 1.4, \eR\ < 3 [330J 



Concluding, using the neutrino oscillation data we are able to find limits to NSNI 
parameters, without assuming any detail about the nature of new physics behind these 
interactions. More work is needed to improve the situation. 
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7 Beyond the three neutrino picture 



7.1 The search for other hght neutrinos 

A neutrino that does not participate in Standard Model interactions (i.e., is sterile) might 
seem of little interest, but this concept includes reasonable theoretical constructs such 
as right-handed neutrinos themselves. We note in passing that in one-family Technicolor 
theories there are also technineutrinos that would couple with the usual strength to the W 
and Z but, because of their Technicolor interactions, are confined and gain large dynamical 
masses of order several hundred GeV; they are therefore not relevant for usual low-energy 
neutrino oscillation experiments. The hypothesis of 'sterility' concerns the weak forces; 
gravity is expected to be felt anyway, and we cannot exclude that the 'sterile' neutrino 
participates in new forces, perhaps, mostly coupled to quarks; or carried by new heavy 
mediators; or that sterile neutrinos have preferential couplings with new particles — say, 
with Majorons. Even putting aside these possibilities, we can probe sterile neutrinos 
by the search for observable effects due to their mixing with the ordinary neutrinos. In 
this section, we will further restrict our attention on 'light' sterile neutrinos (say, below 
10 eV) and discuss the impact on oscillations. We make extensive reference to Ref. [331], 
an updated overview on the phenomenology of one extra sterile neutrinos. 

7.1.1 Issues of theoretical justification 

Many extensions of the Standard Model incorporate particles behaving as sterile neutrinos. 
The main question is: Why are these light and do they have the couplings needed to 
mix with ordinary neutrinos? A recent discussion is in Ref. [95]. Models with mirror 
matter [332, 333] contain mirror neutrinos [334, 335] and offer a straightforward answer: 
ordinary and mirror neutrinos are light for the same reason. It is easy to arrange a 
'communication' term between ordinary and mirror worlds, e.g., due to the operator 
~ z/0z/'0'/Mpi. This leads to long- wavelength oscillations into sterile neutrinos, (see Fig. 
IT7| from [339]). There are many other possibilities. For example, higher- dimensional 
operators in the superpotential in models involving a scalar field with an intermediate 
scale expectation value can naturally lead to small Dirac and Majorana masses of the 
same magnitude, and therefore to light ordinary and sterile neutrinos which can mix [95] . 
Already with mirror matter, the VEV (0') could be different from (0) = 174 GeV, and this 
has important carryings for the phenomenology [334, 335] . Alternatively, one could guess 
on dimensional ground the value TeV^/Mpi as the mass (or mixing) of sterile neutrinos, 
and relate the TeV- value, e.g., to supersymmetry breaking [336] or high GUT theories such 
as Eq [337]. In theories with dynamical electroweak symmetry breaking, the mechanism 
that has been found for light neutrinos predicts also that there are sterile neutrinos with 
masses of order lOO's of MeV to GeV [92,93,441,442]. 

Right-handed neutrino masses of about 1 or 100 TeV lead to exceedingly large masses 
that pose a question for extra dimensions if "large" scales should be absolutely avoided 
and if the Dirac Yukawa couplings are not small on their own. A popular way out is to 
postulate that the masses are Dirac in character. As a benefit, one explains why neu- 
trino masses are small when right-handed neutrinos propagate in the bulk [338]. But 
even if neutrinos turn out to be Dirac particles we would not have an evidence for these 
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theories, since (1) Dirac neutrino masses are possible in conventional 4- dimensional the- 
ories; (2) even in theories with extra dimensions, one can assume that the usual dim.-5 
term is the dominant source of neutrino masses. We just need to fix the scale of mass 
to the desired value of 10^^~^^ GeV, and neutrinos receive (presumably large) Majorana 
masses [339]. To summarize, small neutrino masses arc possible, and even more inter- 
estingly an infinite number of sterile neutrinos, but neutrino masses seem not to be a 
clean signal of low scale gravity theories. Oscillations have special character and there 
are interesting phenomenological constraints [340]. 

We believe that there should be increased attention paid not only toward phenomenol- 
ogy, but also toward theories of (light- or heavy-mass) sterile neutrinos. The fact that 
we do not understand usual fermion masses should be not taken as an excuse to avoid 
confrontation with theory in our view. With these considerations in mind, we will give 
some emphasis to mirror neutrinos, and at the same time point out the interest of under- 
standing the unknown 0{1) coefficients in these theories, that have an important impact 
on the phenomenological implications. 

7.1.2 Phenomenological manifestations 

In the following discussion, we will be concerned mostly with oscillations. However, the 
implications can be also elsewhere. To see that, it is sufficient to recall that when we add 
3 sterile neutrinos we can form Dirac masses, which means that there is no contribution 
to the neutrinoless double beta decay process. 

Terrestrial oscillation experiments Broadly speaking, there are two types of terres- 
trial experiments. The first one includes several disappearance experiments and LSND; 
the second one includes atmospheric neutrinos and long-baseline experiments. The first 
type is sensitive mostly to the mixing of and a sterile state, the other one also to or 
Ur. Both types of experiments probe only relatively large mixing angles, 9s ~ 0.1. Ster- 
ile neutrinos within the sensitivity regions are disfavored if standard cosmology (mostly 
BBN) apphes; further important tests will be done by CMB-I-LSS or BBN data. None 
of these experiments alone requires the existence of sterile neutrinos. A case for sterile 
neutrinos can be made interpreting in terms of oscillations LSND together with solar and 
atmospheric data [341]. The hypothesis that LSND signal is due to a relatively heavy and 
mostly sterile neutrino should be regarded as conservative [342], even though it leads to 
some problems with disappearance in terrestrial experiments, and interesting predictions 
for cosmology (BBN and CMB+LSS spectra). In view of this situation, the test of the 
LSND result is of essential importance. At the same time, we should not forget that 
sterile neutrinos could manifest themselves in other manners. 

Solar and KamLAND neutrinos The solar and KamLAND data can be explained 

well without sterile neutrinos. Even more, the 'LMA' solution received significant confir- 
mations: the sub-MeV energy regions have been probed by Gallium experiments and the 
super-MeV ones by SNO and Kamiokande, and LMA is in agreement with KamLAND. 
Thus we are led to consider minor admixtures of sterile neutrinos, presumably not more 
than 20 %. In many interesting cases sterile neutrinos are invisible at KamLAND but 
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Figure 17: The double degeneracy between mass eigenstates of ordinary and mirror world 
(z/j and ul) is lifted when the small mixing terms are included in the 6x6 mass matrix. 
The new mass eigenstates [u^ and z/~) are in good approximation maximal superpositions 
of Ui and z/-. 

affect the survival probability of solar neutrinos. Quite generally, to test the hypothesis of 
oscillations into sterile states it would be important to improve on (or measure precisely) 
the fluxes from Beryllium and pp-neutrinos. A few selected cases are shown in Fig. 



Neutrinoless double beta decay Massive Majorana sterile neutrinos, mixed with the 
active ones, would participate in mediating neutrinoless double beta decay. In the case of 
one light sterile neutrino, (^)e// is given by the sum of the contributions of four massive 
neutrinos [343]: 



i=l,4 

where rrii denotes the mass of the massive eigenstates and Uei indicates their mixing with 
the electron neutrino. 

Supernovae Supernovae are an ideal place to search for manifestations of sterile neu- 
trinos, either through long wavelength vacuum oscillations of a galactic supernova (whose 
distance is in the kpc scale-SN1987A was at 52 kpc) or through MSW effects (since nuclear 
densities are reached during the collapse); the main trouble is that we are still not able to 
understand supernova explosions theoretically (the data from SN1987A do not contradict 
the simplest current picture for neutrino emission, but they have puzzling features that 
suggest caution). Vacuum oscillations induced by mirror neutrinos (as in Fig. I17|l can led 
to a disappearance of half of the emitted flux [339] . MSW oscillations into sterile neu- 
trinos can produce even more dramatic suppressions; in certain regions of the parameter 
space this can reach 80 % [331]. The last type of oscillations are due to the fact that the 
SN core is deleptonized, and require that the electron neutrino mixes with the sterile one. 
The trouble to verify these predictions is the accuracy of the expectations on the emitted 
neutrino energy, that amounts roughly to a factor of 2; thus it seems that, with present 
knowledge of supernovae, and using the data from SN1987A, we can only safely exclude 



from [331]. 
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Figure 18: Examples of oscillation into sterile neutrinos allowed from present solar neu- 
trino data. In the plots are shown the 3 probabilities of oscillations P(z/e — > i^e) (decreasing 
red curve), P(i^e ^h,t) (increasing blue curve) and P(z/e — > v^) (lower black curve), as a 
function of the neutrino energy. The continuous (dotted) curve are the values during day 
(night). 

the occurrence of dramatic MSW effects. We remark that with a quantitative theory of 
core collapse supernovae, this test of sterile neutrinos would become very powerful. There 
are other interesting effects possibly related to heavy sterile neutrinos, such as r-process 
nucleosynthesis, re-heating of the shock, rocketing of pulsars, not discussed here. 

Big-bang nucleosynthesis and other cosmological probes An often stressed doubt 
concerning cosmology is about neglected or unknown effects. This said, we must recall 
that there has been impressive progress in the last years. The impact on neutrinos can be 
summarized as follows: (1) The number of neutrinos at the photon decoupling is bounded 
to be 3 ± 2. (2) The contribution to the energy density of the universe VLj^K^ is below or 
at about 1 %, or, equivalently, the sum of neutrino masses is below or at about 1 eV. (It 
should be recalled that there is an interesting, not universally accepted claim that the 
cosmology gives not a bound but a value for neutrino masses [344] . This testifies the high 
sensitivity reached by these methods and points to the interest in the value of the bias 
parameter cxg). (3) The effective number of neutrinos at nucleosynthesis time is 3±2, when 
extracted from deuterium abundance, or 2.4±0.7, when extracted from helium abundance. 
These numbers already imply strong bounds on sterile neutrinos, but do not rule out the 
sterile hypothesis for interesting regions of the parameters. One example is given by 
mirror neutrinos, when the new mass splittings are small enough. Another one is given 
by a new sterile state, that has only small mixing with the ordinary neutrinos. A second 
possibility is to have post-BBN phase transition involving the vacuum expectation value 
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of a light scalar field which mixes the active and sterile neutrinos so that at the time of 
BBN, the active and sterile neutrinos are unmixed. There are however strong constraints 
on the nature and interaction of the scalar field [94,345] and an interesting feature of these 
models is that they leave an imprint on the cosmic microwave background that can be 
tested in future experiments such as the Planck satellite mission. Yet another possibility 
is a large electron neutrino asymmetry in the early universe, which can compensate the 
effects of a number of extra neutrino types [346] . 

It has also been pointed out that the production of sterile neutrinos in the Early Uni- 
verse is strongly suppressed in cosmological scenarios for which the reheating temperature 
is as low as few McV [347]. In this case the bounds on the sterile neutrino parameters 
which can be derived from BBN and from the contribution of fli,h^ to the Dark Matter 
are much weaker than in the standard case. 

Ultra-high energy neutrinos Although there is a great interest in the search for ultra- 
high energy neutrinos, the number of reasonable (or even, less reasonable) mechanisms 
that have been discussed to produce them is not large. The reason is that neutrinos 
are produced along with electromagnetic radiation, that can be observed in a variety of 
ways, even when this is reprocessed. Following this line of thought, the astrophysical 
mechanism that can be conceived to overcome such a structure is the concept of a 'hidden 
source'. Another escape way from this constraint involves sterile neutrinos. Indeed, if 
there are ultra-high energy mirror neutrinos, they inevitably oscillate into neutrinos from 
our world on cosmic scales [348]. This scenario can provide intense fluxes of ultra-high 
energy neutrinos, subject only to the observable electromagnetic radiation from their 
interaction with the relic neutrino sea. 

7.1.3 Summciry of what we can leeirn on light sterile neutrinos 

In the view of many, there is something embarrassing in the hypothesis of (light) sterile 
neutrinos, and some of us believe that this is 'a solution searching for a problem'. However, 
history tells that the most prominent characteristic of neutrinos is that they are amazing! 
Said more seriously, we should certainly aim to measure neutrino properties, but we should 
not forget that we could make discoveries. And, when we think to new experiments, we 
should evaluate their potential for the investigation of sterile neutrinos. 

In this section, we recalled that light sterile neutrinos can play a role not only for LSND 
but also in terrestrial oscillations experiment, solar neutrinos, supcrnovae, astrophysics 
and cosmology. There are links between the various observables, but it is not impossible 
to conceive that sterile neutrinos have an important role only in astrophysics or cosmology 
(e.g., in core collapse supernovae or big-bang nucleosynthesis). More measurements and 
theoretical progresses will lead to important tests of the idea that the neutrinos that we 
know are not the full story. 

7.2 What can we learn about four-neutrino mass matrices 

The solar, atmospheric and LSND data require three different (mass)^ sphttings. These 
cannot be accommodated by three neutrino flavors, which provide only two independent 
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Am?. Additional degrees of freedom are necessary in order to understand all this data 
put together. The easiest option is to add a sterile neutrino and interpret the data in 
terms of oscillations of four neutrino flavors. 

The MiniBooNE experiment at Fermilab is crucial for confirming or refuting the LSND 
evidence for neutrino oscillations. If the LSND result is confirmed, a very exciting epoch 
in neutrino physics is just about to begin, as the number of questions that need to be 
answered becomes even larger than in the standard three-flavor case. 

A general 4-neutrino Majorana mass matrix is described by 4 masses, 6 mixing angles 
and 6 CP violating phases, 3 of which would affect oscillations. 

In this case there are 6 possible mass spectra, as shown in Fig. El These can be 
divided in two categories: "3+1" and "2+2". The "3+1" mass patterns are comprised of 
one sterile neutrino separated by Am-LSND from the other three. The group of three is the 
usual group of "active" neutrinos, one pair separated by AttIq and the third separated 
from these by Am\. The "2+2" patterns are comprised of two pairs of neutrinos, one 
separated by AitLq and the other by Am\. The two pairs are separated by AmLgNo- 

Both categories are already very strongly constrained by experiment [349,350]. 

In the "3+1" scenario the 3 states relevant for solar and atmospheric oscillations 
are mostly active and the forth state is almost entirely sterile. This pattern has the 
usual three active flavor scenario as a limiting case, so it agrees very well with all solar 
and atmospheric data. It is however harder to accommodate short baseline neutrinos 
oscillation experiments. This is related to the irony of the fact that LSND is an active 
flavor appearance experiment, showing i/^ — > z/g oscillations, but its solution has to involve 
an almost entirely sterile neutrino, while the other experiments remain unaffected by the 
presence of the sterile state. The bounds on sin^ 2^lsnd coming from KARMEN, CDHS, 
CHOOZ and atmospheric data are rather strong, almost conflicting with the value required 
to explain the LSND signal. The fit to all data is not very good, but the "3+1" scenario 
is not completely excluded at this point. 

In the '2+2" scenario both solar and atmospheric neutrino oscillations involve some 
fraction of conversion into a sterile state. This fractions are now very strongly constrained 
by the atmospheric, solar and reactor data, making the fit to the data to be rather poor 
in the'2+2" case. The global analysis are usually performed by considering three mixing 
angles and neglecting the other three, which are known to be small. Including these 
additional small angles might improve the quality of the fits, but the "2+2" scenario is 
strongly disfavored. 

The presence of the fourth, sterile neutrino also has implications in cosmology and 
cosmological observations impose further constraints on the allowed parameter space, as 
discussed in the previous section. 

The first question regarding four neutrino mass matrices, namely if they are indeed 
necessary to interpret the experimentally observed neutrino oscillation data will soon be 
answered by the MiniBooNE experiment. Assuming the answer is positive, a whole new 
set of questions arises. Just as in the three flavor case, one would like better measurements 
of all Am^'s and mixing angles. Given the much larger number of phases involved in the 
4 neutrino case, the possibilities for observing CP violation in the neutrino sector become 
very rich and maybe more easily accessible [351]. If the LSND signal is confirmed, than 
there must be a state with mass higher than a/ AmLg^D- This would be in the range of 
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sensitivity of future tritium (5 decay experiments like KATRIN, raising the possibility of 
determining the absolute scale of neutrino mass. For Majorana neutrinos, neutrinoless 
double beta decay might also be accessible [343,352]. By combining data from all types 
of experiments, the specific mass pattern could also be determined. 
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Figure 19: Four neutrino mass spectra. 



7.3 Heavy sterile neutrinos 

The general motivation for considering sterile neutrinos has already been discussed. Many 
extensions of the Standard Model imply the existence of more than one sterile neutrino 
with couplings to the active ones. The right-handed neutrinos participating in the seesaw 
mechanism, mirror neutrinos, the neutrinos in extra-dimensional models, and the right- 
handed technisinglet neutrinos in the ETC model discussed in other sections are some 
examples of such sterile neutrinos. 

"Light" (below ~ 10 eV) sterile neutrinos have been discussed in Sec. 17.11 Here we 
concentrate on "heavy" ones, by which we mean sterile neutrinos with masses above ~ 
10 eV, but below ~ IGeV. As noted, the mechanism constructed in Refs. [92,93,441] 
for light neutrinos in Technicolor theories leads to (two) heavy neutrino mass eigenstates 
in this range. We do not discuss here "very heavy" neutrinos (e.g. GUT scale), whose 
properties have been talked about above. Once sterile neutrinos are introduced, there 
is no definitive prediction for either the number or the masses of these light neutral 
fermions. Answering the question of the total number of neutrinos (active and sterile) 
and their masses is fundamental, as it would lead to much progress in understanding 
physics beyond the Standard Model. It is thus very important to address these issues 
from the experimental/observational point of view. 

Heavy sterile neutrinos with couplings to the active ones have profound implications in 
cosmology and astrophysics. They can also be constrained by several types of laboratory 
experiments. We discuss here the present status and future prospects for determining the 
properties of these heavy sterile neutrinos. 
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7.3.1 Laboratory Experiments 

The existence and mixing of heavy sterile neutrinos has many effects on particle and 
nuclear decays. These include contributions to /x"*" e"'"7, /i"*" e~^e~^e~, Kl — > /U^e^, 
and —>■ 7r"'"/x^e^, among others. However, given limits on mixing, these contributions 
are expected to be quite small (see further below on /x"*" e"'"7). 

The two-body leptonic decays of charged pions and kaons, and also measurement of 
the differential decay distribution in fi decay, can be used to search for, and set bounds on, 
the emission of massive neutrinos via lepton mixing [353,354]. The experimental signature 
for the emission of a massive neutrino via lepton mixing in tt^ or K^^ decay would be the 
appearance of an additional peak in the momentum spectrum of the outgoing charged 
lepton ^+ = or e+. The position of the extra peak is determined by the mass of the 
heavy neutrino and the size of the extra peak is proportional to the mixing |?7£/ip, where 

1 = e or n, between the extra state and the neutrino z/^. Initial bounds from retroactive 
data analyses were given in [353], and dedicated searches were carried out in 7r^2 [355], 
-^(u2 [356], Trjg [357], and [358] decays. Because of renewed interest [359], some recent 
searches in 7r^2 ^-^e reported in [360]. Resultant upper bounds on range down to 10~^ 
- 10-'^ in the mass range from several McV to --^ 300 McV. Admixed massive neutrinos 
also affect the observed ratio of branching ratios BR{t:^2) / ^^{''^^2) ^^^^ been 
used to set limits (e.g., [361]). 

The nuclear beta decay spectrum is also affected by the presence of heavy neutrinos 
mixed with v^- The spectrum would have a "kink" at the endpoint energy Emaxi'^h) 
[353,362-364]. The position of the kink determines the mass of the heavy state, m/j, and 
the change in the slope of the spectrum determines the mixing |f/e/iP- Many experiments 
searching for such kinks in the Kurie plots of nuclear beta decays were carried out in the 
1980's and beginning of the 1990's; although a few claimed positive results, these were 
refuted [3]. Resultant upper limits were of the order |C/e/iP ~ 10"^ for masses between 10 
keV and ~ 1 MeV. 

The nuclear transition involving electron capture, e~ + (Z, A) (Z — 1, A) + z/g [365] 
and muon capture, ji" + {Z, A) ^ {Z — 1, A) + z/^ [366] can also be used to search for, 
and put limits on, massive neutrino emission via mixing. 

Assuming that the additional sterile neutrinos are Majorana, there are several | AL| = 

2 transitions and meson and hyperon decays (processes analogous to neutrinoless nuclear 
double beta decay) that can occur. One of these is the nuclear transition fi~ + {A, Z)— > 
H'^ + {A, Z — 2) [367]. Meson decays include n^fi^fi^. A first upper limit on 
the branching ratio for this decay was set in [368], and this has been greatly improved 
by a dedicated search in a recent BNL experiment [369] (see also [370-373]). Analogous 
|AL| = 2 decays of heavy-quark mesons are also of interest. Early searches inchidc one 
by the Mark II detector at PEP for the decays K~ fi^fi^ and — > 7r~/i+yU+ [374] 
and one by the CLEO experiment at CESR for B'^ — > fi^fi^ [375]; current limits are 
given in Ref. [3]. One can also consider |AL| = 2 hyperon decays such as S~ — >• 

and E~ — > pjj,^ n~ . A first upper limit, on BR{E~ pfj,~fi~), was set in Ref. [376]; a 
recent dedicated search reporting a much improved limit is by the HyperCP experiment 
at Fermilab [377]. 

Mixing between heavy and light neutrinos also leads to neutrino decay [378-380]. 
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A number of searches for neutrino decays in accelerator experiments have been carried 
out [381]. Depending on the assumed mass of the neutrino mass eigenstate, various 
decays are possible, including v'e^e~ , — > u'li^e^, Uh u'li^ii". Bounds on 

various combinations of mixing angle factors from these experiments are reported in [381]; 
published limits range down to \Uai\^ < 10"^, i = e,/j, for heavy neutrino masses of 
several hundred McV [3,381]. In addition to weak charged current contributions to these 
neutrino decays, there are also contributions from the weak neutral current, since, in 
the presence of sterile neutrinos, the weak neutral current is not, in general, diagonal in 
mass eigenstates [40,41,382]. Present and future experiments as MiniBooNE and MINOS 
might have the possibility to improve on some of the present bounds on heavy neutrino 
decays [382]. Searches for heavy neutrino production and decay have also been carried 
out at e+e" colliders; see [3] for limits. For masses between a few GeV and mz the best 
bounds come from measurements at the Z pole, where the Z boson could decay into a 
standard neutrino and a sterile one, that would further decay [3] . 

7.3.2 Astrophysics and cosmology 

The existence of sterile neutrinos with even very small mixing to the active ones can have 
dramatic consequences in astrophysics and cosmology. These are discussed by a differ- 
ent working group [9]. Because astrophysical and cosmological observations provide the 
strongest constraints and prospects for future answers regarding heavy sterile neutrinos, 
we do, however, include here an overview of this subject. 
Cosmology 

Massive sterile neutrinos could be produced in the early universe and can provide 
some or even all the required dark matter. Heavy sterile neutrinos can be produced by 
scattering-induced conversion of active neutrinos [383]. These neutrinos are produced 
non-resonantly and they can be a warm dark matter candidate. A different mechanism 
of production of heavy sterile neutrinos appears if there is a non-vanishing initial lepton 
number in the Universe [384]. In this case sterile neutrinos can be produced resonantly 
and the energy spectrum is in this case highly non-thermal. The sterile neutrino can act 
then as a warm, cool or even cold dark matter [385] . 

Cosmological observations impose strong constraints on massive sterile neutrinos [385, 
386] . The radiative decay of such neutrinos to a light neutrino and a photon would affect 
the diffuse extragalactic background radiation, by producing a large number of photons of 
energy of the order rriH- The DEBRA experiment is now constraining the parameter space 
of sterile neutrinos based on this. The Chandra X-ray observatory has a great potential 
of resolving a considerable fraction of the observed X-ray background and consequently 
imposing much stronger constraints or potentially detecting X-ray fluxes from dark matter 
sterile neutrinos in the gravitational potential wells of clusters of galaxies. 

Heavy sterile neutrino decay prior to cosmic microwave background (CMB) decoupling 
increases the energy density in relativistic particles, leading to further constraints on the 
allowed parameter space. 

Big bang nucleosynthesis (BBN) is one of the big successes of the Standard Model 
of cosmology, successfully predicting the primordial abundance of light elements. The 
energy density in the sterile neutrino sea prior to BBN must not be too high in order not 
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to spoil the successful predictions of BBN. Photoproduction of deuterium and ^Li from 
decay of sterile neutrinos after BBN also imposes additional constraints on the sterile 
neutrino parameter space. 

Large scale structure observations are also essential, as they can constrain the nature 
of the dark matter (hot, warm or cold), consequently setting the scale for the mass of the 
sterile state. 

Cosmo logical constraints are illustrated in Fig. 1201 (from [385]). 




Figure 20: Cosmological/astrophysical constraints on heavy sterile neutrinos. 

Supernovae 

Neutrinos play a dominant role in core collapse supernovae. Even small admixtures of 
heavy sterile neutrinos can have profound implications for supernova physics. Too much 
active neutrino conversion into sterile states can lead to too much energy loss to sterile 
neutrinos, contradicting observations from supernova SN1987a [385,386]. The energy 
emitted in sterile neutrinos depends on the mixing angle between the sterile and active 
states in matter. For a long time it has been thought that most of the emission occurs in 
the resonant region, where the effective mixing angle becomes 7r/4. This is not necessarily 
true because of a non-linear effect that limits the resonant emission [385]. The effective 
matter potential for the z/g ^ Vg is given by 

V = Gpp/V2/mNi3Ye - 1 + + 2Y,^ + 2Y,^) (64) 

where Yi = {rii — nj)/nB. For antineutrinos the matter potential changes sign. Due to the 
presence of the Y^, terms, coming from neutral current neutrino-neutrino scattering, this 
effective potential has zero as a fixed point. Once a resonance is reached for, let's say, 
neutrinos, UgS start converting with maximal effective mixing angle into Ug. This decreases 
the Yj, term, driving the parameters off resonance and thus limiting the emission. The 
matter potential is effectively driven to zero on a time scale of less than a second and 
the emission continues with vacuum mixing angle. The parameter space where the sterile 
neutrino emission from supernovae could be relevant is also marked on Fig. 1201 

It is interesting to note that a sterile neutrino in the few keV region could also account 
for pulsar kicks [387]. In the presence of a strong magnetic field, neutrinos are emitted 
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asymmetrically in a supernova core. This asymmetry is washed out for active neutrinos 
which are trapped. If some conversion to sterile neutrinos occurs, these can escape the core 
of the star preserving the initial asymmetry. Only a few percent asymmetry is sufficient 
to account for the initial kick of the star that would explain the unusually high velocities 
of pulsars. 

To summarize, the presence of heavy sterile neutrinos with some (even very small) 
mixing to active neutrinos has numerous implications in astrophysics and cosmology. At 
present, a sterile neutrino with mass of the order of a few keV and very small mixing 
(sin^ 29 10~^) with an active one seems to be allowed by constraints, could account for 
all or some fraction of the dark matter in the Universe, would affect emission of supernovae 
neutrinos, could explain the pulsar kicks and might lead to observable contributions to 
X-ray measurements. In the future, more and more precise cosmological observations can 
constrain very strongly the parameter space for such sterile neutrinos, potentially closing 
the window or leading to detection of some signal. 
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8 Supersymmetry and neutrinos 



Neutrino masses are not the only motivation to extend the standard model. One also 
likes to extend it in order to solve the gauge hierarchy problem. Models of low-energy 
supersymmetry are attractive candidates for the theory of TcV scale physics. In the 
minimal supcrsymmctric extension of the Standard Model (MSSM) neutrinos arc massless. 
Thus, we need to consider supcrsymmctric extensions of the Standard Model that allow 
for neutrino masses. 

There are basically three questions we like to answer when we talk about the relations 
between supersymmetry and neutrinos 

1. Can successful predictions for neutrino masses of non-supersymmetric extensions 
of the Standard Model be retained once these models are supersymmetrized? In 
particular, can supersymmetry help in making such models more motivated? 

2. Are there models where neutrino masses arise only due to supersymmetry? 

3. Are there interesting phenomena in the slepton sector that can shed light on the 
issue of neutrino masses, lepton number violation and lepton flavor violation? 

The questions in the first item were already discussed in previous sections of this 
review. Generically, supersymmetry does not upset model predictions regarding neutrinos, 
and in some cases it in fact helps. For example, in the case of GUT, making the model 
supcrsymmctric helps to achieve coupling unification and thus to make the model realistic. 

In the following we concentrate on the last two items. We briefly describe two frame- 
works where neutrino masses are tightly connected to supersymmetry, or more precisely, 
to supersymmetry breaking, neutrino masses from R-parity violation and from supersym- 
metry breaking. We then discuss two effects, that of sneutrino oscillation and sneutrino 
flavor oscillation, that can help us disentangle the origin of neutrino masses using super- 
symmetric probes. The aspects of supcrsymmctric seesaw models, i.e. the connection to 
decays such as — > 67 have been discussed in previous Sections. 

8.1 Neutrino masses from R-parity violation 

Neutrino masses from R-parity violation have been extensively studied. Here we briefly 
summarize the main results. For a more complete reference list where more details can 
be found, see, for example, [388]. 

Once R-parity is violated there is no conserved quantum number that would distinguish 
between the down-type Higgs doublet and the lepton doublets. (For deflnitions and 
notation see, for example, [389].) Thus, these flelds in general mix. Such mixings generate 
neutrino masses; in fact, they generically produce too large masses. One neutrino gets a 
tree level mass which depends on the mixings between the Higgs and the sneutrinos [390] . 
The other two neutrinos get their masses at the one loop level, and thus their masses are 
smaller by, roughly, a loop factor. 

There are many different one loop contributions to the neutrino masses. (For a com- 
plete list see [391,392].) The "standard" diagrams are those that arise from the R-parity 
violating trilinear couplings in the superpotential, the so-called A and A' couplings. These 
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are the only contributions that are present in the supersymmetric hmit. Once supersym- 
metry breaking is included, there are many more contributions (in fact, also the tree level 
contribution is present only due to supersymmetry breaking). These contribution are 
likely to be much larger than the A and A' loop induced masses. The dominant diagrams 
arc likely to be those that arise due to bilinear couplings [388]. These are the couplings 
that mix the scalar components of the Higgs and the neutrino superficlds. The basic rea- 
son that the A and A' loop induced masses are very small is that they arc proportional to 
the small down-type quark or charged lepton Yukawa couplings. This suppression factor 
is absent in the bilinear induced masses. 

The most attractive feature of R-parity violation models of neutrino masses is that 
they naturally generate hierarchical neutrino masses with large mixing angles. This is 
due to the fact that only one neutrino gets a mass at tree level, while the other neutrinos 
only acquire loop induced masses. Numerically, however, the predicted mass hierarchy is 
in general somewhat too strong. Models with R-parity violation also predict observable 
lepton-number violating processes at possibly observable rates (e.g., [371,376,393]). 

The biggest puzzle posed by R-parity violation models is to understand the smallness 
of the neutrino masses. There must be a mechanism that generates very small R-parity 
violating couplings. There are several ideas of how to do it. For example, the small 
R-parity violation couplings can be a result of an Abelian horizontal symmetry [394] . 

8.2 Neutrino masses from supersymmetry breaking 

The smallness of neutrino masses can be directly related to the mechanism of supersym- 
metry breaking, in particular, to the mechanism that ensures a weak scale ^ parame- 
ter [231,232,395-398]. In general, there is no reason why the MSSM fi parameter is of the 
order of the weak scale. Generically, it is expected to be at the cut-off scale of the theory, 
say the Planck or the GUT scale. Phenomenologically, however, is required to be at 
the weak scale. One explanation, which is known as the Giudice-Masiero mechanism, is 
that a // term in the superpotential is not allowed by a global symmetry. The required 
effective weak scale n is generated due to supersymmetry breaking effects. 

The Giudice-Masiero mechanism can be generalized to generate small neutrino masses. 
It might be that the large Majorana mass term that drives the seesaw mechanism is 
forbidden by a global symmetry. Effective Majorana mass terms for the right-handed 
neutrinos, of the order of the weak scale, are generated due to supersymmetry breaking. 
The same global symmetry can also suppress the Dirac mass between the right and left- 
handed neutrinos. Then, the left-handed neutrinos have very small Majorana or Dirac 
masses as desired. 

The emerging neutrino spectrum depends on the exact form of the global symmetry 
that is used to implement the Giudice-Masiero mechanism. Nevertheless, the feature that 
the left-handed neutrino masses are very small is generic. 

8.3 Sneutrino oscillation 

As already discussed in this report, it is interesting to find out whether neutrinos have 
Majorana masses. In other words, we would like to find out if total lepton number is 
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violated in nature. As already discussed, the most promising way to do it, is to look 
for neutrinoless double beta decay. If supersymmetry is realized in nature we may have 
other probes. Here we describe one example, that of sneutrino-anti-sneutrino mixing and 
oscillation [399]. 

Consider a supersymmetric extension of an extended Standard Model that contains 
Majorana neutrino masses. For simplicity we consider only one neutrino generation. In 
such models, the effect of AL = 2 operators is to introduce a mass splitting and mixing 
into the sneutrino-anti-sneutrino system. This phenomena is analogous to the effect of 
a small AS = 2 perturbation to the leading AS = mass term in the S-system which 
results in a mass splitting between the heavy and light neutral B mesons. The very small 
mass splitting can be measured by observing flavor oscillations. The flavor is tagged in 
S-decays by the final state lepton charge. Since Atub ~ T^, there is time for the fiavor 
to oscillate before the meson decays. 

The sneutrino system can exhibit similar behavior. The lepton number is tagged in 
sneutrino decay using the charge of the outgoing lepton. The relevant scale is the sneutrino 
width. If the sneutrino mass splitting is large, namely when 

_ Anii^ 

xc = —— > 1, (65) 

and the sneutrino branching ratio into final states with a charged lepton is significant, then 
a measurable same sign dilepton signal is expected. Any observation of such oscillation 
will be an evidence for total lepton number violation, namely for Majorana neutrino 

masses. 

The size of the same sign lepton signal is model dependent. It is generically expected 
that the sneutrino mass splittings are of the order of the neutrino masses, as both are 
related to AL = 2 operators. The sneutrino width can vary a lot. When two body decay 
channels are dominant the width is roughly of the order of Fj; ~ aM^ ~ 0{1 GeV), and 
thus Xj> <^ 1 and a very small signal is expected. In models where the two body decay 
channels are forbidden the sneutrino width can be much smaller. For example, in models 
where the stau is the LSP there can be a situation where the sneutrino has only three 
body decay channels allowed. Then may be large enough for the oscillation signal to 
be observed. 



8.4 Sneutrino flavor oscillation 

As we know there are large mixing angles in the lepton sector. An independent probe of 
lepton flavor violation can be provided in supersymmetric models via the slepton sector 
[400-403]. While generally the slepton mixings are not directly related to the lepton 
mixings, both are lepton flavor violating effects. 

Slepton fiavor oscillations arise if the sleptons mass eigenstates are not fiavor eigen- 
states. Experimentally, the signal consists of a lepton fiavor imbalance. For example, 
in a linear collider one can search for a signal like e+e" vv* — * Xe"*"/!", where X is 
the non leptonic part of the final state. In hadron colliders, one can look for signals like 
X~ — > v\x~ — > Xe^[i~ ^ where x~ is a chargino. The oscillation probabilities depend on 
the slepton mass splittings and their mixing angles. In principle, oscillation signals can 
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be used to measure these parameters. Even if this cannot be achieved in practice, just the 
observation of an oscillation signal will provide an independent confirmation for lepton 
flavor non-conservation. 

In many proposals for supersymmetry breaking, a high degree of degeneracy among 
sleptons is predicted. As a result, there is the potential for large flavor mixing among the 
sleptons. This can lead to substantial and observable flavor violating signals at colliders. 
To be readily observable, it is necessary that mass splittings between the states are not 
much smaller than the decay widths, and that the mixing angles are not very small. In a 
large class of supersymmetry breaking models, the splittings can be comparable or even 
larger than the widths, and the mixing angles may be large. In such cases, dramatic 
collider signatures are possible. 
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9 Expectations in Superstring Constructions 



There has been relatively little work on the implications of superstring theories for neu- 
trino masses. However, it is known that some of the ingredients employed in grand uni- 
fied theories and other four-dimensional models may be difficult to implement in known 
types of constructions. For example, the chiral supcrmultiplets that survive in the ef- 
fective four-dimensional field theory are generally bi-fundamental in two of the gauge 
group factors (including the case of fundamental under one factor and charged under 
a U(l)) for lowest level heterotic constructions, or cither bi-fuudamental or adjoint for 
intersecting brane constructions. This makes it difficult to break the GUT symmetry, 
and even more so to find the high-dimcnsional Higgs representations (such as the 126 of 
5*0(10)) usually employed in GUT models for neutrino and other fermion masses. Thus, 
it may be difficult to directly embed many of the models, especially GUT models involv- 
ing high-dimensional representations rather than higher-dimensional operators, in a string 
framework. Perhaps more likely is that the underlying string theory breaks directly to an 
effective four-dimensional theory including the Standard Model and perhaps other group 
factors [404] . Some of the aspects of grand unification, especially in the gauge sector, may 
be maintained in such constructions. However, the GUT relations for Yukawa couplings 
are often not retained [405-407] because the matter multiplets of the effective theory may 
have a complicated origin in terms of the underlying string states. Another difference is 
that Yukawa couplings in string derived models may be absent due to symmetries in the 
underlying string construction, even though they are not forbidden by any obvious sym- 
metries of the four-dimensional theory, contrary to the assumptions in many non-string 
models. Finally, higher-dimensional operators, suppressed by inverse powers of the Planck 
scale, are common. 

Much activity on neutrino masses in string theory occurred following the first super- 
string revolution. In particular, a number of authors considered the implications of an Eq 
subgroup of the heterotic Eg x Es construction [405]- [410]. Assuming that the matter 
content of the effective theory involves three 27's, one can avoid neutrino masses alto- 
gether by fine-tuned assumptions concerning the Yukawa couplings [405]. However, it is 
difficult to implement a canonical type 1 seesaw. Each 27 contains two Standard Model 
singlets, which are candidates for right-handed neutrinos, and for a field which could gen- 
erate a large Majorana mass for the right-handed neutrinos if it acquires a large vacuum 
expectation value and has an appropriate trilinear coupling to the neutrinos. However, 
there are no such allowed trilinear couplings involving three 27's (this is a reflection of 
the fact that the 27 does not contain a 126 of the 5*0(10) subgroup). Eq string-inspired 
models were constructed to get around this problem by invoking additional fields not 
in the 27 [26,407] or higher-dimensional operators [410], typically leading to extended 
versions of the seesaw model involving fields with masses/ vevs at the TeV scale. 

Similarly, more recent heterotic and intersecting brane constructions, e.g., involving 
orbifolds and twisted sectors, may well have the necessary fields for a type 1 seesaw, but 
it is again required that the necessary Dirac Yukawa couplings and Majorana masses for 
the right-handed neutrinos be present simultaneously. Dirac couplings need not emerge 
at the renormalizable level, but can be of the form 

{S[---S',_,)NLHJM^I\ (66) 
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where the S'^ are standard model singlets which acquire large expectation values (d = 3 
corresponds to a renormalizable operator). Similarly, Majorana masses can be generated 
by the operators 

{Si---Sn-2)NN/M^^\ (67) 

Whether such couplings are present at the appropriate orders depends on the underlying 
string symmetries and selection rules, which are often very restrictive. It is also nec- 
essary for the relevant S and S' fields to acquire the needed large expectation values, 
presumably without breaking supersymmetry at a large scale. Possible mechanisms in- 
volve approximately flat directions of the potential, e.g., associated with an additional 
[/(!)' gauge symmetry [95,411], stringy threshold corrections [412,413], or hidden sector 
condensates [414]. 

There have been surprisingly few investigations of neutrino masses in explicit semi- 
realistic string constructions. It is difficult to obtain canonical Majorana masses in in- 
tersecting brane constructions [415], because there are no interactions involving the same 
intersection twice. Two detailed studies [416, 417] of nonsupersymmetric models with a 
low string scale concluded that lepton number was conserved, though a small Dirac mass 
might emerge from a large internal dimension. Large enough internal dimensions for the 
supersymmetric case may be difficult to achieve, at least for simple toroidal orbifolds. 

There are also difficulties for heterotic models. An early study of orbifolds yielded 
no canonical Dirac neutrino Yukawas [418] at low order. Detailed analyses of free fermionic 
models and their flat directions were carried out in [414,419] and [420,421]. Both studies 
concluded that small Majorana masses could be generated if one made some assump- 
tions about dynamics in the hidden sector. In [414, 419] the masses were associated 
with an extended seesaw involving a low mass scale. The seesaw found in [420, 421] 
was of the canonical type I type, but in detail it was rather different from GUT-type 
models. A seesaw was also claimed in a heterotic Z3 orbifold model with Eq breaking 
to SU{3) X SU{3) X SU{3) [422]. A recent study of Zq orbifold constructions found 
Majorana-type operators [423], but (to the order studied) the Si fields did not have the 
required expectation values. 

In [424] a large class of vacua of the bosonic orbifold were analyzed with emphasis 
on the neutrino sector to determine whether the minimal type I seesaw is common, or 
if not to find possible guidance to model building, and possibly to get clues concerning 
textures and mixing if examples were found. Several examples from each of 20 patterns of 
vacua were studied, and the nonzero superpotential terms through degree 9 determined. 
There were huge number of D flat directions, with the number reduced greatly by the F- 
flatness condition. Only two of the patterns had Majorana mass operators, while none had 
simultaneous Dirac operators of low enough degree to allow neutrino masses larger than 
10^^ eV. (One apparently successful model was ruined by off-diagonal Majorana mass 
terms.) It is not clear whether this failure to obtain a minimal seesaw is a feature of the 
particular class of construction, or whether it is suggesting that stringy constraints and 
selection rules might make string vacua with minimal seesaws rare. Systematic analyses 
of the neutrino sector of other classes of constructions would be very useful. 

There are other possibilities for obtaining small neutrino masses in string construc- 
tions, such as extended seesaws and small Dirac masses from higher dimension operators. 
Small Dirac neutrino masses from type I anisotropic compactiflcations have been dis- 
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cussed recently in [425]. The possibility of embedding type II seesaw ideas (involving 
Higgs triplets) in string constructions was considered in [426]. It is possible to obtain a 
Higgs triplet of SU (2) with non-zero hypercharge in a higher level construction (in which 
SU{2) X SU{2) is broken to a diagonal subgroup). In this case, because of the underlying 
SU{2) X SU (2) symmetry the Majorana mass matrix for the light neutrinos should involve 
only off-diagonal elements (often with one of the three off-diagonal elements small or van- 
ishing), with profound phenomenological consequences, including an inverted hierarchy 
and two large mixings. This is a top-down motivation for the texture C in Table |21 
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10 Theories with a TeV-scale 17(1)' 



Many extensions of the Standard Model and MSSM include the existence of additional 
non-anomalous 11(1)' gauge symmetries. These include many superstring constructions 
[427], grand unified theories, Little Higgs models, and models of dynamical symmetry 
breaking (DSB) [428]. In a regular grand unified theory the U{1)' breaking needs to be 
at a large scale, because scalars that can mediate proton decay can have masses no larger 
than the 11(1)' breaking scale. In string theories, the symmetry breaking is usually 

induced by soft supersymmetry breaking effects at the TeV scale [427,429,430], although 
in some cases there is the possibility of breaking along an F and D flat direction at an 
intermediate scale [411] (depending on the sign of a combination of soft mass-squares). 
The Little Higgs and DSB models are at the TeV scale. 

A TeV scale Z' has many interesting phenomenological consequences [431], but here we 
are concerned with neutrino masses. In the intermediate Z'-scale case, higher-dimensional 
operators, involving one or more powers of the fields with intermediate-scale vevs but sup- 
pressed by powers of the Planck mass, can yield naturally small Dirac neutrino masses [95, 
411]. Variants can also lead to mixing between light ordinary and sterile neutrinos, as 
suggested by the LSND results, or even to a type I seesaw. 

Models in which the U{1)' breaking is at the TeV scale generally do not allow a 
canonical type I seesaw model, because the Majorana mass for the right-handed neu- 
trino A^R requires U{1)' breaking (unless the carries no f/(l)' charge). However, a 
number of other possibilities are allowed for the neutrino masses [432], including small 
Dirac masses (e.g., associated with a second U{1)' broken at an intermediate scale), and 
Majorana masses associated with a TeV-scale seesaw [26] or a heavy Higgs triplet (Type 
II seesaw) [226]. The small Dirac mass case involves a strong constraint from Big Bang 
Nucleosynthesis, because the Z' interactions could efficiently produce the right-handed 
components prior to nucleosynthesis, leading to too much "^He [433]. For generic cou- 
phngs of the Z' to the Nr the observed abundance implies a Z' mass larger than around 4 
TeV, stronger than indirect or coUider constraints [434]. This can be evaded or weakened 
if the Nn carries no f/(l)' charge (as can occur naturally in some models involving two 
U{iy factors [431,432]) or if the mass is Majorana. 
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11 Neutrino Masses in Theories with Dynamical Sym- 
metry Breaking 

The source of electroweak symmetry breaking (EWSB) remains unknown, and a dy- 
namical origin of this breaking is an appeahng possibihty. This can be reahzed via the 
formation of a bihnear condensate involving fermions with a new strong gauge interac- 
tion, generically called Technicolor (TC) [435,436]. Indeed, one may recall that in both of 
the well-known two cases in which scalar fields have been used in models of spontaneous 
symmetry breaking, namely the Ginzburg-Landau effective Hamiltonian for superconduc- 
tivity and the Gell-Mann Levy sigma model for hadronic chiral symmetry breaking, the 
scalar fields were not fundamental, and the true underlying physics responsible for these 
respective phenomena involved the formation of bilinear fcrmion condensates (Cooper 
pairs in superconductivity and the (gg) condensate in QCD). In order to communicate 
this symmetry breaking in the Technicolor sector to the standard-model (technisinglet) 
fermions, one embeds the Technicolor model in a larger, extended Technicolor (ETC) the- 
ory [437,438]. To satisfy constraints from fiavor-changing neutral-current processes, the 
ETC vector bosons that mediate generation-changing transitions must have large masses. 
These masses arise from the sequential breaking of the ETC gauge symmetry on mass 
scales ranging from 10'^ TeV down to the TeV level. Precision measurements place tight 
constraints on these models, suggesting that there are a small number of new degrees of 
freedom at the TeV scale and that the Technicolor theory has a slowly running ( "walking" ) 
gauge coupling with large anomalous dimensions [439]. 

Since ETC models do not involve any superheavy GUT-scale mass, there was for a long 
time a puzzle of how one could explain light neutrino masses in these models. A possible 
solution to this puzzle was given in Ref. [92] and studied further in Refs. [93,441,442]. 
This does involve a seesaw, but one of a new type not involving any superheavy GUT 
scale. The resultant formula ~ [M^Y /rriR holds, with the largest Dirac neutrino 
masses of order a few keV and the relevant Majorana neutrino mass of order 0(0. 1) GeV 
to O(IOO) GeV. These Dirac and Majorana neutrino masses are greatly suppressed relative 
to conventional values. This suppression is a natural consequence of the representations 
of the ETC gauge group for the various neutrino fields. These ETC models are not yet 
developed sufficiently to make detailed predictions for leptonic mixing angles, but it seems 
possible to get substantial neutrino mixing. One interesting feature of this mechanism for 
neutrino masses is that there are only two, rather than three right-handed electroweak- 
singlet neutrinos, in contrast, e.g., to 50(10) GUT models. 

The ETC gauge group SU(A^etc) commutes with the Standard Model (SM) group 
GsM- The ETC group gauges the three generations of technisinglet fermions and connects 
them with the technicolored fermions. The ETC gauge symmetry is chiral, so that when 
it becomes strong, sequential breaking occurs naturally. The ETC symmetry breaking 
takes place in stages, leaving the residual exact technicolor gauge symmetry SU(A^t'c')- 
This entails the relation Netc — ^gen + Ntc = 3 -|- Ntc, where Ngen is the number of 
standard-model fermion generations. The choice of Ntc = 2 is required for the mechanism 
of Ref. [92] to work. This thus implies Netc = 5; i.e., one uses an SU (5) etc gauge theory. 
A related SU (5) etc theory had previously been studied in Ref. [440]. The choice Ntc = 2 
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has two other motivations: (a) it minimizes the TC contributions to the electroweak S 
parameter, which is a stringent constraint on TC theories, (b) with a standard-model 
family of technifermions, Ql = {^) j^-, Ll = (^)l' -^R' transforming according 

to the fundamental representation of SU(2)tc, it can yield an approximate infrared fixed 
point and the associated walking behavior. This sequential breaking of the SU{5)etc is 
driven by the condensation of SM-singlet fermions. 

One can explore whether this dynamical neutrino mass mechanism could take place 
in the context of ETC theories in which the strong-electroweak gauge group is extended 
beyond that of the Standard Model. Theories with the left- right symmetry group Glr — 
SU(3)cXSU(2)LxSU(2)flXU(l)B_L [72] and the group G422 = SU(4)psxSU(2)LxSU(2)ij 
[443] are of particular interest here, where B and L denote baryon and lepton number 
and SU{4:)ps combines color and lepton number. Ref. [93] presented a full ETC model 
in which these extended strong-electroweak gauge symmetries can be broken dynamically 
and showed that the mechanism of Ref. [92] can also hold here. Dynamical symme- 
try breaking of Glr has also been studied in Ref. [444]. Further, dynamical symmetry 
breaking of the electroweak symmetry can be triggered by a neutrino condensate [445]. 
ETC theories have many other testable implications. Some recent work has focused on 
ETC contributions to dimension-5 dipole moment operators [148, 149] and dimension-6 
four-fermion operators and their effects [442] . 
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12 Neutrinos in extra dimensions 



The pioneering idea by Kaluza and Klein (KK) [447] that our world may have more 
than four dimensions has attracted renewed interest over the last ten years [448-451]. 
The possible existence of extra dimensions has enriched dramatically our perspectives in 
searching for physics beyond the Standard Model. Obviously, extra dimensions have to 
be sufficiently compact to explain why they have escaped detection so far, although their 
allowed size is highly model-dependent [452]. This means that the derived constraints 
not only depend on the number of the fields sensitive to extra dimensions and their 
transformation properties with respect to those, but also on the geometry and/or the 
shape of the new dimensions. In the latter case, higher-dimensional theories may be 
distinguished between those formulated on a fiat space and those that utilize a warped 
geometry. 

Higher-dimensional theories may also provide interesting alternatives for explaining 
the smallness of the observed light neutrino masses. Their predictions for the light- 
neutrino spectrum can be confronted with recent neutrino oscillation data. In the follow- 
ing, we discuss a generic higher- dimensional neutrino scenario in which a fiat geometry is 
realized. 

The original formulation of higher-dimensional neutrino models [453,454] relies on the 
possible existence of singlet neutrinos that propagate in a higher [1 + (3 + 5)]-dimensional 
space which is usually termed the bulk, where 6 is the number of the additional spatial 
compact dimensions. In this formulation, the ordinary SM particles reside in a (1 + 
3) -dimensional Minkowski subspace, which is called the wall. Hence, the left-handed 
neutrinos and the Higgs bosons live on the wall. The overlap of their wave-functions with 
the bulk neutrinos is suppressed by the volume of the extra-dimensional space {R M f)^^"^ ~ 
Mp/Mp, where R is the common compactification radius, Mp is the fundamental gravity 
scale and Mp pa 10^^ TeV is the usual Planck mass. This volume-suppression factor gives 
rise to effective neutrino Yukawa couplings that are naturally very small, i.e. of order 
Mp/Mp ~ 10~^^, for Mp = 10 TeV, although the original higher-dimensional Yukawa 
couplings of the theory could be of order unity. This suppression mechanism [453,454] 
is a generic feature of these higher- dimensional neutrino models realized on an toroidal 
bulk; it has some dependence on the compactification properties of the bulk neutrinos 
and the structure of the Higgs sector of the theory [455,456]. 

To illuminate the discussion that follows, we consider a minimal extension of the SM 
where one 5-dimensional (bulk) sterile neutrino has been added. Furthermore, the extra 
dimension, say y, is compactified on a /Z2 orbifold. The SM fields are localized on a 
4-dimensional Minkowski subspace at the orbifold fixed point |/ = 0. Different minimal 
models may be constructed depending on the way that lepton number is broken: 

(i) One may add lepton- number violating bilinears of the Majorana type in the La- 
grangian [453], e.g. operators of the form N'^C^^^~^N, where C^^^ = —7173 is the 
charge conjugation operator. 

(ii) Lepton-number-violating mass terms can be generated through the Scherk-Schwartz 
mechanism [457]. This mechanism turns out to be equivalent to (i), after KK 
reduction. 
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(iii) Lepton number can be broken if the Z2-even and Z2-odd two-component spinors 
of the bulk neutrino couple simultaneously to the same left-handed charged lepton 
state. This is only possible if the 3-brane describing our observable world is shifted 
from the /Z2 orbifold fixed point [453,458]. 

(iv) Violation of the lepton number can be achieved by introducing operators of higher 
dimension in the number of fields, e.g. (L$)^/Mf. Such operators may be generated 
through gravity effects [459]. 

The current neutrino oscillation data provide an important test for singling out a 
good candidate model that includes higher-dimensional neutrinos. For example, orbifold 
models with one bulk neutrino [453,455,459-463], based on models of type (i) and/or 
(ii) mentioned above, seem to prefer the Small Mixing Angle (SMA) solution which is 
highly disfavored by recent neutrino data analyses. Alternatively, if all neutrino data are 
to be explained by oscillations of active neutrinos with a small admixture of a sterile KK 
component, then the compactification scale has to be much higher than the brane-Dirac 
mass terms. After integrating out the bulk neutrino of the model, the resulting effective 
light- neutrino mass matrix is of rank 1. Because of this restricted form of the neutrino 
mass matrix, two out of the three active neutrinos are massless. Since only one neutrino- 
mass difference can be formed in this case, it proves difficult to accommodate all neutrino 
oscillation data in these models [461-463]. 

12.1 Three bulk neutrinos 

One way to avoid this problem is to add three bulk neutrinos, one for each generation [460, 
463]. This model, in the absence of CP phases, is characterized by seven parameters: three 
neutrino masses mi 2,3 and three mixing angles for left handed neutrinos as defined earlier 
and the radius of the large extra dimension. Since the three mixing angles are arbitrary, 
the model can easily accommodate the bilarge mixing solution preferred by oscillation 
data. In the diagonal mass basis, the bulk neutrinos are associated with mass eigenstates. 
The mixing of the ith active neutrino with the nth KK mode of the corresponding bulk 
neutrino is given by S,i,n — rriiR/n. It is interesting that all mixings are intimately 
connected with the masses. There are limits on from laboratory data such as CHOOZ- 
Palo- Verde as well as from big bang nucleosynthesis [464] . BBN constraints for one extra 
dimension give ^| < 1.7 x 10~^{eV.R)^'^^ . For a hierarchical pattern, using ^ = m^R, one 
gets R < 0.03 eV~^. The bounds from neutrino data such as solar and atmospheric etc. 
are less stringent [463] and are roughly given by: R < 4 eV~^. 

Among the consequences of this model, two are especially interesting. Both of these 
concern the KK tower of sterile neutrinos. 

(i) In the presence of the infinite tower of states, the magnetic moment of the neutrino 
gets contribution from all the states [465]. For instance in the scattering of a neutrino 
of energy E ^ 10 MeV (corresponding to a reactor neutrino beam) the number of states 
contributing to the magnetic moment is given by [ER)"^ ~ 10^^. Since all the states add 
incoherently, the effective magnetic moment is increased from 10~^°/iB to 10~^^fiB {fJ-B is 
the Bohr magneton). The effect on the differential cross section da/dT, where T is the 
electron recoil energy has recently been calculated and given in Fig. 112.11 [466]. 
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Figure 21: The figures give the contribution of neutrino magnetic moment for the case of 
single Dirac neutrino, for two large extra dimensions (and a comparison between the two) 
to differential cross section ^ (where T is the electron recoil energy) for neutrino electron 
scattering and compares it with the case of one right-handed neutrino ( "Standard Model 
with one right-handed neutrino"). 



(ii) A second consequence of the existence of the KK tower of sterile neutrinos is the 
possibility that when neutrinos travel through dense matter there can be MSW reso- 
nances [455] and give rise to a dip pattern [455,467] in the neutrino survival probability 
corresponding to energies spaced hj E ^ Aml^^^j2V2GFNe (i.e. E,4:E,9E, . . .). The 

dip arises because typically the survival probability goes like e~^^^. Therefore at lower 
energies, there is more suppression which with increase in energy becomes less and less ef- 
fective. Also the resonance condition is not satisfied as energy increases, if it was satisfied 
at lower energies. For the solar neutrinos, such dip structure is quite pronounced [467]. 
In the hierarchical pattern for neutrino masses, this would correspond to ~ 10 MeV for 
densities comparable to solar core. The value of the energy clearly depends on the size 
of the extra dimensions, growing with R~^. This is a very interesting phenomenon which 
could be used to probe this class of extra dimension models. 

12.2 Lepton number breaking in the bulk 

In the three bulk neutrino picture, all the neutrinos are Dirac neutrinos since the model 
has an additional global B — L symmetry. An interesting possibility is the scenario 
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Observable 


he = = hr = h ^ 1 
Lower limit on Mp/h? [TeV] 
5=2 5=3 


5 = 6 


BR(/i ^ 67) 


75 43 


33 


BR(/i eee) 


250 230 


200 


BR(/i -Ti ^ e -Ti) 


380 320 


300 



Table 7: One-loop- level limits on Mp/h? from [468]. 



(iii) described above, where sizable lepton-number violation is induced by shifting the 
?/ = brane by an amount a ~ 1/Mw- In this scenario, the tree-level rank-1 form of 
the effective neutrino mass matrix can be significantly modified through lepton-number 
violating Yukawa terms, thus offering sufficient freedom to describe the neutrino oscillation 
data [458]. 

In addition to constraints from neutrino oscillation data, other experiments can also 
play an important role in constraining higher-dimensional neutrino models. Specifically, 
strong limits on Mp and the Yukawa couplings of the theory may be obtained from the 
non-observation of lepton-flavor-violating decays in muon and tau decays and from the 
absence of /i-e conversion in nuclei [468,469]. Table [7| gives a brief summary of these 
limits. These phenomenological constraints are complementary to those obtained from 
pure theoretical considerations, such as perturbative unitarity [470]. 

Another low-energy experiment of great importance is the neutrinoless double beta 
decay of a nucleus. The recently claimed experimental evidence [31] of an effective 
neutrino mass of order 0.4 eV (see however [32]), combined with information from solar 
and atmospheric neutrino data, restricts the admissible forms of the light-neutrino mass 
hierarchies in 4-dimensional models with 3 left-handed (active) neutrinos. The allowed 
scenarios contain either degenerate neutrinos or neutrinos that have an inverse mass hi- 
erarchy [471]. A positive interpretation of the recently reported Oi//3/3 signal [31] imposes 
additional constraints on model-building. For example, higher-dimensional models that 
utilize the shining mechanism from a distant brane [472] may accommodate an effective 
neutrino mass of 0.4 eV but also predict the emission of Majorons. On the other hand, 
5-dimensional models formulated on a warped geometric space [473] face difficulties to 
account for the observable excess in [31]. 

In the context of /Z2 orbifold models, one has to solve an additional theoretical 
problem. The resulting KK neutrinos group themselves into approximately degenerate 
pairs of opposite CP parities. Because of this, the lepton-number-violating KK-neutrino 
effects cancel each other leading to unobservably small predicted values for the Oi//3/3 
decay. These disastrous CP parity cancellation effects can be avoided by arranging the 
opposite CP parity KK neutrinos to couple to the bosons with unequal strength. 
This feature can naturally be implemented if the ?/ = wall displaced from one of the 
/Z2 orbifold fixed points by an amount of order 1/Mw. A unique prediction of such a 
model [458] is that the effective neutrino mass and the scale of the light neutrino masses 
can be completely decorrelated. 
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13 Other new physics and neutrinos 



13.1 New long range forces 

Long range forces in the context of particle physics originated with the ideas of Yang 
and Lee [474] and Okun [475] who proposed that gauging the baryon number or lepton 
number would give rise to a composition dependent long range force which could be 
tested in the classic Eotovos type experiments [476]. A special class of long range forces 
which distinguish between leptonic flavors have far reaching implications for the neutrino 
oscillations [477, 478] which may be used as a probe of such forces. 

The standard model Lagrangian is invariant under four global symmetries correspond- 
ing to the Baryon and three Lepton numbers {a = e, fi, r). Of these, only three com- 
binations [479] of lepton numbers (i) Lg — L^, (ii) Lg — or (iii) Lfj_ — Lr, can be gauged 
in an anomaly free way without extending the matter content. The existence of neutrino 
oscillations implies that these symmetries have to be broken but the relevant gauge bosons 
can still be light if the corresponding couplings are very weak. It is possible in this case 
to obtain light gauge boson induced forces having terrestrial range (e.g. the Sun-Earth 
distance) without invoking extremely low mass scales [477]. The exchange of such boson 
would induce matter effects in terrestrial, solar and atmospheric neutrino oscillations. For 
example, the electrons inside the sun generate a potential Vlr at the earth surface given 
by 

^ ^^-^^ ""^^ ( iF^) ' ^^^^ 

2 

where a = ^ corresponds to the gauge coupling of the — L^^^. symmetry , is the 
number of electrons inside the sun and -Res is the earth-sun distance ~ 7.6 x 10^^ GeV~^. 
The present bound on the Z-dependent force with range A ~ 10^^ cm is given [476] 
by a < 3.3 x 10~^°. Eq. then shows that the potential Vlr can introduce very 
significant matter-dependent effects in spite of the very strong bound on a. One can 
define a parameter 

^ _ 2EyVLR 
Aim? 

which measures the effect of the long range force in any given neutrino oscillation ex- 
periment. Given the terrestrial bound on a, one sees that ^ is given by ^atm ~ 27.4 in 
atmospheric or typical long-baseline experiments while it is given by S,soiar ~ 7.6 in the 
case of the solar or KamLAND type of experiments. In either case, the long range force 
would change the conventional oscillation analysis. Relatively large value of a suppresses 
the oscillations of the atmospheric neutrinos. The observed oscillations then can be used 
to put stronger constraints on a which were analyzed in [477]. One finds the improved 
90% CL bound . 

aet, < 5.5 X 10"^^ ; aer < 6.4 x 10"^^ , (69) 

in case of the Le — L^^^ symmetries, respectively. 

Although these bounds represent considerable improvement over the conventional fifth 
force bound, they still allow interesting effects which can be used as a probe of such long 
range forces in future long-baseline experiments with super beam or at neutrino factories. 
As a concrete example, let us consider the influence of the Lg — gauge interactions 
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on the long-baseline oscillations of muon neutrinos of 0{GeV) energy. The oscillations of 
these neutrinos are governed by the following 3x3 (mass)^ matrix in the flavor basis: 

= f/*Diag(m?, ml ml)U^ + Diag(Acc + A^r, -A^r, 0) . (70) 

U denotes the (vacuum) mixing matrix for which we adopt the conventional parametriza- 
tion. Acc = 2E^V2GFne ^ (1.04x 10"^^ eV)2E^ describes the conventional MSW matter 
contribution generated by the earth matter (density p ~ 2.8 gm/cm^; electron fraction 
Ye ~ 0.49). The Alr ^ (1.04 x 10"^^ eV)2E,,a52 with 052 denoting the coupling of the 
long range force measured in units of 10"^^. The Alr term dominates over Acc if sat- 
urates the bound in Eq. (jU^ . The matter induced terms in Eq. (fTOj) modify the neutrino 
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Figure 22: The long-baseline neutrino oscillation probability Pen in case of vacuum (solid), 
earth matter effects (dotted) and with inclusion of the long range potential Vlr (dashed). The 
plotted curves correspond to a baseline of 740 km, — 2-5 x 10~^ eV^, Am^i = 7.0 x 10~^ 

eV^, (612,023) = (32°, 45°) , Oe/. = 5.5 x 10"^^ ^nd sin6li3 = 0.05. 

oscillation in a non-trivial manner. This effect is analyzed in the limit of the vanishing 
solar scale and Alr = in [480]. The 23 mixing angle remains unaffected by the matter 
induced contribution but the 13 angle can get resonantly enhanced for a neutrino energy 
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given by 

^ cos 201'^ Amor, ^ , , 

^ ^ ^ 11.8 GeV . 71 

2V2GFne 

This leads to a rise in the oscillation probability Pen, Figure (j^^ . The additional long 
range contribution results in a noticeable shift in the resonance energy as seen from Fig. 
()22|) . We assumed normal hierarchy in this figure. The resonance behavior would be 
absent in case of the inverted hierarchy or in case of the anti-neutrino beam. While more 
detailed study is required to distinguish these cases, it is clear that future observations of 
matter effects in the long-baseline neutrino experiments provide a good probe of additional 
long range forces. 

13.2 Lorentz noninvariance, CPT violation and decoherence 
13.2.1 CPT Violation 

In this section, we discuss neutrino oscillation phenomenology in the presence of CPT 
violation. CPT is a symmetry in any theory that satisfies the three assumptions that are 
normally taken for granted: (1) locality, (2) Lorentz invariance, and (3) hermiticity of 
the Hamiltonian. In particular, it predicts that the mass is common for a particle and 
its anti-particle. Any violation of CPT therefore would have a profound consequence on 
fundamental physics. 

The best limit on CPT violation is in the neutral kaon system, \m{K^) — m{K )| < 
lO'^^rrix = 0.50 x 10^^^ GeV [481]. Having such a stringent bound does not seem a 
sizable CPT violation in neutrino at the first sight. However, the kinematic parameter 
is mass-squared instead of mass, and the constraint may naturally be considered on the 
CPT-violating difference in mass- squared \m^{K^) -m'^{K )\ < 0.25 eV^. In comparison, 
the combination of SNO and KamLAND data leads to the constraint |Am^ — Am?! < 
1.3 X 10-2 eV^ (90% CL) and hence currently the best limit on CPT violation [482]. 

Having seen that the CPT violation in neutrino masses may be of size relevant to 
neutrino oscillation, it is useful to discuss how it may affect the phenomenology. In fact, 
the primary motivation for recent discussions on CPT violation in neutrino oscillation 
has been to reconcile LSND data with other data [483]. It is well known that the LSND 
data is not quite consistent with the other oscillation evidence and limits even if a sterile 
neutrino state is introduced, both for 2-1-2 and 3-1-1 spectra (see [349] for a recent analysis; 
adding more sterile states helps [47]). The main point is that the LSND oscillation 
is primarily an anti- neutrino oscillation z/^ z/g, while the solar neutrino oscillation 
is purely neutrinos v^^t-. It was shown to fit the LSND, solar, and atmospheric 

neutrino data simultaneously without invoking a sterile neutrino at that time [483-485]. 
Phenomenology had been further developed in Refs. [486,487]. 

However, KamLAND data shows Ve ^p.,T oscillation with parameters consistent with 
the solar neutrino oscillation, and the CPT-violation alone cannot explain LSND. A new 
proposal tried to explain LSND and atmospheric anti-neutrino oscillations with a single 
Am^ [487], which was excluded by a global fit in [488]. Currently the best fit to the data 
is obtained by allowing for one sterile neutrino and CPT violation [489]. Because the 
short-baseline experiments that are constraining the interpretation of LSND data with 
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Figure 23: The proposed spectra of neutrinos and anti-neutrinos in [483] and [484]. Ex- 
cluded by KamLAND. 



sterile neutrino are mostly in neutrinos but not in anti-neutrinos, the 3-1-1 spectrum is 
allowed if there is little mixing of the sterile state with others in neutrinos. 




Figure 24: The revised proposal in [487]. Excluded by the analysis taken from [488]. 

Even though arbitrarily changing the neutrino and anti-neutrino masses seems to pre- 
serve Lorentz invariance, interacting theory also violates Lorentz invariance [490]. All 
discussions above assumed Lorentz invariance and hence should be regarded as purely 
phenomenological exercise. One theoretically well-defined way to break CPT is to in- 
troduce a cosmological "matter effect," namely a background number density coupled 
to neutrinos. However, such framework does not explain data consistently [491]. See 
also [492] for a different framework of CPT violation and a recent discussion on the use of 
decoherence and CPT violation [493] . Lorentz invariance violation in the neutrino sector 
can arise via the see-saw mechanism. As discussed in [494] this could explain why it would 
not be seen in the charged lepton sector. 
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Figure 25: The revised proposal in [489] that combines CPT violation and a sterile neu- 
trino. The neutrinos always have 2 + 2 spectrum, while the anti-neutrinos may have either 
3 + 1 or 2 + 2 spectrum. 



13.2.2 Decoherence 

So far, CPT violation as an inequality of masses between particles and antiparticles was 
the only way we understood CPT violation in high energy physics. However, is this CPT 
violation the only way a violation of this symmetry can manifest itself in nature? Such a 
question becomes extremely relevant for the case of LSND, because it is possible that other 
mechanisms leading to CPT violation exist, unrelated, in principle, to mass differences 
between particles and antiparticles. Such additional mechanisms for CPT violation may 
well be capable of explaining the LSND results within a three generation scenario without 
invoking a sterile neutrino (a scenario which, on the other hand, is getting totally excluded 
as new experimental data become available). It is therefore necessary to explore whether 
alternative ways exist to account for the LSND result without invoking extra (sterile) 
neutrino states. 

Quantum decoherence is the key to answer this question. Indeed, quantum decoher- 
ence in matter propagation occurs when the matter subsystem interacts with an 'environ- 
ment', according to the rules of open-system quantum mechanics. At a fundamental level, 
such a decoherence may be the result of propagation of matter in quantum gravity space- 
time backgrounds with 'fuzzy' properties, which may be responsible for violation of CPT 
in a way not necessarily related to mass differences between particles and antiparticles. 

A characteristic example of such a violation occurs in quantum gravity models that 
involve singular space-time configurations, integrated over in a path integral formalism, 
which are such that the axioms of quantum field theory, as well as conventional quantum 
mechanical behavior, cannot be maintained. Such configurations consist of wormholes, 
microscopic (Planck size) black holes, and other topologically non-trivial solitonic objects, 
such as geons etc. Collectively, we may call such configurations space time foam. 

It has been argued that, as result, a mixed state description must be used {QG-induced 
decoherence) [495], given that such objects cannot be accessible to low-energy observers, 
and as such must be traced over in an effective field theory context. As a consequence of 
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that CPT invariance in its strong form must be abandoned in a foamy quantum gravity 
theory. Such a breakdown of CPT symmetry is a fundamental one, and, in particular, 
implies that a proper CPT operator may be ill defined in such QG decoherence cases. 

Some caution should be paid regarding CPT violation through decoherence. Prom a 
formal view point, the non-invertibility of the S-matrix, which implies a strong violation 
of CPT, does not preclude a softer form of CPT invariance, in the sense that any strong 
form of CPT violation does not necessarily have to show up in any single experimental 
measurement. This implies that, despite the general evolution of pure to mixed states, it 
may still be possible in the laboratory to ensure that the system evolves from an initial 
pure state to a single final state, and that the weak form of CPT invariance is manifested 
through the equality of probabilities between these states. If this is the case, then the 
decoherence-induced CPT violation will not show up in any experimental measurement. 

In the parametrization of [495] for the decoherence effects, one uses three decoher- 
ence parameters with dimensions of energy, a, /3, 7, where the positivity of p, required 
by the fact that its diagonal elements express probability densities, implies a.j > 0, 
and 0:7 > If the requirement of a completely positive map p{t) is imposed in 

the two generation case, then £ becomes diagonal, with only one non vanishing entry 
occupied by the decoherence parameter 7 > [496]. Following this approach, for a 
three generation scenario, we will assume for the 9x9 decoherence matrix C: [jCn^,] — 
Diag (0, —71, —72, —73, —74, —75, —76, —77, —78) in direct analogy with the two-level case 
of complete positivity [496,497], although there is no strong physical motivation behind 
such restricted forms of decoherence. This assumption, however, leads to the simplest 
possible decoherence models, and, for our phenomenological purposes in this work, we 
will assume the above form, which we will use to fit all the available neutrino data. It 
must be clear to the reader though, that such a simplification, if proved to be successful 
(which, as we shall argue below, is the case here), just adds more in favor of decoherence 
models, given the restricted number of available parameters for the fit in this case. 

In order to check these models, we have performed a comparison (as opposed to a 

fit) to SuperKamiokande sub-GeV and multi GeV data (40 data points) , CHOOZ data 
(15 data points) and LSND (1 datum), for a sample point in the vast parameter space of 
our extremely simplified version of decoherence models. Let us emphasize that we have 
not performed a x^-fit and therefore the point we arc selecting (by "eye" and not by x) 
is not optimized to give the best fit to the existing data. Instead, it must be regarded 
as one among the many equally good sons in this family of solutions, being extremely 
possible to find a better fitting one through a complete (and highly time consuming) scan 
over the whole parameter space. 

Cutting the long story short, and to make the analysis easier, we have set all the 
7j in the neutrino sector to zero, restricting this way, all the decoherence effects to the 
antineutrino one where we have assumed for the sake of simplicity, 7r = ^=2-10~^^- 
E and ^ — ^ — 1 ■ 10~'^/E , where E is the neutrino energy, and barred quantities 
refer to the antineutrinos, given that decoherence takes place only in this sector in our 
model. All the other parameters are assumed to be zero. All in all, we have introduced 
only two new parameters, two new degrees of freedom, 77 and 73, and we shall try to 
explain with them all the available experimental data. (For further details we refer the 
reader to [493]). 
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At this point it is important to stress that the inclusion of two new degrees of freedom is 
not sufficient to guarantee that one will indeed be able to account for all the experimental 
observations. We have to keep in mind that, in no-decoherence situations, the addition 
of a sterile neutrino (which comes along with four new degrees of freedom -excluding 
the possibility of CP violating phases) did not seem to be sufficient for matching all the 
available experimental data, at least in CPT conserving situations. 

In order to test our model we have calculated the of the 56 data points mentioned 
above for different scenarios. The results with which we hope all our claims become 
crystal clear are summarized in the table, were we present the comparison for the 
following cases: (a) pure decoherence in the antineutrino sector, (b) pure decoherence 
in both sectors, (c) mixing plus decoherence in the antineutrino sector, (d) mixing plus 
decoherence in both sectors, and (e) mixing only - the standard scenario: 



Model 


without LSND 


including LSND 


(a) 


980.7 


980.8 


(b) 


979.8 


980.0 


(c) 


52.2 


52.3 


(cl) 


51.1 


51.6 


(e) 


53.9 


60.7 



From the table it becomes clear that the mixing plus decoherence scenario in the an- 
tineutrino sector can easily account for all the available experimental information, includ- 
ing LSND. It is important to stress once more that our sample point was not obtained 
through a scan over all the parameter space, but by an educated guess, and therefore 
plenty of room is left for improvements. On the other hand, for the mixing-only/no- 
decoherence scenario, we have taken the best fit values of the state of the art analysis and 
therefore no significant improvements are expected. At this point a word of warning is in 
order: although superficially it seems that scenario (d), decoherence plus mixing in both 
sectors, provides an equally good fit, one should remember that including decoherence ef- 
fects in the neutrino sector can have undesirable effects in solar neutrinos, especially due 
to the fact that decoherence effects are weighted by the distance traveled by the neutrino, 
something that may lead to seizable (not observed!) effects in the solar case. 

One might wonder then, whether decohering effects, which affect the antineutrino sec- 
tor sufficiently to account for the LSND result, have any impact on the solar-neutrino 
related parameters, measured through antineutrinos in the KamLAND experiment. In 
order to answer this question, it will be sufficient to calculate the electron survival proba- 
bility for KamLAND in our model, which turns out to be Pp^^p^ iKamLANo— -63, in perfect 
agreement with observations. It is also interesting to notice that in our model, the LSND 
effect is not given by the phase inside the oscillation term (which is proportional to the 
solar mass difference) but rather by the decoherence factor multiplying the oscillation 
term. Therefore the tension between LSND and KARMEN data is naturally eliminated, 
because the difference in length leads to an exponential suppression. 

Having said that, it is now clear that decoherence models (once neutrino mixing is 
taken into account) are the best (and arguably the only) way to explain all the observations 
including the LSND result. This scenario, which makes dramatic predictions for the 
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upcoming neutrino experiments, expresses a strong observable form of CPT violation in 
the laboratory, and in this sense, our fit gives a clear answer to the question asked in 
the introduction as to whether the weak form of CPT invariance is violated in Nature. It 
seems that, in order to account for the LSND results, we should invoke such a decoherence- 
induced CPT violation, which however is independent of any mass differences between 
particles and antiparticles. 

This CPT violating pattern, with equal mass spectra for neutrinos and antineutri- 
nos, will have dramatic signatures in future neutrino oscillation experiments. The most 
striking consequence will be seen in MiniBooNE, according to our picture, MiniBooNE 
will be able to confirm LSND only when running in the antineutrino mode and not in 
the neutrino one, as decoherence effects live only in the former. Smaller but experimen- 
tally accessible signatures will be seen also in MINOS, by comparing conjugated channels 
(most noticeably, the muon survival probability). Higher energy neutrino beams or long- 
baseline experiments, will have significant deviations from the non-decoherence models, 
as our effects scale with energy and distance traveled, being therefore the best tool to 
explore decoherence models. 

If the neutrino masses are actually related to decoherence as a result of quantum grav- 
ity, this may have far reaching consequences for our understanding of the Early stages 
of our Universe, and even the issue of Dark Energy that came up recently as a result 
of astrophysical observations on a current acceleration of the Universe from either dis- 
tant supernovae data or measurements on Cosmic Microwave Background temperature 
fluctuations from the WMAP satellite. Indeed, decoherence implies an absence of a well- 
defined scattering S-matrix, which in turn would imply CPT violation in the strong form. 
A positive cosmological constant A > will also lead to an ill definition of an S-matrix, 
precisely due to the existence, in such a case, of an asymptotic-future dc Sitter (infiation- 
ary) phase of the universe, with Hubble parameter ~ implying the existence of a 
cosmic (Hubble) horizon. This in turn will prevent a proper definition of pure asymptotic 
states. 

We would hke to point out at this stage that the claimed value of the dark en- 
ergy density component of the (four-dimensional) Universe today, A ~ 10^^^^ M|,, with 
Mp ~ 10^^ GeV (the Planck mass scale), can actually be accounted for (in an amusing 
coincidence?) by the scale of the neutrino mass differences used in order to explain the 
oscillation experiments. Indeed, A ~ [{Ara^f/Mp^M^ ~ 10"^^^ M|, for Aw? ~ 10"^ eV^, 
the order of magnitude of the solar neutrino mass difference assumed in oscillation ex- 
periments (which is the one that encompasses the decoherence effects). The quantum 
decoherence origin of this mass then would be in perfect agreement with the decoherence 
properties of the cosmological constant vacuum, mentioned previously. 
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14 NuTeV Physics 



The NuTeV experiment [498] at Fermilab has measured the ratios of neutral to charged 
current events in muon (anti)neutrino-nucleon scattering: 

(t[v^N 12+X) r 



(73) 



where 

aju^N ^ /i+X) 1 
a{u^N ^ /i-X) ^ 2 ' 

and has determined the parameters g1 and gj^ [499] to be 

gl = 0.30005 ±0.00137, 

gl = 0.03076 ± 0.00110 . (74) 

The Standard Model (SM) predictions of these parameters based on a global fit to non- 
NuTeV data, cited as [c/|]sm = 0.3042 and [^|]sm = 0.0301 in Ref. [498], differ from 
the NuTeV result by 3a in gj^. Alternatively, if the SM is fit to the NuTeV result, the 
preferred range of the Higgs mass is 660 GeV < niH (90% C.L.) [500], well above the 
value of rriH ~ 90 GeV preferred by the non-NuTeV global fit [501]. 

The significance of the NuTeV result remains controversial [502] , and a critical exami- 
nation of the initial analysis is ongoing. Several groups are evaluating potential theoretical 
uncertainties arising from purely Standard Model physics which might be comparable to 
or larger than the quoted experimental uncertainty of the NuTeV result. Candidate 
sources of large theoretical uncertainty include next-to-leading-order (NLO) QCD cor- 
rections [503], NLO electroweak corrections [504], and parton distribution functions (es- 
pecially as involves assumptions about sea-quark asymmetries) [505]. The effect of the 
former has been estimated to be comparable in size to the NuTeV experimental uncer- 
tainty, while the latter two might give rise to effects comparable in size to the full NuTeV 
discrepancy with the Standard Model. Elucidation of the actual impact of these effects 
on the NuTeV result awaits a reanalysis of the NuTeV data. However, it remains a dis- 
tinct possibility that the discrepancy with the Standard Model prediction is genuine and 
that its resolution lies in physics beyond the Standard Model. Indeed, as Chanowitz has 
emphasized [500], the precision electroweak data indicate new physics whether anoma- 
lous data are excluded from global fits (since the preferred Higgs mass is then well below 
the direct search limit) or included in the fits (in which case anomalous data themselves 
demand a new physics explanation). 

Note that the NuTeV value for g^ in Eq. ()74|) is smaller than its SM prediction. This 
is a reflection of the fact that the ratios Ri, and Rp were smaller than expected by the SM. 
(The gj^ term is smaller than the gj^ term by an order of magnitude and is insignificant.) 
Thus, possible new physics explanations of the NuTeV anomaly would be those that 
suppress the neutral current cross sections over the charged current cross sections, or 
enhance the charged current cross sections over the neutral current cross sections. Two 
classes of models have been proposed which accomplish this task. 
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The first class comprises models which suppress R^y and Rp with the introduction of 
new neutrino-quark interactions, mediated by leptoquarks or extra U{1) gauge bosons 
(Z"s), which interfere either destructively with the Z-exchange amplitude, or construc- 
tively with the 1^-exchange amphtude [502] . In order to preserve the excellent agreement 
between the SM and non-NuTeV data, the new interactions must selectively interfere 
with the u^N ip^N) scattering process, but little else. This severely restricts the types of 
interactions that may be introduced. 

Ref. [502] proposes a model in which the Z' couples to B — 3L^. This model must be 
fine-tuned to avoid Z-Z' mixing [506] which would disrupt, among other things, lepton 
universahty at the Z-pole. Fitting the NuTeV anomaly requires 

^ 3 TeV . (75) 
9z' 

Bounds from direct Z' searches at the Tevatron and LEP limit the possible range of Mz' 
to Mz' > 600 GeV for gz' ~ 1, or 2 GeV < Mz' < lOGeV for gz' - 10"=^. 

The Z' in the model proposed in Ref. [507] does not couple the neutrinos and quarks 
directly, since the gauged charge is — L^. Rather, it is a tunable Z-Z' mixing in 
the model which is responsible for suppressing the neutral channel cross section. The 
same mixing violates lepton universality on the Z-pole and prevents the mechanism from 
completely mitigating the NuTeV anomaly. Z' masses in the range 60 GeV < Mz' < 
72 GeV, or Mz' > 178 GeV brings the theoretical value of gl within 1.6a of the NuTeV 
value while keeping lepton universality violation within 2a. 

In general, models in this class are constrained strongly by lepton universality, be- 
cause ui is the SU{2)l partner of New interactions which respect the SU{2)l gauge 
symmetry cannot affect neutrino couplings alone: they necessarily affect couplings of the 
charged leptons. Nevertheless, they provide possible explanations of the NuTeV anomaly, 
and predict a fiavor-selective gauge boson in the several 100 GeV to TeV range, well 
within reach of the LHC 

Models of the second class suppress the Zuu coupling by mixing the neutrino with 
heavy gauge singlet states (neutrissimos, i.e. right-handed neutrinos) [508-511]. For in- 
stance, if the SU (2) l active is a linear combination of two mass eigenstates with mixing 
angle 9, 

Vf, = (cos6')t'iight + (sin6')z/heavy , (76) 

then the Zv^v^ coupling is suppressed by a factor of cos^ Q (assuming the heavy states 
are too massive to be created on-shell). Likewise, the WfiPf^ coupling is suppressed by 
cos 6*. Although both the numerators and denominators of Ri, and Rp are suppressed in 
such a model, the suppression of the numerators exceeds that of the denominators, and 
the ratios are therefore diminished. More generally, if the ZuiUi coupling {i — e, /i, r) is 
suppressed by a factor of (1 — 6^), then the W£ui coupling is suppressed by (1 — £^/2), 
and Ru and Rp are suppressed by (1 — e^). 

The effect of such suppressions of the neutrino-gauge couplings is not limited to NuTeV 
observables alone. In addition to the obvious suppression of the Z invisible width by a 
factor of [1 — (2/3) (sg + + £t)], all SM observables will be affected through the Fermi 
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constant Gp which is no longer equal to the muon decay constant G^: 

Gf = (l + . (77) 

This shift in Gp will destroy the excellent agreement between the SM and Z-pole ob- 
servables. However, since Gp always appears in the combination pGp in neutral current 
amplitudes, the agreement can be recovered by absorbing the shift in Gp into a shift in 
p, or equivalcntly, in the oblique correction parameter T [512]. Indeed, it was shown in 
Ref. [510], that the Z-pole, NuTeV, and W mass data can all be fit with the oblique correc- 
tion parameters S, T, U, and a flavor universal suppression parameter s = = = Sr, 
the best fit values given by 

S = -0.03 ±0.10, 
T = -0.44 ±0.15, 
U = 0.62 ±0.16, 

e = 0.0030 ± 0.0010 , (78) 

for a reference SM with ttih = 115 GeV. Therefore, for this class of models to work, 
neutrino mixing with heavy gauge singlet states must be accompanied by new physics 
contributions to S, T, and U. The values of S and T can be accommodated within the 
SM by simply increasing the Higgs mass to hundreds of GeV, but the W mass requires 
a large and positive U parameter which cannot be generated within the SM. Thus, the 
models are not complete until some mechanism is found which explains the W mass. But 
then, if the SM is fit to the W mass alone, the preferred Higgs mass is far below direct 
search limits [500], which could be an indication that the W mass requires new physics 
regardless of NuTeV. 

At first blush, the preferred value of e above is also problematic. This implies a large 
mixing angle, 6 = 0.055 ± 0.010, if interpreted as due to mixing with a single heavy state. 
The commonly accepted seesaw mechanism [67-70] relates the mixing angle to the ratio 
of the neutrino masses: 

^ e'' . (79) 

^hcavy 

Choosing mught ~ 0.1 eV and mheavy ^ 100 GeV (mheavy > Mz is needed to suppress 
Finv) we find the mixing angle orders of magnitude too small: 6 ~ 10~^. However, this 
result does not mean that it is impossible to have a large enough mixing angle between 
the light and heavy states. As pointed out in Ref. [509], in models with more than one 
generation, the generic mass matrix includes enough degrees of freedom to allow us to 
adjust all the masses and mixings independently. Concrete examples of models with 
large mass hierarchies AND large mixing angles can be found in Refs. [511,513]. What 
is sacrificed, however, is the traditional seesaw explanation of the small neutrino mass: 
i.e. since the Majorana mass M in the neutrino mass matrix should be of the order of 
the GUT scale, the neutrino mass mught ~ rr? jM is naturally suppressed if the Dirac 
mass m is comparable to that of the other fermions. An alternative mechanism is used in 
Ref. [511]. There, an intergenerational symmetry is imposed on the neutrino mass texture 
which reduces its rank, generating naturally light (massless) mass eigenstates. 
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Abandoning the seesaw mechanism also frees the masses of the heavy states from being 
fixed at the GUT scale. Indeed, in the model discussed in Ref. [511], the assumption that 
neutrinos and up-type quarks have a common Dirac mass implies that the masses of the 
heavy state could be a few TeV, well within the reach of the LHC. Without quark-lepton 
unification mheavy could be even lighter, rendering them accessible to Tcvatron Run II. 

Because of the large mixing angles between the light and heavy states in this class of 
models, flavor changing processes mediated by the heavy states may be greatly enhanced 
[40,511,513,514]. As a result, stringent constraints can be placed on the models from the 
experimental limits on — > 67, r /I'j, [145] /i-e conversion in nuclei [515,516], muonium- 
antimuonium oscillation [517,518], etc. For instance, the MEGA limit on — > 67 leads 
to the constraint [511] 

£e£/. ~ . (80) 

Therefore, lepton universality among the ee must be broken maximally. Ref. [519] shows 
that it is possible to fit the Z-pole, NuTeV, and lepton universality data while satisfying 
this condition. 

The MEG (Mu-E-Gamma) experiment at PSI [520] plans to improve upon the MEGA 
limit by about two orders of magnitude. The MECO (Muon on Electron Conversion) 
experiment at Brookhaven [521] aims to improve the limits on fi — e conversion in nuclei by 
three orders of magnitude. Further constraints can be obtained from muon g — 2 [522,523], 
and the violation of CKM unitarity [524-526]. 

The NuTeV anomaly, even if it does not ultimately endure sustained scrutiny, stirs 
us to look past orthodoxies in our model-building (seesaw, SUSY, GUTs,...) and to 
ask broadly what is permitted by the data. The neutrino mixing solution is relatively 
conservative in its use of the neutrino sector to address the NuTeV question. Nonetheless, 
it makes interesting predictions about new particles at LHC, can be probed by a wide 
range of neutrino oscillation experiments, precision measurements and rare decay searches, 
and introduces an alternative to the seesaw paradigm. Whether this or another solution 
resolves the NuTeV anomaly, the NuTeV result serves to focus the imagination of the 
theorist on the opportunities presented by the experiments. 
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15 Conclusions 



In this report, we have presented a brief review of the present^ knowledge of neutrino 
physics and what we can learn from the planned experiments in the next decade. Three 
very important measurements that are guaranteed to have a significant impact on the 
search for physics beyond the Standard Model are: (i) the rate of /3/?o,/, which will inform 
us not only whether the neutrino is a Major ana or Dirac particle but may also provide 
information about the neutrino masses; (ii) the value of ^13, which will considerably narrow 
the field of flavor models and (iii) the sign of the Am^g, which determines the neutrino 
mass hierarchy and will also help guide our understanding of flavor physics. Within the 
three neutrino picture, more precise measiuements of the solar and atmospheric mixing 
angles will be helpful in discriminating among various new physics possibilities. 

Important though somewhat model-dependent constraints can be drawn from exper- 
imental searches for charged lepton flavor violating processes, such as // — > 67 or // — > e 
conversion in nuclei, and from searches for nonzero electric dipole moments of leptons. 
Keep in mind that the matter-antimatter symmetry of the Universe may have its expla- 
nation in the very same mechanism that generates the small neutrino masses, and that 
we may be able to test this hypothesis with enough low-energy information. 

Beyond the three neutrino picture, a very important issue is the status of the LSND 
result and whether the existence of light sterile neutrinos can be inferred from terrestrial 
neutrino oscillations experiments. The results of MiniBooNE, assuming they conflrms 
those from LSND, have the potential to revolutionize our current understanding of neu- 
trinos. Even if MiniBooNE does not confirm the LSND result, sterile neutrino effects can 
still be present in different channels at a more subdominant level, as has been suggested 
in several theoretical models. 

Another important issue in neutrino physics is the magnetic moment of the neutrino, 
which is expected at to be nonzero but very small within the standard picture of eV sized 
neutrino masses and in the absence of new physics at the TeV scale. Thus, evidence for 
a nonzero neutrino magnetic moment close to the current astrophysical limit of 10~^^ /ib 
would have to be interpreted as evidence of TeV scale new physics such as TeV scale left- 
right models, horizontal models, or large extra dimensions. Other unique probes of TeV 
scale physics are provided by neutrino oscillation experiments, thanks to their sensitivity 
to non-standard neutrino interactions. 

Finally, one can use results in neutrino physics to test the limits of the assumptions 
on which the Standard Model is based, including Lorentz and CPT invariance. 
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